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5.0 GLOBAL STRUCTURE OF ADD CODE 


This section of the manual is intended for the special user who wishes to 
modify the ADD code or adopt the ADD code to a diffferent computer. Section 5.0 
provides an overview of the code in terms of the principal tasks. These principal 
tasks are clearly labeled in the main program ALTMi'J and agree with the tasks listed 
in the Global Task Chart Section 5.2 and the Global Tree Structure Section 5.3. 
Internal flags, as opposed to input options, are described in Section 5.4. These 
flags are set by the code and control the calculation flow between different sub- 
routines. Flags used only within a subroutine are described in the detailed sub- 
routine descriptions in Section 7.0. 

Special problems associated with machine specific language are also described 
in this Section. The operation of a general I/O routine, which uses UNIVAC library 
I/O routine NTRAN, and a description of the data files is given in Section 5.5. This 
general I/O routine has been developed to allow NTRAN compatibility with ANSI stan- 
dard DEFINE FILE for easy conversion. The ADD code also uses a standard spline 
smoothing routine ICSVKU provided by International Mathematical and Statistical 
Libraries, Inc. The use of this subroutine is described in Section 5.6. All 
labeled C0MM0N block and DIMENSI0N statements are set by INCLUDE statements. Block 
sizes and EQUIVALENCE statements are set by PARAMETER statements. The use of these 
in the code is described in Section 5.7. 
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5.1 Main Program ALTMN 


Object 

Main program for ADD code 
Options 

See Input Options Section 4.2 
Variables 

See C0MM0N BL0CKS in Sect. 6.1. 

Theory 

The main program executes the major tasks according to the input options I0PT 
as shown by the global task chart described in Section 5.2. 
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5.2 Global Task Chart 


The major tasks are called in sequence by the main program ALTMN which is the 
only purpose of the main program. These tasks are executed, according to the input 
options, by calling only the subroutine listed on the Global Task Chart. The major 
tasks are labeled on the Task Chart as they are labeled by comment cards in the main 
program. Each major task may be executed independently and therefore may be over- 
layed. However, since it is not necessary to overlay every task, more than one sub- 
routine (task) may be in the same UNIVAC MAP segment of the overlay. The UNIVAC MAP 
segment is given in the brackets following the calling subroutine listed on the 
Global Task Chart. 


11-3 



5.2 Global Task Chart 




Task 

Subroutine (segment) 

Print Title Page 

0UTPUT (INITO) 

Read Input Variables 

REDINP (R00TM) 

Print Input Variables 

ECOINP (INITO) 





Set Up Problem 
Calculate Duct Geometry 
Normalize Geometric Variables 

Calculate, Store, Read Coordinates 
Set Up Blade Geometry 
Construct Slots in Duct 
Check Inlet Flow Conditions 
Calculate Inlet Flow Conditions 
Normalize Flow Variables 


GDUCT (INIT) 

FN0RM (INITO) 
SL0TA (C0RE) 

C00RST (C0RE) 

BLDG0M (F0RCESEG) 

SLOTA (C0RE) 

CKINPT (INITl) 

WFITER (INITO) 

FN0RM (INITO) 


Print Flow Set Up Conditions 
Print Default Input Parameters 
Print Calculated Coordinates 

WRTINP 

WRTGDC 

(R00TM) 

(R00TM) 

Calculate Inviscid Flow Field 
Calculate Potential Flow Field 

P0IS0N 

(PSD 

Calculate Rotational Flow Field 

CALINV 

(INITl) 

Calculate Blade Forces 

FORCE 

(F0RCESEG) 

Calculate Wall Bleed 

WBLEED 

(INITl) 
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5.2 GLOBAL TASK CHART (Cont'd) 


No. 

Task 

Subroutine (segment) 




7.0 

Calculate Viscous Flow Field 

S0LVI (C0MP) 

8.0 

Print Output Summary 


8.1 

Print Neglected Terms 

PRT0UT (NGADD) 

8.2 

Print Slot Flow Conditions 

SL0TA (C0RE) 

8.3 

Print Mean Flow Summary 

WRTSUM (WRTP) 

8.4 

Print Calculated Blade Forces 

F0RCE (F0RCESEG) 

8.5 

Print Boundary Layer Parameters 

BLPARM (ENDIN) 

8.6 

Print Solution on Output Line 

F0UTP (S0UTP) 


II-5 











5.3 Global Tree Structure By Task 


An overall view of the ADD code is shown in the Global Tree Structure shown on 
the following pages. This tree structure is arranged by tasks and labeled as they 
appear in the main program ALTMN. Column 1 shows the main program which calls the 
subroutines in column 2. Each subroutine in column 2 executes a major task and 
calls the subroutines in column 3. For any given calculation not all branches of 
the tree are called. As an example, if the coordinates have been calculated 
and stored on a data file (I0PT9=3) , then C00RST (Task 2.3) is not executed. C00R 
is a general subroutine which reads and interpolates the data on the coordinate 
file and BLKRED is a general subroutine for reading data blocks. Thus these two 
subroutines are used throughout the code. Subroutine C00R can also call C00R1 which 
is an approximate coordinate calculation used when I0PT9=O. The user who wishes to 
replace the UNIVAC general I/O routine NTRAN with standard ANSI FORTRAN DEFINE FILE 
need only modify subroutine BLKRED. The IMSL routines are only called by subroutine 
SM00TH and the user may replace the IMSL routines by only modifying this routine. 
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5.3 Global Tree Structure By Task 


0.0 

Print Title Page 



ALTMN 

OUTPUT 


1.0 

Read Input Variables 



ALTMN 

REDINP 

EC0INP 

SL0TA 



2.0 

Set Up Problem 





2.1 Calculate Duct Geometry 



ALTMN 

GDUCT 

SM00TH 

ICSYKU 

ICSFKY 

ICSFKU 

UERTST 



INITQ 

BLKRED 

BLKRED 


2.2 Normalize Geometric Variables 

ALTMN 

FN0RM 

TRBLDE 

AMU 




SL0TA 






Z.3 Calculate, 

Store, Read Coordinates 


ALTMN 

C00RST 

TPRINT 

C00R1 

XT 

XH 

CDS 

R0BRTS 

DR0BRT 

GBLADE 

FTHIK 



C00R4 

SM00TH 

CPLXl 

ICSVKU 




FINTG 

FCPLX 



C00R2 

C00R3 

CPLXl 

FETA 

FINTG 

FCPLX 




GBLADE 

C00R5 

FTHIK 
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5.3 Global Tree Structure By Task (Cont'd) 



2.4 

Set Up Blade Geometry 


ALTMN 

BLDG0M 

TRBLDE 

ALINE 

RDL0R 



FLINE 

CROSSl 

BLINE 

CLINE 

DLINE 

CR0SS2 

ELINE 

CR0SS1 


i 

SLETE 

TRLETE 

TRBLDE 

ALINE 

RDC0R 


2.5 

Construct Slots 

In Duct 

1 — 1 

1 


ALTMN SL0TA 



2.6 

Check Inlet Flow Conditions 


ALTMN 

CKINPT 

1 

TPRINT 

C00R 

C00R1 

XT 

XH 

CDS 

R0BRTS 

DR0BRT 

GBLADE 




BLKRED 




! 

WKTCKI 

GBLADE 

FTHIK 

2.7 

Calculate Inlet 

Flow Conditions 


ALTMN 

WRITER 

C00R 

BLKRED 

GBLADE 

FTHIK 




C00R1 

! 

XT 

XH 

CDS 

R0BRTS 

DR0BRT 

GBLADE 



FL0WIN 

UBLAS 

AMU 

CFC0LE 

BPLUSR 

FC0LES 

AMF 

BPLUSR 







II-8 
















5.3 Global Tree Structure By Task (Cont'd) 


2.8 Normalize Flow Variables 


ALTMN 

FN0RM 

TRBLDE 

AMU 



3.0 Print 

Flow Set Up Conditions 



3.1 

Print Input Parameters 



ALTMN 

WRTINP 




3.2 

Print Calculated 

Coordinates 



ALTMN 

WRTGDC 

C00R 

BLKRED 

GBLADE 

FTHIK 




C00R1 

XT 

XH 

CDS 

R0BRTS 

DR0BRT 

GBLADE 

4.0 Calculate Inviscid Flow Field 



4.1 

Calculate Potential Flow Field 



ALTMN 

P0IS0N 

TPRINT 

INITQ 

INTFRE 

BLKRED 




P01SCF 

TPRINT 

BLKRED 

FETA 

BLKWRT 




P0IS 

READPF 

BLKRED 




WRITPF 

BLKWRT 



SCURVA 

READPG 

BLKRED 




BLKRED 

BLKWRT 

QINTER 

BLKRED 





BLKWRT 
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5.3 Global Tree Structure By Task (Cont'd) 


A. 2 Calculate Rotational Flow Field 


ALTMN 

CALINV 

TPRINT 

C00R 

BLKRED 

GBLADE 

C00R1 



BLKWRT 

WRTCAL 



5.0 Calculate Blade Force 


6.0 Calculate Wall Bleed 


WBLEED 



BLKRED 

WRTBLD 



FTHIK 

XT 

XH 

CDS 

R0BRTS 

DR0BRT 

GBLADE 



7.0 Calculate Viscous Flow Field 


ALTMN 


S0LVI 


TPRINT 

BLKRED 

C00R 

F0RCL 

0RTFUN 

SLTFL0 


STRESI 

TURB 

MINVRT 

BLKWRT 

FC0RCT 

TURB2Q 

FAVER 

WRTS0V 


BLKRED 

GBLADE 

GBLADE 

AMU 

CFC0LE 

FC0LES 

TURB 

AMU 


AMU 
SUB LAY 

AMU 

BAMF0R 


BPLUSR 

AMU 
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5.3 Global Tree Structure By Task (Cont’d) 


8.0 

Print Output Summary 




8.1 

Print Neglected Terms 



ALTMN 


PRT0UT 

BLKRED 

C00R 

BLKRED 

GBLADE 

C00R1 

XT 

XH 

CDS 

R0BRTS 

DR0BRT 

GBLADE 




VARFUN 




8.2 

Print Slot Flow Conditions 



ALTMN 

SL0TA 





8.3 

Print Mean Flow Summary 



1 

ALTMN 

1 

1 

WRTSUM 

BLKRED 

QPTT-IAX 




8.4 

Print Calculated Blade Force 


ALTMN 


F0RCE 

GBLADE 

FTHIK 





BLKRED 

C00R 

BLKRED 

GBLADE 

FTHIK 





C00R 

XT 

XH 

CDS 

R0BRTS 

DR0BRT 

GBLADE 


8.5 

Print Boundary 

Layer Parameters 


1 

ALTMN 


BLPAEM 

BLKRED 

AMU 

C00R 

BLKRED 

GBLADE 

XH 





C00R1 

XT 

CDS 

R0BRTS 

DR0BRT 

. 



, 



GBLADE 
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5.3 Global Tree Structure By Task (Cont'd) 


■ 

8.6 Print Solution on Output Line 



ALTMN 

FOUTP 

INITFV 

TPRINT 

ALINE 

BLINE 

CR0SS1 




BLINE 

FLINE 

ALINE 

BLINE 

CR0SS1 




RDC0R 

BLINE 

CLINE 

DLINE 

CR0SS2 

ELINE 

BLKRED 



REC0RD 

C00R 

BLKRED 

GBLADE 

FTHIK 



BLKRED 

WRTF0U 








5.4 List of Internal Flags 


Flag Name 


Function 


N0PT1 = 0 
= 1 

N0PT2 = 1 
= 2 

N0PT3 = 1 
= 2 

N0PT4 = 1 
= 2 

N0PT5 = 0 
> 0 

N0PT6 

N0PT7 

N0PT8 

N0PT9 = 0 
= 1 

N0PT1O 
N0PT11 
N0PT12 
N0PT13 
N0PT14 
N0PT15 = 0 
= 1 

N0PT16 
N0PT17 
N0PT18 = 0 
= 1 

N0PT19 
N0PT2O = 0 
= 1 

N0PT21 = 0 
= 1 

N0PT22 


Set up initial flovi' 

Calculate flow station J 

Adiabatic wall 

Prescribed wall temperature 

Calculate and store approximate coordinates 

Calculate local approximate coordinates 

Duct with centerbody 

Duct without centerbody 

Continue calculation 

Fatal error - stop calculation 

Force backpressure iteration counter 

Block I/O counter in CALINV 

Block I/O counter in S0LVI 

Complete complex function calculation 

Shortened complex function calculation 

Not used 

Not used 

Not used 

FCPLX Flag 

Not used 

Integrate along S 
Integrate along n 
S0LVI station counter 
Not used 
Turbulent flow 
Laminar flow 

Weight flow iteration counter 
Optimize step size KDS 
Fixed KDS 

Greater than critical Reynolds number 
Less than critical Reynolds number 
Not used 
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5.5 General Input/Output for Data Files 


The file assignments are described in Table I. This table shows the unit number, 
the names of the arrays stored on the file, the block (record) length in words, 
the number of blocks (records), and the subroutine generating the file. All unit 
numbers are set by parameter statements (in brackets). Thus, changing the unit 
numbers requires changing only the paramater statements. All arrays are single 
precision except AAF or unit LDRUM which is double precision. The data stored on 
unit JDRUM and unit KP0IS are similiar. Unit JDRUli has coordinate data for a non- 
uniform mesh used by subroutine S0LVI to obtain the solution of the viscous flow 
field. Unit KP0IS has the coordinate data for a uniform mesh used by subroutine 
P0IS to obtain the compressible axisymmetric potential flow solution. The number 
of blocks (records) stored on the files is set primarily by the number of stream- 
wise coordinate stations calculated by subroutine C00RST. The maximum numbers of 
streamwise stations (IS) is set by a parameter statement such that JL _< IS. The 
number of words in a block (record) is set primarily by the number of streamlines 
KL calculated by subroutine C00RST. The maximum number of streamlines 1ST is set 
by a parameter statement such that KL _< 1ST. Subroutine P0IS has provision for 
read/write of NST blocks (records) where NST is set by a parameter statement. 

All data files are read/written by a general I/O routine BLKRED, with a second 
entry point BLKWRT, which uses a UNIVAC Library I/O routine NTRAN. This subroutine, 
BLKRED, has been specially developed to allow NTRAN compatability with ANSI standard 
F0RTRAN DEFINE FILE. Thus, the user who wishes to use DEFINE FILE need change only 
one subroutine, BLKRED, which is described in detail in Section 7. In addition to 
the general routine BLKRED, special I/O routines have been developed which still use 
BLKRED. Subroutine C00R reads two successive blocks of coordinate data and inter- 
polates to obtain the fine grid used by subroutine S0LVI when KDS > 1. Subroutines 
READPF, READPG, WRITPF, which also use BLKRED, is specially developed to read NST 
blocks (records). The logic for counting blocks (records) is built into these sub- 
routines and described in Section 7. 

When using the global iteration option I0PT14 > 0 or for post processing of data 
such as drawing CALC0MP plots, the corresponding record numbers for the inviscid 
solution stored on file MDRUM, the viscous solution stored on file NDRUM, and the 
coordinate data stored on file JDRUM must be obtained. This is shown in Table II. 

The parameters controlling the record counting are shown on the top of the table. 

The mesh size parameters JLXKL is stored on the coordinate file along with scale 
factors. The remaining parameters controlling the record counting are stored on 
the first record on file NDRUli. For I0PT14 = 0, the viscous solution is stored 
only for each Jth coordinate station regardless of the value of KDS used in the 
calculation. Therefore, set KDS = JKDS = 1. For I0PT14 > 0, the viscous solution 
is stored for every interpolated station on the fine grid (JKDS = 1, KDS). There- 
fore, to recall data at a Jth coordinate station set JKDS = 1 and KDS equal to the 
value used in viscous flow calculation. In some cases the solution stops before 
completion. N0PT7 and N0PT8 are the number of records stored on files MDRUM and 
NDRUM respectively. Therefore, 

N0PT7 j< I0PT16 - I0PT15 - 1 

N0PT8 _< KDS (I0PT16 - I0PT15) + KDS - 1 
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Table I 


Table of File Assignments 


UNIT 

NO. 

ARRAY 

NAME 

BLOCK LENGTH 
(WDS) 

NO. 

BLOCKS 

SUBROUTINE 
WRITING BLOCK 

*8(NDRUM) 

F(NEQ,3,IST) 

FPARM(15) 

3915 

JL-2 

S0LVI 

9(JDRUM) 

JSTEP 

Q(19,IST) 

RHS(IO) 

RMS(IO) 

RTS (10) 

DSTEP 

QPARM(9) 

2514 

JL 

C00RST 

lO(CDRUM) 

FF(17,2,IST) 

4420 

1 

F0RCE 

ll(LDRUM) 

AFF(LNGTO)** 

6600 

5 

S0LVI 

12 (LF0RC) 

F0RC 

780 

1 

F0RCE 

19(KP0IS) 

JSTEP 

Q(19,1ST) 

RHS(IO) 

RMS (10) 

RTS (10) 

DSTEP 

QPARM(9) 

2514 

JL 

C00RST 

22 (MDRUM) 

FIV(NEQ,3,IST) 

FIPARM(15) 

3915 

JL-2 

CALINV 

23(NFDRM) 

F(1ST4) 

134 

JL+NST-2 

P0ISCF 

24(NPDRM) 

P(IST) 

130 

JL+NST-2 

P0IS 

25 (NGDRM) 

BLK(NGIST) 

400 

JL+NST-2 

P0ISCF 


* KDS*JL-2 blocks when lOPTlA > 0, JL _< IS 
**Double precision variables 
IS = 100 
NST = 25 
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Table II 


JL 

KL 

*KDS 

N0PT7 

N0PT8 

I0PT15 

I0PT16 

J 

JRECFI 

JRECFV 

**JKDS 


JRECFI 

JRECFI 

JRECFI 


JRECFV 

JRECFV 


* KDS = 
**JKDS ■ 


Table of Record Counting 


= QPARM(3) 

= QPARM(4) 

= FPARM(IO) 
= FPARM(12) 
= FPARM(13) 
= FPARM(14) 
= FPARM(15) 


, No. streamwise coordinate, JL j< IS 
, No. streamline coordinate, KL _< 1ST 
, No. steps/station fine grid 
, No. MDRUM records 
, No. NDRUM records 
, First J station 
, Last J station 

, Record no. for coordinates, J = 1, JL 
, Record no. for inviscid solution 
, Record no. for viscous solution 
, Fine grid station counter, JKDS = 1, KDS 


FIV(I,N,K) Array 

J-I0PT15-1 , N = 3, for J > I0PT15 + 1 
1 , N = 2, for J = I0PT15 + 1 
1 , N = 1, for J = I0PT15 


F(I,N,K) Array 

KDS*(J-I0PT15)+JKDS-1 , N = 3, J > I0PT15 

1 , N = 2, J = I0PT15 


1 for I0PT14 = 0 
1 for Jth record 


11-16 



5.6 IMSL Routines 


The ADD code uses a standard least squares cubic spline routine provided by 
International Mathematical and Statistical Libraries, Inc. This routine is called 
ICSVKU and is described in the IMSL library dated November 1975. Subroutine ICSVKU 
computes a least squares approximation to a given set of points by a cubic spline 
with a given number of knots. Since the successful use of splines for purposes of 
providing a smooth approximation to a given set of points depends strongly on the 
proper placement of the knots, ICSVKU starts with a given set of knots and varies 
them one by one in order to determine the knot locations that minimizes the least 
square error. Subroutine ICSVKU also uses subroutines ICSFKU, ICSFKV and VERTST. 

The least square spline routine ICSVKU is called by only subroutine SM00TH 
which is described in Section 7. Subroutine SM00TH, however, uses ICSVKU to spline 
fit the second derivative of the data points. Therefore, for the user who does not 
have the IMSL library or wishes to replace ICSVU, a brief description of the sub- 
routine is given on the following page. 
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Subroutine ICSVKU (Arg. List) 


Object 

Calculate least squares spline fit 
Options 


None 

Arg. List 


X(I) 


F(I) 

F 

I 

NX 

N 



XK(K) 

\ 

NXK 

M 

Y(K) 

Y 

K 

C(K,3) 

"kl 

1C 

NXK-1 

ERR0R 


WK 


lER 



Vector containing NX abscissa points 
Vector containing NX ordinate points 
Number of elements in X and F 
Vector containing NXK knot locations 
Number of knots 

Vector of NXK-1 spline coefficients 
Matrix of NXK-1 x 3 spline coefficients 
Dimension of spline coefficient row 
Cubic spline error 
Work area of dimensions NS=(NXF+6) 

Error parameter 


Theory 

The vectors (X^, 1 = 1» NX) containing the points to be fitted with a cubic 

spline are input along with the vector (Xj^, K = 1, NXK) which is the initial guess 
for the knot locations. The subroutine returns the spline coefficients (Y^., C^^ 

L = 1, 3, k = 1, NXK-1) and the new locations of the knots K = 1, NXK) where 

the spline coefficients are given in the form 

F(T) = ((Cj^ 3+D + Cj^2> + D + Cj^^) + D + \ (D 

D = T-% (2) 
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where 




< T < 


'\y 

X 


K+1 


(3) 


For this subroutine the end points are fixed and only the interior knots are moved. 
Subroutine SM00TH treats Eq. (1) as an equation for the second derivative. There- 
fore, Eq. (1) is integrated analytically to obtain the smoothed curve with continuous 
fifth derivatives. 

Reference 


deBoor, Carl and John R. Rice: Cubic Spline Approximation II - Variable Knots. 

Computer Sciences Department TR21 Purdue University, April 1968. 
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6.0 GLOBAL VARIABLES 


6.1 List of Labeled COMMON BLOCK Variables 


Name 

Object 

AC0NS 

Input blade geometry variables 

AC0NSO 

Blade geometry parameters 

ADPS 

Stored radial coordinate for slot interpolation 

AINV 

Variables for CALINV I/O 

AKK 

Complex coordinate flags 

AMATRX 

Variables for viscous flow matrix equations 

AVER2 

Variables for mass flow average calculation 

BCPLX 

Variables for complex coordinate calculation 

BSL0T 

Variable for slot geometry 

BTHIK 

Blade thickness variables 

CCPLX 

Variables for complex coordinate calculation 

CINF 

Parameters for Poisson equation 

C0NST 

Calculated gas properties 

C0r.Z 

Coordinate Functions 

C0RE2 

Calculated wall variables 

CSL0T 

Dimensionless slot geometry variables 

DERIV 

Strut force functions 

DREDl 

Coordinate variables 1/0 

DRED2 

Coordinate variables I/O 

DSKB 

Parameters for I/O routine 

DSL0T 

Parameters for slot geometry 

DUC0UT 

Calculated wall coordinates 
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6.1 List of Labeled C0MMON BL0CK Variables (Cont’d) 


DUCTIN 

FC0R 

Five 

FLAGS 

FL0WI 

F0RS 

F0RS2 

FUNG 

INTINP 

MATRXA 

NGTERM 

REALIN 

SPECFD 

SPCGD 

SPI0 

SIRES 

STRMES 

SVARB 

TITLIN 

TRBL 

TURBS 

TURB2 


Input duct and flow data 
Mass Flow weighted flow errors 
Inviscid flow variables I/O 
Internal flags 
Initial flow variables 
Calculated blade force variables 
Blade force variables 
Viscous flow variables 

Mesh parameters, input options, debug options 
Variables for viscous flow matrix equations 
Variables for neglected terms calculation 
Input flow properties, gas properties 
Variables in Poisson equation 

Variables for streamline curvature calculation 
Initial flow variables, mass average variables 
Geometric functions 

Poisson solver stretching parameters 
Mesh variables, flow variables 
Title card 

Variables for two equation turbulence model 
Turbulent viscosity, conductivity 
Parameters for automatic step size algorithm 
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6.1 List of Labeled C0MM0N BL0CK Variables (Cont'd) 

WLINE Variables for locating arbitrary line 

XF0U Parameters for output data line calculation 

ZF0U Variables for output data line calculation 
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6.2 List of C0MM0N Block Variables 

LIST OF VARIABLES IN C01>IM0N/AC0NS/ 
Input Blade Geometry Variables 


C0NSTI(1,L) 

= 

y 

9 

Independent variable (stacking line) (ft) 

C0NSTI(2,L) 

= 

a 

s 

9 

Stagger angle from blade face (deg.) 

C0NSTI(3,L) 

= 

c 


Blade chord (ft) 

C0NSTI(4,L) 

= 

t/c 


Blade thickness/chord 

C0NSTI(5,L) 

= 


9 

Circular arc camber (deg.) 

C0NST1(6,L) 

= 

1 

X 


Dependent variable (ft) 

ISHAPE 



9 

Blade shape index 

NB 

= 

Nb 


Number of blades 

NUM 

= 

N 


Number of input points 

0MEGZI 

= 

fiz 

9 

Rotor speed (rpm) 

RCLOI 

= 

’^CLO 


Radial location of stacking line (ft) 

THCLI 

= 

®CLO 


Rotation of stacking line (deg.) 

ZCLOI 

= 

^CLO 


Axial location of stacking line (ft) 

L 

= 

l.N 
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LIST OF VARIABLES IN C0MM0N/AC0NS0/ 

A 

Blade Geometry Parameters 


C0NST(1,L) 

= 

y’/rr . 

Independent variable 

,CpNST(2,L) 

= 

“s 

Stagger angle to axis (deg) 

C0NST(3,L) 

= 

c/r^ , 

Chord 

C0NST(4,L) 

= 

t/rr . 

Blade maximum thickness 

C0NST(5,L) 

= 

<^c 

Circular arc camber (deg) 

C0NST(6,L) 

= 

x'/rj. , 

Dependent variable 

C0NST(7,L) 

= 


Streamline distance to center line (CL) 

C0NST(8,L) 

= 

s , 

Streamwise coordinate of center line (CL) 

C0NST(9,L) 

= 

n , 

Normal coordinate of center line (CL) 

ALPHS 

= 

CXs(n) . 

Stagger angle to axis (deg) 

ALPSH 

= 

“SH 

Stagger angle to axis hub (deg) 

ALPST 

= 

“sT ’ 

Stagger angle to axis tip (deg) 

CH0RD 

= 

c(n)/rj. , 

Chord 

CH0RDH 

= 


Chord hub 

CH0RDT 


Ci/rr , 

Chord tip 

GAP 

= 

g/rr . 

Gap 

QfMEGZ 

= 


Rotor speed 

PHIC 

= 

‘f'c » 

Circular arc camber angle (deg) 

RCLH 

= 

’^an^^r ’ 

Radius CL hub 

RCLT 

= 

’^CLT^’^r » 

Radius CL tip 

RCL(? 

= 

’^CLO^’^r ’ 

Radius to stacking line origin 

RLEH 

= 

^LEH^^r ’ 

Radius leading edge hub 

RLET 

= 

^LET^’^r ’ 

Radius leading edge tip 
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LIST OF VARIABLES IN C0MI40N/AC0NS0/ (Cont'd) 


S(5ld 

— 

a 

THCL 

= 

®CL 

THIKM 

= 

t(n)/c(n) 

THIKN 

= 

t(n)/r^ 

ZCLH 

= 

^CLH'^^r 

ZCLT 

= 

‘^CLT'^’^r 

ZCL0 

= 


ZLEH 

= 

^LER/’^r 

ZLET 

= 

^LET^’^r 


Blade solidity 

Rotation of stacking line, deg 

Thickness/ chord 

Blade maximum thickness 

Axial location CL hub 

Axial location CL tip 

Axial location stacking line origin 

Axial location leading edge hub 

Axial location leading edge tip 


* 


See Subroutine GBLADE 
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LIST OF VARIABLES C0MM0N/ADPS/ 

Store Radial Coordinate for Slot Interpolation 
DPSI(K) = R (K) , Radial coordinate (dimensionless) 

K = 1,KL 
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LIST OF VARIABLES IN C0MM0N/AKK/ 


Complex Coordinate Flags 


Ml 

y 

Flags for debug printout 

M2 

» 

See subroutine FCPLX 


11-27 



LIST OF VARIABLES IN C0MM0N/AMATRX/ 
Variables for Viscous Flow Idatrix Equations 


AD(I,J) 

= 

■»ij 

> 

Elements of D matrix (see S0LVI) 

ADI(I,J) 

= 



Elements of matrix (see S0LVI) 
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LIST OF VARIABLES IN C0MM0N/AINV/ 


CINP(1,K) 

CINP(2,K) 

CINP(3,K) 

CINP(4,K) 

K 


n 

n. 


1,KL 


Variables for CALINV I/O 

, Static pressure (P/Pj.) 

, Total pressure (Pq/P^.) 

, Swirl angle (deg) 

, Total temperature (T^/Tj.) 
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LIST OF VARIABLES IN C0I>1M0N/AVERZ/ 


ASH 

AST 

DASHl 

DASTl 

DENTPl 

DPSIl 

DQSHl 

DQSTl 

DTHEOl 

ENTP 

ENTPO 

QSH 

QST 

THEO 

THEOO 


AA^ 


AI 

A'l' 

AQt 

A0 

o 


I 

I 

o 

Qh 

Qt 

®o 

0 

oo 


Surface area hub (dimensionless) 

Surface area tip (dimensionless) 

Increment wall area hub (dimensionless) 
Increment wall area tip (dimensionless) 
Increment in Entropy (dimensionless) 

Increment in mass flow (dimensionless) 
Increment wall heat flux hub (dimensionless) 
Increment wall heat flux tip (dimensionless) 
Increment in total temperature (dimensionless) 
Entropy (dimensionless) 

Inlet entropy (dimensionless) 

Total heat flux hub (dimensionless) 

Total heat flux tip (dimensionless) 

Total temperature (dimensionless) 

Inlet total temperature (dimensionless) 
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LIST OF SYMBOLS IN C0MM0N/BCPLX/ 
Variables for Complex Coordinate Calculation 


A(1,I) 

= 

A. 

1 

, Source strength (dimensionless) 

A(2.I) 

= 


, Location of pole (dimensionless) 

A(3,I) 

= 


, Wall angle change (deg) 

A(4,I) 

= 

r . 
1 

, Radius in z plane (dimensionless) 

A(5,l) 

= 

^i 

, Angle in z plane (radians) 

A( 6 ,I) 

= 

^i 

, Relative X distance in z plane 

(dimensionless) 

A(7,I) 

= 

Y 

1 

, Relative Y distance in z plane 

(dimensionless) 

B(l.I.K) 

= 

ASX 

s 

, Change in coordinate X (dimensionless) 

B(2,J,K) 

= 

ASYg 

, Change in coordinate Y (dimensionless) 

B(3,I,K) 

= 

ASC 3 

, Change in coordinate C (dimensionless) 

B(4,I,K) 

= 

ASr; 

s 

, Change in coordinate p (dimensionless) 

X(1,K) 

= 

S 

, Streamwise coordinate (dimensionless) 

X(2,K) 

= 

n 

, Normal coordinate (dimensionless) 

X(3,K) 

= 

X 

, X coordinate in z plane (dimensionless) 

X(4,K) 

= 

Y 

, Y coordinate in z plane (dimensionless) 

X(5,K) 

= 

K 

, coordinate in w plane (dimensionless) 

X( 6 ,K) 

= 

n 

, coordinate in w plane (dimensionless) 

X(7,K) 

= 


, Streamwise derivative of ? (dimensionless) 

X( 8 ,K) 

= 


, Streamwise derivative of n (dimensionless) 
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LIST OF SYMBOLS IN C0MM0N/BCPLX/ 


X(9,K) 

= 

X(S+AS) 

> 

X(10,K) 

= 

Y(S+AS) 

9 

X(11,K) 


5(S+AS) 

, Coordinates at Station S+AS 

X(12,K) 

= 

ti(S+AS) 

9 

X(13,K) 

= 

V 

, Metric scale coefficients (dimensionless) 

X(14,K) 

= 

^ss 

, Second derivative of 5 (dimensionless) 

X(15,K) 

= 

^sn 

, Cross derivative of ? (dimensionless) 

X(16,K) 

= 

V 

n 

, Normal derivative of V (dimensionless) 

X(17,K) 

s 

V 

s 

, Streamwise derivative of V (dimensionless) 


11-32 



LIST OF VARIABLES IN C0MM0N/BLS0T/ 


ALSLTKN.L) 

HSL0TI(N,L) 

P0SLTI(N,L) 

T0SLTI(N,L) 

W0SLTI(N,L) 

ZSL0TI(N,L) 


Variables for Slot Geometry 


* 

“SL 

* 

Slot swirl angle (deg) 

= 

"^SL 


Slot height (ft) 

= 

^OSL 

y 

Slot total pressure (psfa) 

s 

’^OSL 

> 

Slot total temperature (deg R) 

s 

^OSL 

f 

Slot weight flow (Ib/sec) 

= 

^SL 

f 

Slot axial location (ft) 


1=1 Tip wall 
=2 Hub wall 
-=1,15 Slot Number 
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VARIABLES IN C0MM0N/BTHIK/ 


KBLADE 



, Number of points 

XK(I) 

= 

^I 

, Fractional chordwise distance 

YK(I) 

= 

^I 

, Thickness/Chord 


1=1, KBLADE 
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LIST OF VARIABLES IN C01M0N/CCPLX/ 


Variables for Complex Coordinate Calculation 


NPTS 


N 

P 

9 

Number of singularities in complex 
transformation 

NS0URC 

= 

Ns 

9 

Number of sources in z plane 

^RDERl 

= 

0l 

9 

Absolute magnitude of largest term 

0RDER2 

= 

°2 

9 

Absolute magnitude of largest term 

0RDER3 

= 

°3 

9 

Absolute magnitude of largest term 

SLO 

= 

^LO 

9 

Length of streamwise coordinate 

XDS 

= 

dS 

9 

Step size for complex integration 

XDN 

= 

dS 

9 

Step size for complex integration 
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LIST OF VARIABLES IN C0MM0N/CINF/ 


Parameters for Poisson Equation 


AMINF 

m 


Freestream 

Mach number 

AN(K) 

= 


Transverse 

coordinate 

DEDN(K) 

= 

(dn/dn)j^ 

Transverse 

coordinate stretching 

D2EDN(K) 

= ■ 

(d^n/dn^)j^ 

Transverse 

coordinate stretching 

PINF 

= 

^OO 

Freestream 

static pressure (dimensionless) 

PSIKL 

= 

'p 

^00 

Freestream 

stream function (dimensionless) 

PTINF 

= 

^00 

Freestream 

total pressure (dimensionless) 

RH0INF 

= 

Poo 

Freestream 

density (dimensionless) 

R0TINF 

= 

Pqoo 

Freestream 

total density (dimensionless) 

TEMINF 

= 

®co 

Freestream 

static temperature (dimensionless) 

TTINF 

= 

®Ooo 

Freestream 

total temperature (dimensionless) 

UINF 

= 


Freestream 

velocity (dimensionless) 

UPINF 

= 

^(()00 

Freestream 

tangential velocity (dimensionless) 

USINF 

= 

Us» 

Freestream 

streamwise velocity (dimensionless) 

UUO 

= 

Uo 

Reference velocity (dimensionless) 

WO 

= 

Vo 

Reference metric coefficient 
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LIST OF VARIABLES IN C0MM0N/C0NST/ 


Gas Properties 


ACHI 

— 

X 

, Clauser constant (0.016) 

AKAPPA 

= 

K 

, von Karman constant (0.41) 

APLUS 

= 


, van Driest constant (26.0) 

CPR 

= 

s 

, Specific heat at constant pressure 

(5997.0 ft2/sec2/deg R) 

CVR 

= 

C 

V 

, Specific heat at constant volume 

(3283.0 ft^/sec^/deg R) 

EP 

= 

e' 

, 2.7182818 (exponential) 

GAMMA 

= 

Y 

, Ratio of specific heat (1.4) 

GASR 

= 

R 

9 2 

, Gas constant (1714.0 ft^/sec /deg R) 

PI 

= 

IT 

, 3.1415926 

PRESR 

= 

P 

r 

, Reference static pressure (psfa) 

PRL 

= 

P^^L 

, Prandtl number laminar 0.70 

PRT 

= 

Pr^ 

, Prandtl number turbulent 0.72 

RH0R 

= 

Pr 

, Reference density (slugs/ft^) 

SNDR 

= 

^r 

, Reference speed of sound 1116.0 (ft/sec) 

TEMPR 

= 

T 

r 

, Reference temperature (deg rankin) 

TI 

= 


, (0.1745329 radians/deg) 

VISCR 

= 

Pr 

, Reference molecular viscosity 

(0.370 X 10"^) 
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LIST OF VARIABLES IN C0MM0N/C0RE/ 


Coordinate Functions 


Q(1,K) 

= 

R 

, Radius (dimensionless) 

Q(2,K) 

= 

Z 

, Axial distance (dimensionless) 

Q(3,K) 

= 

3R/3n 

, Normal derivative of radius (dimensionless) 

Q(A,K) 

= 

3R/3S 

, Streamwise derivative of radius 

(dimensionless) 

Q(5.K) 

= 

cos^O 

, Axisymmetric flow angle 

Q(6,K) 

= 

V 

, Metric scale coefficient (dimensionless) 

Q(7,K) 

= 

3V/3n 

, Curvature of streamline (dimensionless) 

Q(8,K) 

= 

3V/3S 

, Curvature of potential line (dimensionless) 

Q(9,K) 

= 

X 

, Distance along streamline (dimensionless) 

Q(10,K) 

= 

Y 

, Physical distance across duct (dimensionless) 

Q(11,K) 

= 

Y/Y^ 

, Fractional distance across duct (dimensionless) 

Q(12,K) 

= 

A 

, Area between adjacent streamlines 

(dimensionless) 

Q(13,K) 

= 

G 

, Gap between blade surfaces (dimensionless) 

Q(14,K) 


3G/3n 

, Normal derivative of blade surface 

(dimensionless) 

Q(15,K) 

= 

3G/3S 

, Streamwise derivative of blade surface, 

(dimensionless) 

Q(16,K) 

= 

3n/3n = X 

, Transform of normal coordinate (dimensionless) 

Q(17,K) 

= 


, Second derivative (dimensionless) 

Q(18,K) 

= 

n 

, Normal coordinate (dimensionless) 

Q(19,K) 

= 

T1 

, Transformed normal coordinate (dimensionless) 

K 


1,KL 

, Number of streamlines (dimensionless) 
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LIST OF VARIABLES IN C0MM0N/C0R2/ 


Caluclate Wall Values 


RHS(l) 

ss 


» 

Radius ID wall (dimensionless) 

RHS(2) 

= 

dRjj/dZ 

> 

Derivative (dimensionless) 

RHS(3) 

= 


9 

Input curvature at station S (dimensionless) 

RHS(4) 

= 


9 

Axial distance ID wall (dimensionless) 

RHS(5) 

= 


9 

Metric coefficient ID wall (dimensionless) 

RHS(6) 

= 

dVjj/dS 

9 

Derivative (dimensionless) 

RHS(7) 

= 

dVy/dn = -K 

9 

Derivative (dimensionless) 

RHS(8) 

= 


9 

Arc length ID wall (dimensionless) 

RHS(9) 

= 

O 

“h 

9 

Mass flow bleed (dimensionless) 

RHS(IO) 

= 


9 

Wall temperature (dimensionless) 

RMS(l) 

- 

^^I 

9 

Input radius ID wall (dimensionless) 

RMS(2) 

= 

^HI 

9 

Input axial distance ID wall (dimensionless) 

RMS(3) 

= 

\l 

9 

Input arc length ID wall (dimensionless) 

RMS (4) 

= 

-%I 

9 

Input wall curvature ID wall (dimensionless) 

RMS(5) 

= 

Rjl 

9 

Input wall radius OD wall (dimensionless) 

RMS(6) 

= 

Zti 

9 

Input wall axial length OD wall (dimensionless) 

RMS (7) 

= 

x^i 

9 

Input wall arc length OD wall (dimensionless) 

RMS (8) 

= 

-^TI 

9 

Input wall curvature OD wall (dimensionless) 

RMS(9) 





RMS(IO) 





RTS(l) 

= 


9 

Radius OD wall (dimensionless) 

RTS (2) 

= 

dRp/dZ 

9 

Derivative (dimensionless) 
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C0MM0N/C0R2/ (Cont'd) 


RTS (3) 

a: 

-ICpj(S) 

, Input curvature at station S (dimensionless) 

RTS (4) 

= 

^T 

, Axial Distance OD wall (dimensionless) 

RTS(5) 

= 


, Metric coefficient (dimensionless) 

RTS (6) 

= 

dV^/dS 

, Derivative (dimensionless) 

RTS (7) 

= 

dV^/dn = -K 

, Derivative (dimensionless) 

RTS (8) 

= 


, Arc length OD wall (dimensionless) 

RTS (9) 

= 

O 

nXj 

, Mass flow bleed (dimensionless) 

RTS (10) 

= 

0 

T 

, Wall temperature (dimensionless) 
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LIST OF VARIABLES IN C0MM0N/CSL0T/ 


ALSL0T(N,2) = 
HSL0T(N,L) 
JLSL0T(N) 
LISZE(N,L) 

LSL0T(L,N) 
POSL0T(N,L) = 
TOSL0T(N,L) - 
ZSL0T(N,L) 


SL 


H 


SL 


0 

1 


OSL 


OSL 


^SL 


Dimensionless Slot Geometry Variables 

, Slot swirl angle (rad) 

, Slot height (dimensionless) 

, Number of slots 

, Small slot 

, Large slot 

, Slot index 

, Slot total pressure (dimensionless) 

, Slot total temperature (dimensionless) 

, Axial slot location (dimensionless) 


N = 1 Tip wall 
= 2 Hub wall 

L = 1,15 Slot index 
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LIST OF VARIABLES IN C0MMto/DERIV/ 


Strut Force Functions 


DF(1,K) 

= 

[H /XV]'^ 

® K-1/2 

9 

Streamwise 

strut force/volume (dimensionless) 

DF(2,K) 

= 


9 

Tangential 

strut force/volume (dimensionless) 

DF(3,K) 

= 


9 

Total pressure loss/volume (dimensionless) 

DF(4,K) 

= 

k-1/2 

9 

Coordinate 

distortion (dimensionless) 

DF(5,K) 

= 

[H /XV]-^"^ 
s K-1/2 

9 

Streamwise 

strut force/volume (dimensionless) 

DF(6,K) 

= 


9 

Tangential 

strut force/volume (dimensionless) 

DF(7,K) 

= 

[$„/XV]'^“^ 

B k-1/2 

9 

Total pressure loss/volume (dimensionless) 

DF(8,K) 

= 

[x]^-i 

k-1/2 

9 

Coordinate 

distortion (dimensionless) 


where 


X = dn/dn 



LIST OF VARIABLES IN C0MM0N/DSKB/ 


Parameters for 1/0 Routine 


DSKL0C (UNIT) 

= 

N 

where: 

UNIT 

= 

Unit 

number (See PARAMETER list) 

N 

s 

Last 

block number read 
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LIST OF VARIABLES C0MM0N/DRED1/ 
C0MM0N/DRED2/ 

Coordinate Variables I/O 

BL0CK(I) 1=1, ISL Coordinates at J+1 

BL0CK1(I) 1=1, ISL Coordinates at J 

See EQUIVALENCE statement. 

where for J+1 station we have 


JSTEP 

J+1 

9 

Station or block number 

Q2(I,J) E 

Q(I,J) 

9 

See C0MM0N/C0RE/ 

RHS2(I) E 

RHS(I) 

9 

See C0MM0N/C0R2/ 

RMS2(I) E 

RMS (I) 

9 

See C0MM0N/C0R2/ 

RTS2(I) E 

RTS (I) 

9 

See C0MM0N/C0R2 

DSTEP 

DS 

9 

Step Size in streamwise direction 

QPARM(l) = 

r 

r 

9 

Reference radius 


QPARM(2) 


QPARM(3) = JL , No. streamwise stations 

QPARM(4) = KL , No. streamlines 

QPARM(5) 

QPARM(6) 

QPARM(7) 

QPARM(8) 

QPARM(9) 

QPARM(IO) 
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LIST OF VARIABLES IN C0MM0.N/DUC0UT/ 
Calcualted Wall Coordinates 


R(1.1,J) 

— 


, Radius of hub (dimensionless) 

R(2,1.J) 

= 


, Radius of tip (dimensionless) 

R(1,2,J) 

= 

Sj(Zj) 

, Mass flow tip bleed (dimensionless) 

R(2.2,J) 

= 


, Mass flow of hub bleed (dimensionless) 

R(1,3,J) 

= 

w 

, Wall temperature of tip (dimensionless) 

R(2,3,J) 

= 

W 

, Axial distance of hub (dimensionless) 

R(1,3,J) 

= 

Zt 

, Axial distance tip (dimensionless) 

R(2.3,J) 

= 


, Axial distance hub (dimensionless) 
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LIST OF VARIABLES IN C0MM0N/DUCTIN/ 

Input Duct and Flow Data 

AINI(L) = AINPUT(I,J,K) , See C0MM0N/SPI0/ 

BINI(L) = BINPUT(I,J,K) , See C01>IM0N/SPI0/ 

I = 1,5 

J = 1,2 

K - 1,KLL 

L = 5*(K-1)+I+(J-1)*5*KLL 

DUCTI(I) I = 1,15 , Arbitrary duct geometry parameters 

RDII(L) = R(1,1,L)*RADR, Tip wall coordinates (ft) 

RD2I(L) = R(1,2,L)*RADR, Hub wall coordinates (ft) 

ZDIJ(L) = R(4,1,L)*RADR, Tip wall coordinates (ft) 

ZD2J(L) = R(4,2,L)*RADR, Hub wall coordinates (ft) 

L = 1, JL 
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LIST OF VARIABLES IN C0MM0N/FC0R/ 


Mass Flow Weighted Flow Errors 


EAP 

= 

^U()> 

, Error 

in swirl velocity 

EAS 

= 

^Us 

, Error 

in streamwise velocity 

EEN 

= 

^I 

, Error 

in Entropy 

EPO 

= 

^TTO 

, Error 

in total pressure 

ESI 

= 


, Error 

in mass flow 

ETO 

= 

P 

^TTO 

, Error 

in total pressure 
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LIST OF VARIABLES IN C0MM0N/FIVC/ 





Inviscid Flow Variables I/O 

FIV(1,L,K) 

= 

'p 

, Stream function (dimensionless) 

FIV(2,L,K) 

= 

u 

s 

, Streamwise velocity (dimensionless) 

FIV(3,L,K) 

ss 

^<p 

, Tangential velocity (dimensionless) 

FIV(4,L,K) 

= 

n 

, Static pressure (dimensionless) 

FIV(5,L,K) 

= 

I 

, Entropy (dimensionless) 

FIV(6,L,K) 

= 

® 

, Static temperature (dimensionless) 

FIV(7,L,K) 

= 

P 

, Density (dimensionless) 

FIV(8,L,K) 

= 

M 

, Mach number 


FIV(9,L,K) = 
FIV(10,L,K) = 


L=1 

@ 

J-1 

station 


L=2 

@ 

J 

station 


L=3 

@ 

J+1 

station 


K=1,KL 



streamlines 


FIPARM(l) 

= 

^r 

9 

Reference density (slugs/ff^) 

FIPARM(2) 

= 

T 

r 

9 

Reference temperature (deg R) 

FIPAEM(3) 

= 

P 

r 

9 

Reference pressure (psfa) 

FIPARM(4) 

= 

g 

9 

Gravitational constant (ft/sec^) 

FIPARM(5) 

■= 

P 

r 

9 

Reference viscosity (slugs/f t/sec) 

FIPARM(6) 

= 

C 

P 

9 

Specific heat constant pressure (f f^/sec^/deg) 

FIPARM(7) 

= 

C 

V 

9 

Specific heat constant volume (f t^/sec^/deg) 

FIPARM(8) 

= 

R 

9 

r\ O 

Gas constant (f t'^/sec^/deg) 
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C(6MM0N/FIVC/ (Cont'd) 


FIPARM(9) 

Pr,p 

, Prandtl number (turbulent) 

FIPASM(IO) 



FIPARM(ll) 

Ur 

, Reference velocity (ft/sec) 

FIPARM(12) 

N0PT7 

Number of stations stored in inviscid solver 

FIPARM(13) 



FIPARM(14) 

I0PT15 

First station 

FIPARM(15) 

I0PT16 

Last station 


11-49 



LIST OF VARIABLES IN C0MM0N/FLAGS/ 

Internal Flags 

N0PT0 where 0 = 1,17 , Flags to regulate calculation 

*See section 6.3 List of Flags for Options 
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LIST OF VARIABLES IN C0MM0N/FL0WI/ 


FG(1,K) 

FG(2,K) 

FG(3,K) 

FG(4,K) 

FG(5,K) 

FG(6,K) 

K 


0 . 


M 

Po 

U 

1, KL 


Internal Flow Variables 

, Inlet swirl angle (deg) 

, Inlet stagnation pressure (dimensionless) 

, Inlet stagnation temperature (dimensionless) 

, Inlet Mach number (dimensionless) 

, Inlet stagnation density (dimensionless) 

, Inlet magnitude of velocity (dimensionless) 

, Number of streamlines 
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LIST OF VARIABLES IN C0MM0N/F0RS/ 


Calculated Blade Force Variables 


F0RC(1,K) 

= 


, Streamwise force/area (dimensionless) 

F0RC(2,K) 

= 

-0- 

, Swirl force/area (dimensionless) 

F0RC(3,K) 


~s 

, Streamwise force/span (dimensionless) 

F0RC(4,K) 

= 

~<t> 

, Swirl force/span (dimensionless) 

F0RC(5,K) 

= 

^B 

, Blade dissipation/area (dimensionless) 

F0RC(6,K) 

= 

V 

, Blade dissipation/span (dimensionless) 

K 

= 

1,KL 

, Number of streamlines (dimensionless) 
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LIST OF VARIABLES IN C0MM0N/F0RS2/ 


Blade Force Variables 


FF(1,I,K) 

— 

M 

, Inviscid Mach number (dimensionless) 

FF(2,I,K) 

= 

A 

n 

, Inviscid static pressure (dimensionless) 

FF(3,I,K) 

= 

A 

0 

, Inviscid static pressure (dimensionless) 

FF(4,I,K) 

= 

A 

, Inviscid total temperature (dimensionless) 

FF(5,I,K) 

= 

A 

, Inviscid total pressure (dimensionless) 

FF(6,I,K) 

= 

A 

P 

, Inviscid density (dimensionless) 

FF(7,I,K) 

= 

A 

Us 

, Inviscid streamwise velocity (dimensionless) 

FF(8,I,K) 

= 

A 

, Absolute swirl velocity (dimensionless) 

FF(9,I,K) 

= 

A 

4> 

, Relative swirl velocity (dimensionless) 

FF(10,I,K) 


A 

R 

, Radius (dimensionless) 

FF(11,1,K) 

= 

A 

a 

, Absolute angle to axis (deg) 

FF(12,I,K) 

= 

A 

B 

, Relative angle to axis (deg) 

FF(13,I,K) 

= 

A 

I 

, Inviscid flow entropy (dimensionless) 

FF(14,I,K) 

= 

A 

u 

, Magnitude of absolute inviscid flow velocity 

(dimensionless) 

FF(15,1,K) 

= 

A 

, Loss coefficient (dimensionless) 

FF(15,2,K) 

= 

A 

, Blade entropy rise (dimensionless) 

FF(16,1.K) 

= 

A 

, Stream function (dimensionless) 

FF(17,1,K) 

= 

^L 

, Lift coefficient (dimensionless) 

FF(17,2,K) 

= 


, Drag coefficient (dimensionless) 

I 

= 

1 

, Upstream of blade row 

I 

= 

2 

, Downstream of blade row 

K 

= 

1, KL 

, Number of streamlines 
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LIST OF VARIABLES IN C0MM0N/FUNC/ 


Viscous Flow Variables 


F(1,1.K) 

= 

¥ 

, Stream function J=1 (dimensionless) 

F(2,1,K) 

= 

®o 

, Total temperature J=1 (dimensionless) 

F(3,1,K) 

= 

"o 

, Total pressure J=1 (dimensionless) 

F(4,1,K) 

= 

Ll 

, Entropy rise along (dimensionless) 

F(5,1,K) 

= 

A©o 

j 

, Total temperature rise along 

(dimensionless) 

F(6,1,K) 

= 


, Not used 

F(7,1,K) 

= 


, Not used 

F(8.1,K) 

= 


, Turbulent kinetic energy at J (dimensionless) 

F(9,1,K) 

= 

e" 

, Turbulent dissipation at J (dimensionless) 

F(10,1,K) 

= 

y^/p 
E r 

, Effective viscosity at J (dimensionless) 

F(1,N,K) 

= 

Y 

, Stream function (dimensionless) 

F(2,N,K) 

= 

Us 

, Streamwise velocity (dimensionless) 

F(3,N,K) 

= 


, Tangential velocity (dimensionless) 

F(4,N,K) 

= 

n 

, Static pressure (dimensionless) 

F(5,N,K) 

= 

I 

, Entropy (dimensionless) 

F(6,N,K) 

= 

® 

, Static temperature (dimensionless) 

F(7,N,K) 

= 

P 

, Density (dimensionless) 

F(8,N,K) 

= 

^ns 

, Streamwise stress (dimensionless) 

F(9,N,K) 

= 

M 

-e- 

, Tangential stress (dimensionless) 

F(10,N,K) 

= 

Q 

, Heat flux (dimensional) 

N = 2 

for station J 



= 3 for station J+1 


K = 1, KL 
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C^MMto/FUNC/ (Cont'd) 


FPAEM(L) 

= 

P 

r 

, Reference density (slugs/ft^) 

FPARM(2) 

= 

T 

r 

, Reference temperature (deg R) 

FPARM(3) 

= 

P 

r 

, Reference pressure (psfa) 

FPARM(4) 

= 

g 

, Gravitational constant (ft/sec^) 

FPARM(5) 

= 


, Reference viscosity (slug/f t/sec) 

FPARM(6) 

= 

S 

, Specific heat (f t^/sec^/deg R) 

FPARM(7) 

= 

Cv 

? 2 

, Specific heat (ft /sec /deg R) 

FPARM(8) 

= 

R 

2 2 

, Gas constant (ft /sec /deg R) 

FPARM(9) 

= 


, Turbulent Prandtl number 

FPAEM(IO) 

= 


, Not used 

FPARM(ll) 

= 

u 

r 

, Reference velocity (ft/sec) 

FPARM(12) 

= 

N0PT7 

, No. blocks stored in Inviscid solver 

FPARM(13) 

= 

N0PT8 

, No. blocks stored in viscous solver 

FPARM(IA) 

= 

I0PT15 

, First station 

FPARM(15) 

= 

10PT16 

, Last station 
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LIST OF VARIABLES IN C0MM0N/INTINP/ 


Mesh Parameters, Input Options, Debug Options 


IDBG0 

0 = 1,20 

9 

Debug options see section 4.3 

I0PT0 

0 = 1,20 

9 

Input options see section 4.2 

JFIRST 


9 

First streamwise station 

JL 


9 

No. of streamwise stations 

JLAST 


9 

Last streamwise station 

JLEDG 


9 

Blade leading edge 

JLPTS 


9 

No. wall input data points 

JSEP 


9 

No. stored blocks viscous solution 

JTEDG 


9 

Blade trailing edge station 

KDS 


9 

No. of steps/station 

KL 


9 

No. streamlines 

KLL 


9 

No. of input flow data points 
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LIST OF VARIABLES IN C0MM0N/MATRXA/ 

Variables for Viscous Flow Matrix Equations 
AE(I,J,K) = E , Element of E matrix in S0LVI 

AZ(I,K) = Zj , Element of Z vector in S0LVI 

BA(I,K) = B^ , See S0LVI 

FV(I,J,K) = F(I,J,K) , Solution from global iteration S0LVI 
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LIST OF VARIABLES IN C0MM0N/NGTERM/ 

Variables for Neglected Terms Calculation 
VARBL(I,J,K) , See subroutine VARFVN 

UN(K) = , Normal velocity (dimensionless) 

K = 1,KL 
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LIST OF VARIABLES IN C0MM0N/REALIN/ 


Input Flow Properties, Gas Properties 


ACHI 

X 


Clauser constant 

AKI 

K 

* 

Von Kasmen constant 

ALPl 

“l 

9 

Inlet swirl angle (deg) 

AMSl 

^Sl 

9 

Inlet mach number 

ANH 


9 

Power law ID wall 

ANT 

’^T 

9 

Power law OD wall 

API 

A+ 

9 

Van Driest constant 

CPRI 

C 

P 

9 

O 

Specific heat constant pressure (f t^/sec^/deg) 

CVRI 

Cv 

9 

Specific heat constant volume (f f^/sec^/deg) 

DDS 


9 

Distortion parameter 

DSHI 


9 

Displacement thickness ID wall (ft) 

DSTI 

S 

9 

Displacement thickness OD wall (ft) 

PRESO 

Po 

9 

Total pressure (psfa) 

PRLI 


9 

Prandtl number laminar 

PRTI 


9 

Prandtl number turbulent 

TEMPO 


9 

Total temperature (deg R) 

VISCRI 


9 

Reference viscosity (slug/ft/sec) 
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LIST OF VARIABLES IN C0MM0N/SPCFD/ 


F(K,J) 

P(K,J) 


Variables in Poisson Equation 


1/(PV)^ 


Coefficients of Poisson equation 
Stream function 
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LIST OF VARIABLES IN C0MM0N/SPCGD/ 


Variables for Streamline Curvature Calculation 


F(K) 

= 

i/(Pv)k 

, Coefficient of Poisson equation 

G(K) 

= 

v/(gP)j^ 

, Coefficient for velocity 

P(K) 

= 

Pk 

, Density ratio (p/Pj.) 

T(K) 


®K 

, Temperature ratio (T/Tj.) 

V(K) 

= 


, Metric coefficient 
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AINPUT(1,J,K) 

AINPUT(2,J,K) 

AINPUT(3,J,K) 

AINPUT(4,J,K) 

AINPUT(5,J,K) 

J 

J 

K 

AVE(1,J) 

AVE(2,J) 

AVE(3,J) 

AVE(4,J) 

AVE(5,J) 

AVE(6,J) 

AVE(7,J) 

AVE(8,J) 

AVE(9,J) 

AVE(10,J) 


LIST OF SYMBOLS IN C0MM0N/SPI0/ 

Initial Flow Variables, Mass Average Variables 


= Y/Y^ 

= n 

o 

= n 

= a 

= ®o 

= 1 

= 2 

= 1,KLL 

= 



P 

M 



Spanwise location (dimensionless) 

Total pressure (Ib/ft^ abs) 

O 

Static pressure (Ib/ft abs) 

Swirl angle (deg to axis) 

Total temperature (deg R) 

Upstream 

Downstream of blade row (dimensionless) 
Number of spanwise stations 
Mass flow (dimensionless) 

Average streamwise velocity (dimensionless) 
Average swirl velocity (dimensionless) 
Average entropy (dimensionless) 

Average entropy (dimensionless) 

Average static temperature (dimensionless) 
Average density (dimensionless) 

Average Mach number (dimensionless) 

Average total pressure (dimensionless) 
Average total temperature (dimensionless) 
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LIST OF VARIABLES IN C0MM0N/SPI^/ 


BINPUT(2,J,K) 

BINPUT(3,J,K) 

BINPUT(4,J,K) 

BINPUT(5,J,K) 

J 

J 

K 


n 

a 

0 

o 

1 

2 

l.KLL 


2 

Total pressure (Ib/ft abs) 

2 

Static pressure (Ib/ft abs) 
Swirl angle (deg to axis) 
Total temperature (deg R) 
Inlet flow 
Exit flow 

Number of spanwise stations 
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LIST OF VARIABLES IN C0MM0N/STRES/ 
Geometric Functions 


G/vl'^ 

JK-] 


J 

K-1/2 


G(XV) 


I bV 


[XV 

3n 

1 

3R 

XR 

3n 

1 

3R 

[XR 

3S 

V 

_3_ 


XG 9n 


J 

K-1/2 

J 

K-1/2 

J 

K-1/2 

J 

K-1/2 

(f)| 


G(9+I,K) 


1 — 

— (-) - 

1 3Vl 


1 XG 

3n W ' 

XY 3nl 

K-1/2 

[V., 

1_ (G'l - 

1 3_R^ 

J 

IXG 

dn \v > 

XR 3n 

K-1/2 

G(I 

,K) (3 J-1, 

1=1,9 



where 


X = dn/dn 
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LIST OF VARIABLES IN C0MM0N/STRMES/ 


Poisson Solver Stretching Parameter 
BP0IS = B , Stretching parameter 

BP0ISI = Bj , Input stretching parameter 
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LIST OF VARIABLES IN C0MM0N/SVARB/ 
Mesh Variables, Flow Variables 


AMACHR 


Mr 

, Reference mach number 

AMACHl 

= 

Ml 

, Average inlet mach number 

APRESl 

= 

^1 

, Average inlet static pressure (psfa) 

APROl 

= 

^01 

, Average inlet total pressure (psfa) 

AREAR 

= 

r2 

r 

2 

, Reference area (ft ) 

AREAl 

= 

^1 

2 

, Inlet area (ft ) 

ATEMPl 

= 

■^1 

, Average inlet static temperature (deg R) 

DELO 

= 

6'Vr 

o r 

, Average displacement thickness 

BETA 

= 

An 

, Transverse step size 

DR0UGI 

= 

d 

, Equivalent sand roughness (in. microns) 

DS 

= 

As 

, Streamwise step coarse grid 

DSS 

= 

ds 

, Streamwise step fine grid 

DSH 

= 

^H^^r 

, Inlet displacement thickness ID wall 

DST 

= 


, Inlet displacement thickness OD wall 

DYNPl 

= 

^1 

, Average inlet dynamic pressure (psfa) 

DZ 

= 

Az/r 

r 

, Step' size in axial distance (dimensionless) 

ERRBKN 

= 


, Error in back pressure iteration 

ERRPIN 

= 

e 

P 

, Error in normal static pressure 

GMRl 

= 

(Y-l)Mp 


GMR2 

= 

2 

yM 

r 
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LIST OF VARIABLES IN C0MM0N/SVARB/ 


RADR 

— 

r 

r 

REY 

= 

^-r 

RH0O 

= 


SL 

= 

"l 

USR 

= 

^r 

USTARH 

= 

u*/u 

USTART 

= 

u*/u, 
T ' 

WFLI 

= 

w 

I 

WFLO 

= 

w 

Z1 


^1 


Reference radius (ft) 

Reynolds number 
Not used 

Streamwise coordinate length 
Reference velocity (ft/sec) 
Friction velocity ID wall 
Friction velocity OD wall 
Input weight flow (Ib/sec) 
Calculated weight flow (Ib/sec) 
Axial length of duct (ft) 
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LIST OF VARIABLES IN C0MM0N/TETLIN/ 


TITLE(12) 


Input title card 
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LIST OF VARIABLES IN C0MM0N/TRBL/ 


AMUE(K) 

~ 


, Effective viscosity 

AMUL(K) 

= 

y/pr 

, Molecular viscosity 

AMUT(K) 

= 


, Turbulent viscosity 

TRB(1,1.K) 

= 

J-1 

\ 

2 

, Turbulence kinetic energy (k/u ) 

r 

TRB(1,2,K) 

= 

K 

2 

, Turbulence kinetic energy (k/u^) 

TRB(1,3,K) 

= 


2 

, Turbulence kinetic energy (k/u^ 

TRB(2,1,K) 

= 

K 

, Turbulence dissipation 

TRB(2,2,K) 

= 

E-' 

K 

2 

, Turbulence dissipation (d r/p^ (r^/u,^.) e) 

TRB(2,3,K) 

= 

(4)^ 

2 

, Turbulence dissipation (p^/y^/ (r^/u^) e) 

Cl 



, Constant 

C2 


C2 

, Constant 

SIGE 


°E 

, Constant 

SIGK 



, Constant 

NTVRB2 



, Flag 

CDAMP 


Cm 

, Constant 
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LIST OF VARIABLES IN C0MM0N/TVRBS/ 


Turbulent Viscosity, Conductivity 


DHF(1,K) 


i 

j-i 

Not used 


DHF(2,K) 

DPF(1,K) 


(Pp/u-) , 

K+1/2 

> 

Effective 
Not used 

viscosity 

DPF(2,K) 

= 

k- 

\ Prrp ’ 

Effective 

conductivity 
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LIST OF VARIABLES IN C0MM0N/TURB2/ 


DDSTRM 

= 

(A6*) 

u 

DEL0R 

= 

(AX/6*) 

DSIGM 

= 

m 

EVSVSM 

= 

® m 


m 


Minimum change in 6 
Maximum step size 
Maximum change in wall stress 
Maximum error in Us 
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LIST OF SYMBOLS IN C0MM0N/WLINE/ 


BLNE(1,L) 

BLNE(2,L) 

BLNE(3,L) 

BLNE(4,L) 

BLNE(5,L) 

JKLNE(1,L) 

JKLNE(2,L) 

LEAST 


Variables for 






Locating Arbitrary Line 
Radius (dimensionless) 

Axial location (dimensionless) 
Streamline distance (dimensionless) 
Streamwise coordinate (dimensionless) 
Normal coordinate (dimensionless) 

J index corner point 
K index corner point 
Number of input points 


1, LEAST 
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LIST OF SYMBOLS IN C0MM0N/ZF0U/ 


Variables for 


FI(1,L) 

~ 


> 

FI(2,L) 

= 

Us 

> 

FI(3,L) 

= 


9 

FI(4,L) 

= 

n 

9 

FI(5,L) 

= 

I 

9 

FI(6,L) 

= 

0 

9 

FI(7,L) 

= 

P 

9 

FI(8,L) 

= 

Fns 

9 

FI(9,L) 

= 

^n(|) 

9 

FI(10,L) 

= 

0 

9 

FI(11,L) 

where 

L = l.NLAST 

Un 

9 

Output 


Output Data Line Calculation 

Stream function (dimensionless) 
Streamwise velocity (dimensionless) 
Swirl velocity (dimensionless) 
Static pressure (dimensionless) 
Entropy (dimensionless) 

Static temperature (dimensionless) 
Density (dimensionless) 

Streamwise stress (dimensionless 
Tangential stress (dimensionless) 
Heat flux (dimensionless) 

Normal velocity (dimensionless) 

data points 


11-73 



LIST OF SYMBOLS IN C0MM0N/XF0U/ 


Parameters for Output Data Line Calculation 


FIl(l) 

— 


> 

Stream function (dimensionless) 

FI1(2) 

= 

Us(Sj,Nl) 

9 

Streamwise velocity (dimensionless) 

FI1(3) 

= 


9 

Swirl velocity (dimensionless) 

FI1(4) 

= 

n (Sj,Nl) 

9 

Static pressure (dimensionless) 

FI1(5) 

= 

I (Sj,Nl) 

9 

Entropy (dimensionless) 

FI1(6) 

= 

0(Sj,Nl) 

9 

Static temperature (dimensionless) 

FI1(7) 

= 

P 

9 

Density (dimensionless) 

FI1(8) 

= 

^ns 

9 

Streamwise stress (dimensionless) 

RI1(9) 

= 

^nif) 

9 

Tangential stress (dimensionless) 

FIl(lO) 

= 

Q (Sj,Nl) 

9 

Heat flux (dimensionless) 

FIl(ll) 

= 

Un (Sj,Nl) 

9 

Normal velocity (dimensionless) 

FI2(N) 

where 

= 

FIl(N) at Sj^j^ 

• \ 



Sj = streamwise coordinate at station J 
Sj_j_^ = streamwise coordinate at station J+1 

= normal coordinate for output data point L 
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7.0 DETAILED DESCRIPTION OF ADD CODE 


Name 

ALINE 

ALTMN 

AMF 

AMFIJIR 

AMINLT 

AMU 

BAMFQ(R 

BLDGQIM 

BLINE 

BLKDAT 

BLKRED 

BLPARM 

BPLUSR 

CALINV 

CASC 

CDS 

CFCGfLE 

CKINPT 

CLINE 

Ctf(/R 

COfQ^RST 

CgfQfRl 


7.1 List of Subroutines and Functions 

Object 

Find intersection of output data line with wall 
Main Program (see Section 5.1) 

Calculate Mach number from area ratio 
Find stream thrust average 
Compute Mach number from mass flow 
Compute molecular viscosity (function) 

Compute choked weight flow and stream thrust 
Locate blade centerline in (n,s) coordinates 
Calculate neighboring points on output line 
Set program constants in BLOCK DATA 
Read/write data blocks on data files 
Computer boundary layer parameters 
Rough wall integration constant 
Calculate inviscid flow solution 
Calculate cascade performance 

Calculate Roberts mesh distortion parameter (function) 
Calculate Coles' friction coefficient (function) 

Check input data for radial equilibrium 
Calculate (n,s) coordinates for point (R,Z) 

Interpolate coordinates in streamwise direction 
Controls flow of coordinate calculation 
Computer approximate coordinate functions 
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Name 

CClCiR2 

CCiCSRS 

C(iCiRA 

C0C1R5 

CPLXl 

CRCSSSI 

CR(^SS2 

DEFFIL 

DLINF 

DR0BRT 

DUMI0 

EC0INP 

ELINE 

ERPIN 

EWL0SS 

FAVER 

FC0LES 

FC0RCT 

FCPLX 

FETA 

FINTG 

FINVIS 


Obj ect 

Calculate wall curvature and arc length 

Compute coordinate functions 

Find Schwartz Christoffel parameters 

Interpolate wall curvature at station S 

Real and imaginary parts Schwartz-Christof f el transform 

Find intersection of two straight lines 

Search coordinate rectangle for crossing point 

Define random access data files (IBM version) 

Store coordinates of output data line 

Computation of derivative of Robert's transformation (function) 
Read/write random access data files (IBM version) 

Write input data 

Extrapolate for end points on output data line 
Compute error in normal pressure gradient 
Calculate end wall loss 

Solve for mass flow weighted average flow 
Compute Cole's velocity profile (function) 

Correct solution based on mass flow average 

Evaluate Schwartz Christoffel transformation 

Calculate distorted mesh to be used in POIS 

Integrate equations of Schwartz Christoffel transformation 

Calculate upstream/downstream force variables 
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Name 

Obj ect 



FLINE 

Find coordinates of output data line 


FLCSWIN 

Set initial flow conditions 



fncSrm 

Normalize input variables 



FCiRCE 

Compute blade force/span 



fcSrcl 

Computer blade force/volume 



F0UTP 

Interpolate solution on output data 

line 


FTHIK 

Calculate blade thickness 



GBLADE 

Calculate blade geometry 



GDUCT 

Prescribe duct shape 



INIT 

Set all C0MM0N block variables to zero 


INITF 

Initialize data file parameters for 

inviscid 

solution 

INITFV 

Initialize data file parameters for 

viscous 

solution 

INITQ 

Initialize data file parameters for 

Q array 


INTFRE 

Initialize freestream conditions 



IZER«S 

Set integer variables to zero 



MINVRT 

Invert MxN matrix 



0RTFUN 

Set up coordinate functions used in 

SfSLVI 


CSUTPUT 

Print title page 



P0IS 

Solve Poisson equation 



P«JISCF 

Setup coefficients for POIS 



P^IS0N 

Calculate axisymmetric streamline curvature 


PRT0UT 

Calculate and print terms neglected 

in ADD code equations 
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Name 

QINTER 

QPTMAX 

RDC0R 

READPF 

READPG 

REC0RD 

REDINP 

ReJBRTS 

SCURVA 

SLETE 

SL0TA 

SLTFL0 

SM00TH 

S0LVI 

STRESI 

SUBLAY 

TPRINT 

TRBLADE 

TRLETE 

TURB 

TURB2Q 

UBLAS 


Obj ect 

Interpolate curvature from P0IS to mesh to S0LVI mesh 

Calculate maximum total and dynamic pressure 

Read J and J+1 coordinate block 

Read P and F files in NIST word blocks 

Read variables for curvature calculation 

Calculate record numbers 

Read input parameters and variables 

Calculate distorted mesh using Robert's transformation(function) 

Calculate curvature from potential flow solution 

Find indices for strut control surface 

Construct duct with slots 

Compute slot inlet flow 

Smooth wall contour 

Integrate equations of motion 

Initial stress and heat flux 

Initialize K,E turbulence model in inner layer 
Print computer time 

Transform from stacking plane to duct plane 
Find leading and trailing edge in stacking plane 
Compute turbulent viscosity using algebraic turbulence model 
Calculate turbulent viscosity using K,E turbulence model 
Blasius' velocity profile (function) 
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Name 

UMACH 

VARFUN 

WBLEED 

WFITER 

VJRITPF 

WRTBLD 

WRTCAL 

WRTCKI 

WRTF0U 

WRTGDC 

WRTINP 

WRTS0V 

WRTSUM 

XH 

XT 

ZER0 


Object 

Calculate inviscid velocity 

Store variables to be used by PRT0UT 

Calculate perforated wall bleed 

Inlet weight flow iteration 

Store updated potential flow solution 

Write output for WBLEED 

Write output for CALINV 

Write output for CKINPT 

Print solution on output data line 

Write output summary - coordinate calculation 

Write input parameters and variables 

Print solution at selected stations 

Write output summary papers 

Calculate arc length hub wall 

Claculate arc length tip wall 

Set real variables to zero 
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7.2 Detail Description of Subroutines 


Subroutine ALINE (IWALL, JBLN, RX, ZX, XX) 


Obj ect 

Find intersection of output data line with the wall 
Options 

IWALL =1 OD WALL 

=2 ID WALL 

List of Symbols 

JBLN , 

(R1,Z1,R2,Z2) = (Rj^,Z^,R2,Z2) 

(RB1,ZB1,RB2,ZB2) = (R^^ , , ^2 , Z^) , 

RX,ZX = R,Z 

XX = X 

Theory 


Intersection occurs between JBLN, JBLN+1 

Neighboring points on wall (dimensionless) 

Neighboring points on input line (dimensionless) 

Coordinates of intersection (dimensionless) 

Distance along wall to intersection (dimensionless) 


The wall coordinates are searched for the intersection of the wall and output 
data line. This intersection must occur within the pair of points (RB1,ZB1) and 
(RB2,ZB2). If this intersection does not occur within the coordinate mesh, 
DIAGN0STIC No. 39 is printed and the program stops. If an intersection of the wall 
and the extrapolated output data line occurs, this intersection (RX,ZX) is printed. 
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Subroutine AMF (AA, AMI, AMG, Ail, ACPC, ACPI) 


Object 

Calculate Mach Number from Area Ratio 
Options 

AMG < 1 Subsonic Solution 
> 1 Supersonic Solution 


List 

of Symbols 



AA 

•= a/a^ 

Area ratio 


ACPC 

= c 

PC 

Pressure coefficient 

(compressible) 

ACPI 

II 

0 

M 

Pressure coefficient 

(incompressible) 

AM 

= M 

Mach number 


AMG 

= Mg 

Initial guess for Mach number 

AMI 

= M, 

Mach number at inlet 



Y , Ratio of specific heat 

Pq/P , Total to static pressure ratio 

Theory 

For one dimensional isentropic flow the area ratio is related to the Mach 
number by (see Ref. 1), 



M 


r 1 + Y-i M 


2 ^ Y+l 
2(y-1) 


L 1 + Yzi J 
2 ^ 


( 1 ) 


Eq. (1) has a minimum at M = 1. This subroutine first checks to see if the area 
ratio is less than this minimum. If the area ratio is less than the minimum it 
returns M = 1.0. If the area ratio is greater than the minimum, then it solves for 
M given A 2 /A and M^ using Newton's method. Either the subsonic or supersonic root 
is possible dependent on the value of the input AMG. Once the Mach number itera- 
tion converges, the pressure is given by 
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( 2 ) 




and the pressure coefficients are given by 



(3) 

(4) 


References 

1. Shapiro, A. H.: The Dynamics and Thermodynamics of Compressible Fluid Flow, 

Vol. 1, Ronald Press Co., New York, 1953. 
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Subroutine AMF0R(AW,AF,AA,AM,APR,APRO) 


Object 






Find stream 

thrust average properties 

Options 



None 


List 

of Symbols 


AA 

= A 

Area (a/(r^^)) 

AF 

II 

Stream thrust (f,p/(Pj. 

AM 

= ft 

Stream thrust average Mach number 

AN 

= N 

Mach number function 

AW 

% 

= w 

Reduced mass flow (M^W/(p^u^r^^) 

APR 

II 

Stream thrust average static pressure (p/p^) 

APRO 

II 

o 

Stream thrust average total pressure (p/p^) 


0o 

Mass flow average total temperature (T /T ) 

Theory 



The stream thrust and reduced mass flow (Ref. 1) can be related to one 
dimensional average flow parameters by the relations 


Then 
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Subroutine AMF0R (Cont'd) 


Eq. (3) has a minimum at M = 1 . If the input N is less then this minimum, 
the subroutine returns DIAGNOSTIC NO. 03. If N is greater than the minimum, 

Eq. (3) is solved for ft using Newton's method. The stream thrust average static 
and total pressure is then given by 


M^]/A 

(4) 

M 

(5) 


References 


1 . 


Shapiro, A. H. : The Dynamics and Thermodynamics of Compressible Fluid Flow 

Vol. 1, Ronald Press Co, New York 1953. 
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Subroutine AMINLT (WFL,P0,T0,ALI,AMI,P1) 


Obj ect 

Compute Mach number form mass flow 
Options 



None 


List 

of Symbols 


AINLET = A 

Area (ft^) 

ALI 

= “l ’ 

Inlet swirl angle (deg) 

AMI 

= M , 

Inlet Mach number 

AN 

= D , 

Mach number function 

PO 

= ^o ’ 

Total pressure (psfa) 

PI 

= P 

Static pressure (psfa) 

TO 

= To , 

Total temperature (deg. R) 

WFL 

= V , 

Weight flow (Ib/sec) 


"o ’ 

Total pressure ratio (Pq/P^) 


®o . 

Total temperature ratio (T^^/T 


Theory 

For average inlet conditions, the Mach number is related to the weight flow 
by the relation (see Ref. 1), 


- 1/2 

D = v/(gPj.C^Ilo®o 


y+1 


( 1 ) 
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Subroutine AMINLT (Cont’d) 


The Mach number function D(M) has a maximum at M = 1. If the specified weight flow 
is greater than this maximum, the subroutine returns DIAGN0STIC No. 2. If it is less, 
then Eq. (1) is solved for M using Newton's method. 
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Function AMU(T) 


Ob.i ect 

Compute Molecular Viscosity 
Options 

None 

List of Symbols 

AMU = v/Pj., Ratio of Molecular Viscosity (dimensionless) 

T =0 , Static temperature Ratio (dimensionless) 

Theory 

The molecular viscosity is computed according to Sutherland's formula (Ref. 1). 
The working fluid is assumed to be air. Accordingly, 

y_ 3/2 1 + 198.0/Tr 
yr 0+198.0/Tj. 

References 

1. Schlicting, H. : Boundary Layer Theory. 6th Ed., McGraw-Hill, New York, 1968. 
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Subroutine BAMF0R 


Obj ect 


Calculate choked weight flow and stream thrust 
Options 

None 

List of Symbols 


AM = ft 
APR = n 
APRO = 

0 

ARO = A 
AW - W 
DAPRO= An 

0 

DFO = AF 

T 

DTO = a'^ 

o 

'Xj 

FO = F^ 
STRT = f^ 
WC = w 

c 

WO = W 


0 0 


, Stream thrust average Mach number 
, Stream thrust average static pressure (P/P^.) 
, Stream thrust average total pressure (Pj/Pj.) 
, Area (a/Vj.^) 

, Reduced mass flow (M^w/|'p^u^r^^) 

, Mean deviation total pressure 

, Mean deviation stream thrust 

, Mean deviation total temperature 

, Stream thrust (f 'p/ (P^r^^) ) 

, Stream thrust (lb) 

, Choked weight flow (Ib/sec) 

2 

, Mass flow (w/(p^u^r^ )) 

, Mass flow average total temperature (T /T ) 

o r 


Theory 

v V ^ 

The stream thrust average quantities M, H, are calculated by subroutine 
AMF0R. The remaining quantities are given by 


W = /"pUdA 


( 1 ) 


Subroutine BAMF0R (Cont'd) 


'V 


% 

AF„ 


= / |ii + yMj.^ pu^j 
y | n + yMr 


dA 


+ yM.2 Pu2 - F^/A 


dA 


Ae 




■pu jeo ■ 0ol 


AH 


W. 


("o - \l) ' "ol 


Y+1 


•'•Vi 


where the "bar" denotes mass flow weighted average conditions 


( 2 ) 

(3) 

(A) 

(5) 

( 6 ) 
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Subroutine BLDG0M 


Object 

Locate blade centerline in (n,s) coordinates 
Options 


I0PT2 =0 No blades in duct 

= 1 Blades in duct 

List of Symbols 
C0MM0N BL0CK Variables 
Theory 

This subroutine uses subroutine FLINE to find the centerline in the (n,s) 
coordinates. Thus the input data from the blade stacking plane is transformed 
to the duct plane using subroutine TRBLD and stored in the input data line 
data block BLME (1,2) used by subroutine FLINE. The output (n,s) coordinates 
are then stored in the blade data array C0NST(I,L). At the completion of this 
calculation, the location of the upstream and downstream blade force calculation 
surfaces are determined by calling subroutine SLETE. 
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Subroutine BLINE (L,RBl,ZBl,RB2,ZB2) 


Object 

Calculate neighboring points on output line 
Options 
None 
Variables 

L , Point number 

RBl.ZBl = , Point at L-1 

RB2,ZB2 = R 2 ,Z 2 , Point at L 

Theory 

The points are read from an input table. 
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Subroutine BLKDAT 


Object 

Set program constants in BL0CK DATA 
Options 


None 



List of Symbols 



ACHI 

= X 

> 

0.016 (dimensionless) Clauser constant 

AKAPPA 

— K. 

» 

0.41 (dimensionless) Prandtl constant 

APLUS 

= A+ 

> 

26.0 (dimensionless) Van Driest constant 

CPR 

= Sr 

> 

5997. (ft^/sec^/deg R) Specific heat 

CVR 

^vr 


4283 (ft^/sec^/deg R) Specific heat 

EP 

= e 

> 

2.7182818 (dimensionless) 

GAMMA 

= Y 

» 

1.4 (dimensionless) Ratio of specific heats 

GASR 

= R 

9 

1714. (f t^/sec2/deg R) Gas constant 

GRAVR 

= g 

9 

32.2 (ft/sec^) Acceleration of gravity 

PI 

= IT 

9 

3.1415926 (dimensionless) 

PRESR 

- Pr 

9 

2117. (Ib^/ft^) Reference pressure 

PRL 

= Pr 

T 

9 

0.72 (dimensionless) Laminar Prandtl number 

PRT 

L 

= Prx 

9 

0.90 (dimensionless) Turbulent Prandtl number 

RH0R 

= Pr' 

9 

0.00238 (slugs/ft^) Reference density 

SNDR 

= C 

r 

9 

116.0 (ft/sec) Reference speed of sound 

TEMPR 

= T 

r 

9 

519.0 (deg R) Reference temperature 

TI 

= t 

9 

0.01745329 (dimensionless) 

VISCR 

= y 

9 

0.37 X 10"^ (slugs/ft/sec) Reference viscosity 
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Subroutine BLKDAT (Cont'd) 


Values of Parameter’s defined in C0MM0N/PARAM/ are also included for IBM 
and CDC computer programs. 
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Subroutine BLKRED (UNIT, RECSIZ, ADDR, BEGREC, NRECS) 


Object 

Reads NREC 'records' from file 'UNIT' beginning with record BEGREC. NRECS 
records are stored as a single block, beginning at ADDR. The Univac 1100 library 
I/O routine NTRAN is used in this subroutine. 


Variables 


UNIT 
RECSIZ = 
BEGREC = 
NRECS = 
ADDR 


logical unit// 

record size in words 

first record to read 

// of logical records to read 

beginning address to store the NRECS 

records read 


(Integer) 

(Integer) 

(Integer) 

(Integer) 

& RECSIZ 


Theory 

BLKRED (with entry BLKWRT) was developed to allow NTRAN compatibility with ANSI 
standard DEFINE FILE I/O operations.. In particular, a call to BLKRED with NRECS = 1 
is identical to a random access fortran read. 

In order to simulate DEFINE FILE I/O, it is necessary for BLKRED to maintain 
a list of pointers into the various disk files. The pointer list, DSKLOC, is in a 
common block | UNITS | which must be allocated in a static (root) segment. The loca- 
tion pointer and read size are used to position the disk for I/O access. After the 
I/O access, the pointer is positioned accordingly. 

BLKRED will issue a diagnostic message and cause program termination if either 
of two abnormal conditions are detected. 

1) the record // is negative 

2) NTRAN returns on error status less than zero (see UNIVAC FORTRAN V 
library routine NTRAN description on UNIVAC ASCII FORTRAN routine 
NTRAN^ description) 
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Subroutine BLPAEM 


Object 

Calculate boundary layer parameters 
Options 


N0PT4 


I0PT11 


1 Duct with centerbody 

2 Duct without centerbody 

0 Wall boundary OD wall 

1 Freestream boundary 


List of Symbols 


BLP(1,I) 

BLP(2,I) 

BLP(3,I) 

BLP(4,I) 

BLP(5,I) 

BLP(6,I) 

BLP(7,I) 

BLP(8,I) 

BLP(9,I) 

BLP(10,I) 

BLP(11,I) 

BLP(12,I) 

BLP(13,I) 

BLP(14,I) 


U_ , Inviscid edge :streamwise velocity (u_ /u^) 
°0 ^e r 

Ilg , Viscous edge static pressure (Pg/P^) 

, Viscous edge total pressure (Pqo/P ) 

’-’e ^ r 

0^ , Viscous edge static temperature (T^/T^) 

, Viscous edge Mach number 

U , Viscous edge streamwise velocity (u„/u_) 
s e r 

e 


P , Viscous edge density (p^/pj.) 
~e 

6/r^ , Edge of boundary layer 

* 

6/r , Displacement thickness 

e/r^ , Momentum thickness 
» Shape factor 

R , Reynolds number (p u„6/u„) 
e e ® 6 

e 

p , Inviscid edge density (p /p ) 

“e e r 

6,/r , Swirl displacement thickness 

4> r 
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Subroutine BLPARM (Cont'd) 


BLP(15,I) = 6,/r„ , Swirl momentum thickness 

<P ^ 

BLP(16,I) = , Viscous edge swirl velocity (u^ /u^) 

BLP(17,I) = , Inviscid edge swirl velocity (u,^^/Uy) 

WS “ streamwise vorticity (r^w^/u^) 

WP = , Viscous tangential vorticity (r j.o)|j,/u^) 

WV = fi , Viscous total vorticity (rj.u/u^) 

WSI = n , Inviscid streamwise vorticity (r*-u)c/u ) 

s r 

WPI = fijjj , Inviscid tangential vorticity (rj.o)^/u^) 

WI = , Inviscid total vorticity (r^,o)/u^) 

EWP = e , Vorticity test 


Theory 


g/g^ , Gap/Gap at wall 


For internal flows with thick boundary layers, normal pressure gradients, 
and swirl; the following boundary layer parameters may be defined. 



11-96 



Subroutine BLPARM (Cont'd) 



The edge of the boundary layer is determined by comparing the viscous and inviscid 
vorticity distributions where vorticity is given by 


. = 9 

^s 



V 

(5) 

; = vi_ 


(6) 

® r 8n 


'■ - (“/ 


(7) 


Then the edge of the boundary layer is defined by 


~ = 0.1 ( 8 ) 

u/ 5 
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FUNCTI0N BPLUSR (AKPULS) 


Object 


Calculate rough wall integration constant 
Options 

None 

List of Symbols 


AKPLUS = 
BPLUSR = B"''' 

* 


, Roughness Reynolds number (dimensionless) 
, Law of wall constant (dimensionless) 

, Universal velocity (dimensionless) 

, Universal distance (dimensionless) 


Theory 


The law of the wall as shown in Ref. 1 can be written 


U+ = i |ln y^ + B+ 
K 


( 1 ) 


where for a smooth wall 


and for a rough wall 


B+ = 2.05, K+ < 3.33 


B"*" = 3.4 - In K"*" > 3.33 


( 2 ) 

(3) 


References 

1. Schlicting, H.: Boundary Layer Theory 6th Ed., McGraw Hill Book Co., 1968. 
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Subroutine CALINV 


Object 

Calculate inviscid flow field solution. 
Options 

Calculate flow from J = JFIRS, JLAS 
IF (I0PT15.NE.O) JFIRS = I0PT15 
IF (I0PT16.NE.O) JLAS = I0PT16 
Calculate only for I0PT1 = 3 or 4 
N0PT5 4- 0 Error exit 


List of 

Symbols 




Same as 

CKINPT 




BINP(1, 

K) 

= P 

o 

t 

Total pressure (psf) 

BINP(2, 

K) 

= P 

f 

Static pressure (psf) 

BINP(3, 

K) 

= a 

S 

Swirl angle (deg) 

BINP(4, 

K) 

= T 

0 

9 

Total temperature (deg R) 

ITERAL 



9 

Swirl angle iteration number 

ERRA 


= E 

a 

9 

Local error in swirl angle 

ERRAM 


" ^nd 

9 

Maximum error in swirl angle 

Theory 






Given the swirl angle a, the analysis is identical to subroutine CKINPT. 
The calculation of the inviscid flow field requires also the solution of the 
angular momentum equation which is given by 


(RU ) = (R U ) 

(j) K I ((.I K 


( 1 ) 


where RU is the angular momentum at an arbitrary station and is the inlet 

angular momentum which is given. An outer iteration loop is then programmed to 
solve Eq. (1) to get the swirl angle a. With a known, the inner iteration loop 
is the same as subroutine CKINPT. This solution is obtained for each streamwise 
station J = JFIRS, JLAS. 
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Subroutine CALINV (Cont'd) 


References 


1. Anderson, 0. L. and D. E. Edwards: Extensions to an Analysis of Turbulent 

Swirling Compressible Flow in Axisymmetric Ducts UTRC Report R81-914720-18, 
February 1981. 
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Subroutine CASC (ARG. List) 


Object 

Calculate cascade performance 
Options 


N0PT5 = 0 
> 1 

I0PT2 = 2 
= 3 


ISHAPE 


continue calculation 

Error - stop calculation 

Not used - predict C , C 

i-< D 

Predict ol 2 , 

Blade shape index see input data 


= a. 


= a. 


List of Symbols 

Arg^ument—lis^ _(see_ 

ALIND = a . , 

1 

ALPS = Ug , 

ALPI 
ALP 2 

AMACH = M 

BCH0RD = c 

CL 
CD 

CX21 
GAM 
I0PT2 
ISHAPE 
N0PT5 


= Ct 


= C. 


D 

Usz/Usi, 


= Y 


Fi 8^ 1 ) 

Induced angle of attack (deg) 
Stagger angle (deg) 

Inlet angle (deg) 

Exit angle (deg) 

Inlet Mach number 
Chord (dimensionless) 

Life coefficient 
Drag coefficient 
Streamwise velocity ratio 
Ratio of specific heats 
Force option 
Shape option 
Error flag 
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Subroutine CASC (Cont'd) 


PHIC 

RH0CX 

S0LD 

TM 

WIND 

ZL0SS 

Cl Data 


, Circular arc camber angle (deg) 


= (PU„) o, /(PU ) , , Mass flow ratio 


= a 
= t 
= W. 

= z 


B 


, Solidity (dimensionless) 

, Thickness (dimensionless) 

, Induced velocity (dimensionless) 
, Loss coefficient 


Correlation (see Fig. 2) 


ALPMAX 

ALPO 

AMCR 

CB 

CBE 

CLE 

CLMAX 

CLP 

CLO 

DCLS 

CDLO 

EPC 

PHICE 


M 


= a 


= M 


cr 


= C 


Be 


= C 


Le 


= C 


LM 


= n 


= C 


LO 


= a 


= (dCL/da), 


= E 


= <t> 


ce 


. Cl(u^) = C^ 

> CL(oto) = 0 
, Critical Mach number 
, Design lift coefficient 
, Reference design lift coefficient 
, Reference lift coefficient 
, Maximum lift coefficient 
, Correlation parameter 

• >:lW> ■ 

, Correlation parameter 
, Lift curve slope 

, Deviation from ideal lift curve slope 
, Reference camber angle (deg) 
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Subroutine CASC (Cont'd) 


Cj) Data Correlation (see Fig. 3) 


ANG20 

“20 

9 

Angle to minimum C^q (deg) 

CDBP 

" "CBD 

y 

Power correlation parameter 

CDP 

’^CD 

y 

Power correlation parameter 

CDBO 

Sbo 

y 

Minimum drag (loss bucket) 

CDO 

“ ^DO 

y 

Minimum drag 

DMGl 

= Aa.2 

y 

Angle to 2 - C^^^ 

DANG2 

— Aa 

1 

y 

Angle to 2 - C^q 

DCBO 

^""DBO 

y 

Correlation parameter 

DCDO 

= "So 

y 

Correlation parameter 

f Ll 

to 02 , Zg (see Fig 

. 4) 

FS 

= S 

y 

Streamwise force coefficient 

FP 


y 

Tangential force coefficient 

P02P01 


y 

Total pressure ratio 

P2P1 

= P2/P1 

y 

Static pressure ratio 

T1 

- 

y 

Temperature ratio 

T2 

- Wh 

y 

Temperature ratio 

WS 

= s 

9 

Streamwise induced velocity (dimensionless) 

WP 

= 

9 

Tangential induced velociy (dimensionless) 

Cascade 

Correlation 

(see |^i£. 5) 

AKDELS 


9 

Shape parameter 

AMSIG 

= 

9 

Camber parameter 

AKDELT 

= St 

9 

Thickness parameter 

AlOO 

= So 

9 

Incidence angle (deg) 
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Subroutine CASC (Cont'd) 


AN 

= n 

> 

Power 


AKIT 

= K 

it 

> 

Thickness parameter 


AIMO 

= i 

mo 

) 

Minimum loss incidence 

angle (deg) 

AINCD 

= i 

9 

Incidence angle (deg) 


B 

= b 

9 

Power law correlation 


D 

= D 

9 

Diffusion parameter 


DELO 

= 5 

o 

9 

Deviation angle (deg) 


DELOO 

= 6 

oo 

9 

Deviation angle (J>^ = 0 

(deg) 

ZLOSM 

= Z 

BM 

9 

Minimum loss 



Theory 

The analysis used in this subroutine is described in detail in Ref. 1. The 
experimental data correlations for the cascade analysis was obtained from Ref. 2. 
Experimental data for the back pressure correction was obtained from Ref. 3. The 
lift and drag coefficient correlation parameters are based on thin airfoil theory 
described in Ref. 4. Lift and drag data were obtained from Ref. 4. The conversion 
of lift and drag to exit angle and loss is based on resolution of force vectors, 
(see Ref. 1). Since the angle of attack must account for induced velocity effects, 
the concept of momentum induced velocity was introduced (see Ref. 5). This in- 
duced velocity can be related to cascade back pressure effects (see Ref. 1). 

References 


1. Egolf, T. A., 0. L. Anderson, D. E. Edwards, and A. J. Landgrebe: An Analysis 

for High Speed Propellor-Nacelle Aerodynamic Performance Prediction. Vol. 1, 
Theory and Initial Application. UTRC Report R79-912949-19 , NASA Contract 
NAS3-20961, June 1979. 

2. Johnson, I. A. and R. 0. Bullock: Aerodynamic Design of Axial Flow Compressors. 

NASA SP-36, 1965. 

3. Mensing, A. E.: Investigation of the Effect of Back Pressure on Cascade 

Performance. UTRC Report R-0837-2, 1960. 

4. Abbot, I. H. and A. E. Von Doenhoff: Theory of Wing Sections, Dover 

Publications, 1959. 

5. McCormick, B. W.: Aerodynamics of VSTOL Flight, Academic Press, New York, 

1967. 
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Function CDS (DDS , BETA) 


Object 

Calculate Roberts' mesh distortion parameter 
Options 

None 

Input Variables 

DDS = Ratio of mesh distortion at wall An/An 

DETA = An , Mesh size at boundary - uniform mesh 

An , Distorted mesh size at wall 

Output Variables 


CDS C 


Roberts' mesh parameter 


Theory 

Let 

Then Roberts' 


C = 1/2 + € 


transformation can be written 


H-€ yZAi7-i) 


I + </> 
An = A 17 /DDS 


( 1 ) 

( 2 ) 

(3) 

(4) 


Eq. (2) through (4) can be solved iteratively for e as follows: 

e = 0 


j An + € 

^ ^ i+€-An 

€ = [<^ 2A17-1 - i]'' 


(5) 

( 6 ) 

(7) 
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Subroutine CDS (Cont'd) 


Convergence occurs when 


V +1 V 
€ - € 


< IE-04 


( 8 ) 
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Function CFC0LE (DEL, DSTR, H12, RET) 


Object 


Calculate Coles' wall friction coefficient 


Options 


List of Symbols 


, Von Karmen constant (.41) 


BPLUS 


, Constant of integration (2.05) 


CFC0LE 


, Friction coefficient 


, Boundary layer thickness (dimensionless) 
, Displacement thickness (dimensionless) 

, Shape factor 
, Reynolds number 


= U /U' 


, Velocity ratio 


Coles' skin friction law, Ref. 1, can be written 




= 


Eq. (1) can be solved for U /U by successive substitution. 


References 


1. Coles, D. E. : The Turbulent Boundary Layer in a Compfessible Fluid, Rand 

Report R-403-PR, 1962. 
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Subroutine CKINPT 



Options 

I0PT1 

i 4 

N0PT5 

^ 0 

I0PT5 

= 2 

IJ 

= 1 

IB 

= 2 

INEX 

= 1 

INEX 

= 2 

IDBG13 

= 1 

WFL0W 

= 0 

WFL0W 

> 0 

List of Syi 

nbols 

ERR 

= e 

EPS 

“ ^o 

FG(2, K) 

= M2) 

FG(4, K) 

= M2) 

ITER 

= V 

PHMAX 


PSll, PS12 

V 

= 

WF 

= W^ 


Do not calculate 
Error exit 

F0RCE data equals FL0WIN data 

Inlet flow data 

Force data 

Upstream data 

Downstream data 

Debug printout 

Static pressure check only 

Pressure check and weight flow iteration 

, Error in interaction 
, Minimum error 
, Equation 4 
, Equation 5 
, Iteration number 
, Maximum pressure possible (psfa) 

, Upper and lower bound air stream function (dimensionless) 
, Weight flow iteration (Ib/sec) 


II-lll 



Subroutine CKINPT (Cont'd) 


WFL0W 

= W 

WMAX 

= W 


max 

WMIN 

= w 


mxn 

XL 


XM 


XU 


XI 

= x,v,: 

PSIHT,PSIT 

II 

Theory 



Input weight flow (Ib/sec) 

Maximum weight flow possible (Ib/sec) 
Minimum weight flow possible (Ib/sec) 
Lower bound on X (dimensionless) 

X for choked flow (dimensionless) 

Upper bound on X (dimensionless) 
Iterative values for X (dimensionless) 
Value of stream function (dimensionless) 


Input data for the total pressure, static pressure, swirl angle, and total 
temperature must satisfy the continuity equation, and the radial momentum equa- 
tion. If these equations are not satisfied, the static pressure is adjusted. 

The solution of these equations can be obtained by a transformation of variables, 
Let 


y-i 


( 1 ) 


and 


0(17) = 2 ^ ^ cos^ “ ITr ^ 

' Xv on XR on ■> 

and then the radial momentum equation becomes 


( 2 ) 


dx 
d t 


[a (77) 


y-i 

r dv 


(in Rq)] X = a ( 77 ) 


( 3 ) 


which is an ordinary first order "linear" equation. The solution is given by 


( 77 ) = exp 


a (y) d Y + 


y-i 



( 4 ) 


V 

<t (77) = Jq(y) ^(y) dY ( 5 ) 

0 
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Subroutine CKINPT (Cont'd) 


X = (Xq+ <I>(tj)) /(^ {77) 


( 6 ) 


where is the hub static pressure ratio. The continuity equation becomes 


_ n 0 G COS g ^J/2 

dr) ’ Mf Vy^ 


(7) 


and 


w (77) = 2 77- Cr g 'J' {77) 

The constant X is determined by the boundary condition 


( 8 ) 


w ( I ) = w 


(9) 


using an iteration scheme described below. 
First let us examine the function 


f (X) = X (l-X)'^^ 


( 10 ) 


f(X) has a maximum at 


X 


_ 7^ 
- x+i 


( 11 ) 


Hence from Eq. (1) 



2 

r + 1 


Y 

r- 1 


( 12 ) 


Equation (12) is precisely the condition for choked flow when M = 1. Substitu- 
tion of Eq. (12) into Eq. (10) yields an a-priori condition for the maximum 
weight flow possible (i.e., the choked flow condition). For a subsonic solution 
we have the condition 
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Subroutine CKINPT (Cont'd) 


Furthermore, it is noted that we can find a-priori such that 

^ <Xl <Xo <Xu <1 

f (Xj <f (Xm) 
f (Xq) > 0 


( 13 ) 


(14) 


by substituting Eq. (6) into Eq. (10). Thus and X^ are the bounds for 

choosing subsonic solutions with no reverse flow. The iteration scheme then 

consists of narrowing the bounds of X, and X until convergence occurs. This 

Ij u 

procedure is illustrated in Fig. 17. 


I/ + 1 1/ * • 1 IV V \ 

X =X, + -p iT Xp -X, 


I ' 2 1 / 

(15) 

If X,*' =X" j '^,''='1' 

(16) 

If {'^*'■^'>'1') X2"=X*' j'F2*' = 'I' 



and is obtained by integrating Eq. (7) with X^ = X ' substituted into 
Eq. (6). Convergence occurs when 


ix"^'-x" 


< €, 


(17) 


Once X is known, the static pressure is obtained from (1) and substituted 
for the input static pressure. 

References 


1. Anderson, 0. L. and D. E. Edwards: Extensions to an Analysis of Turbulent 

Swirling Compressible Flow in Axisymmetric Ducts. UTRC Report R81-914720-18, 
1981. 
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Subroutine CLINE (JJ, KK, RB, ZB, SB, ANB, XB) 


Object 

Calculate (n,s) coordinates for point (R,Z) 
Options 
None 

List of Symbols 


JlK 


, Corner point on grid 

RB, ZB 

= R,Z 

, Dependent coordinates (dimensionless) 

SB, ANB 

= s,n 

, Independent coordinates (dimensionless) 

XB 

= X 

, Distance from inlet (dimensionless) 

Theory 



For 

a first 

order expansion we have 


R R 1 + 2 R 


'J,K 




K ''J + l, K + l ”J,K, ^J,K+|J AS 


Ifr, r, TO-n 

T L^J.K + I ^J + l,K + l“^J,K ~^J + I,KJ An 




^ ~ Zj.K 2 [^J + I.K + 2j -Zj ,< + |] 


AS 


+ ![? 4--7 _-7 ■» 

2 L J.K + I ^J+I,K+I ^J.K ^J+I,kJ 


( 1 ) 


( 2 ) 


These two equations can be solved for s,n and X may be obtained by interpolation. 
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Subroutine C00R(JS,KS) 


Object 

Interpolates coordinates in streamwise direction. 
Options 

I0PT9=O Computes approximate coordinates 
I0PT9^ Coordinates stored on drum 


List of Symbols 



DSTEP 

= AS 

, Streamwise step size (dimensionles 

,s) 

DX 

= AX 

, Interpolation between JS and JS-1 

stations 

DXIffiXT 


, Interpolation for next step 


JDRUM 

= 

, Drum unit number 


JS 

= 

, station stored on drum 


KS 

= 

, station interpolated between 

JS 

ZHUB 

= Zjj 

, Axial station hub (dimensionless) 


ZNEXT 


, Next axial station (dimensionless) 


ZTIP 

= Zip 

, Axial station tip (dimensionless) 


Theory 




Let 

the streamwise 

coordinate S be given by 



S = AS(JS-1) + dS-(KS-l) 


( 1 ) 


where 


AS = SlAJLAST-i) 
dS = AS/KDS 


( 2 ) 


(3) 
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Then if station 1 is at JS-1 and station 2 at JS, a simple linear interpolation of 
the coordinates may be made from those stored on the dram. ZNEXT is the axial 
location of the KS+i station. 
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Subroutine C00RST 


Ob.j ect 

Controls flow for coordinate calculation 
Options 

I0PT13=O No slots in duct 
=1 Slots in duct 

I0PT14=O ten streamlines are calculated 

>0 I0PT14 streamlines are calculated 


List of Symbols 



DSTEP 

= AS = DS 

5 

Strearawise step size (dimensionless) 

DX 

= ATI 

9 

Normal coordinate step size (dimensionless) 

DYl 

= 6Y 

9 

First derivative (dimensionless) 

DY2 

= 

9 

Second derivative (dimensionless) 

ICK 

= 0,1 

9 

Flag; no overlap, overlap 

II 

= 1,2 

9 

Flag; start integration, continue integration 

ISLTl, 

ISLT2 

= 

9 

First slot number, second slot number 

IW1,IW2 

= 

9 

First slot wall, second slot wall 

JC0UNT 

= 

9 

Strearawise station counter 

KLH0LD 

- 

9 

Number of streamlines to interpolate 

KN 

= 

9 

Number of streamlines to integrate 

MSL0T 

= 

9 

Slot counter 

NSL0T 

= 

9 

Total number of slots 

RA(I,J) 

= R(I,1,J) 

9 

Temporary storage for wall coorindate 

X 

= X 

9 

Normal coordinate (dimensionless) 


II-U3 



XM 

= X 

, Interpolation distance (dimensionless) 

X2 

= X2 

, Midpoint of three point difference (dimensionless) 

Y 

= T 

, Function to be interpolated (dimensionless) 

Y1 

= y(Xr)=Yi 

, Known values of Y (dimensionless) 

Y2 

= Y(X2)=Y2 

, Known values of Y (dimensionless) 

Y3 

= Y(X3)=Y3 

, Known values of Y (dimensionless) 

Z 

= Z 

, Axial distance (dimensionless) 

ZZ1,ZZ2 

= Z]_, Z 2 

, Location of adjacent slots (dimensionless) 

Theory 




This subroutine controls the calculation flow for the coordinates according to 
flow chart Fig. 6. The basic calculation scheme with slots is to calculate the 
streamlines through successively larger ducts and' storing only those coordinates 
satisfying the condition (Z]_^Zs:Z 2 ) as shown in Fig. 7. 

In addition, it was determined that only KN streamlines need be calculated by 
integration, the remainder up to KL streamlines may be calculated using a linear 
interpolation. 
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FIG. 6. FLOW CHART FOR SUBROUTINE COORST 
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RADIAL COORDINATE R 


AXIAL COORDINATE Z 



RG. 7. CALCULATING DUCTS WITH SLOTS 
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Subroutine C00R1 (KSS , JSS) 


Object 

Complete approximate coordinate functions. 

Option 

N0PT3 = 1 Construct coordinate 

= 2 Calculate coordinates at KSS, JSS 


List of 

Symbol! 

3 



R0(1.K) 

= 

R 

9 

Radius (dimensionless) 

R0(2,K) 

= 

dR/dZ 

9 

Derivative of radius (dimensionless) 

R0(3.K) 

= 

d^R/dZ^ 

> 

Second derivative of radius (dimensionless) 

R0(4,K) 

= 

Z 

9 

Axial distance (dimensionless) 

R0(5,K) 

= 

V 

9 

Metric scale coefficient (dimensionless) 

R0(6,K) 

= 

dV/dS 

9 

Derivative of metric scale coefficient (dimensionles 

R0(7,K) 

= 

Yt 

9 

Distance across duct (dimensionless) 

R0(8,K) 

= 

s 

9 

Streamwise coordinate (dimensionless) 

R0(9,K) 

= 

0 

m 

9 

Mass flow bleed (dimensionless) 

R0(1O,K) 

= 

w 

9 

Wall temperature (dimensionless) 

R0S(I) RO(I,K) 





Note; Three arrays defined where (j> takes on the value H, M, T 

4) = H, Hub wall 
= M, Mean line 
= T, Tip wall 

KSS , Coarse grid station 

JSS , Fine grid station 


s) 
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Subroutine C00R1 (Cont'd) 



The approximate coordinates are calculated by geometric construction using 
the relations derived in Ref. 1. These approximate coordinates are used as a first 
guess by subroutine C00R4. 


References 


1. Anderson, 0. L.: User's Manual for a Finite Difference Calculation of Turbulent 

Swirling Flow in Axisymmetric Ducts with Slot Cooled Walls. Vols. 1 & 2. 
USAAMRDL-TR-74-50, 1974. 
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Subroutine C00R2 


Object 

Calculate wall curvature and arc length 


Options 


None 



Variables 



RM(1,J) 


.Radius ID wall (dimensionless) 

RM(2,J) 


.Axial length ID wall (dimensionless) 

RM(3,J) 


.Arc length ID wall (dimensionless) 

RM(4,J) 


.Curvature ID wall (dimensionless) 

RM(5,J) 

T 

.Radius OD wall (dimensionless) 

RM(6,J) 

1 

Zti 

.Axial length OD wall (dimensionless) 

RM(7,J) 

Xti 

.Arc length OD wall (dimensionless) 

RM(8,J) 

Kti 

.Curvature OD wall (dimensionless) 

RM(9,J) 



RM(10, J) 



Theory 



The arc 

length is 

given by 


dXj=yrdR);+(dz)/ 


Xj =/dXj 


(dZ)j = - Zj 


( 1 ) 

( 2 ) 

(3) 
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Subroutine C^0R2 (Cont'd) 


and the curvature is given by 


dz d^R dR d^z 


dx "d^ dx "dx? 



3/2 


where the derivatives are obtained using a central difference formula. 


( 4 ) 
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Ob.iect 


Subroutine CgS0R3(KSS ,JSS) 


Compute Coordinate Functions 
Options 

L0P=1 Compute initial constants 
=2 Integrate one step 
=3 Do not integrate 


List of Symbols 



AO 

-.Aq 

5 

Inlet height W plane (dimensionless) 

ALO 

= “o 

9 

Inlet angle W plane (dimensionless) 

ANO 

= "o 

9 

Inlet height Z plane (dimensionless) 

SLO 

— q O 

- Sp 

9 

Initial coordinate length (dimensionless) 

RO 

= ^o 

9 

Inlet radius Z plane (dimensionless) 

VO 

= Vo 

9 

Inlet metric scale coefficient (dimensionless) 

SO 

= si 

9 

Inlet strearawise coordinate (dimensionless) 

XBO 

= Xo 

9 

Real part of C^^ (dimensionless) 

YBO 

= Xo 

9 

Imaginary part of (dimensionless) 

ZETAO 

= §o 

9 

Constant of integration real (dimensionless) 

ETAO 

= ^o 

9 

Constant of integration imaginary (dimensionless) 

Theory 




The 
to get R 

Schwartz 

,Z,V,S,N 

-Christoffel transformation (Ref. 1) is integrated along streamlines 
once the location of the poles are known. 


References 

1. Anderson, 0. L.: User’s Manual for a Finite Difference Calculation of Turbulent 

Swirling Flow in Axisymmetric Ducts with Slot Cooled Walls. Vols. 1 & 2, 
USJ!JVMRDL-TR-TU-50, 19Ti+. 
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Subroutine C00R4 


Ob.lect 

Find Schvartz-Christof fel Parameters 
Options 
None 

List of Symbols 


AX(1,J) 

AX(2,J) 

AX(3,J) 

AX(4,J) 

AY(1,I,J) 

AY(2,I,J) 

ay(3,i,J) 

AY(4,I,J) 

AY(5,I,J) 

AY(6,I,J) 

AY(7,I,J) 

AY(8,I,J) 

AERR(1,J) 

AERR(2,J) 

AO 

ALO 

ANO 

SLO 


= X- 


H 


= Xr 


= a. 


H 


= a, 


T 

,v 
H 

V 
H 

= Vh 


= S 
= X- 


= R. 
= S 
= X, 


= Vrr 


= Rj 


v +1 

= Sh - Sh^ 


= s 


v+1 




= h. 


= Or- 


= n. 


= Sl 


, Distance along hub wall (dimensionless) 

, Distance along tip wall (dimensionless) 

, Wall angle hub (deg) 

, Wall angle tip (deg) 

, Strearawise coordinate hub (dimensionless) 

, Distance along hub wall (dimensionless) 

, Metric scale coefficient hub (dimensionless) 

, Radius hub (dimensionless) 

, Strearawise coordinate tip (dimensionless) 

, Distance along tip wall (dimensionless) 

, Metric scale coefficient tip (dimensionless) 

, Radius tip (dimensionless) 

, Error in (dimensionless) 

, Error in (dimensionless) 

, Inlet height W plane (dimensionless) 

, Inlet angle W plane (dimensionless 
, Inlet height Z plane (dimensionless) 

, Initial guess of coordinate length (dimensionless) 
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RO 

= ^o 

, Inlet radius in Z plane (dimensionless) 

VO 

0 

> 

II 

, Inlet metric scale coefficient (dimensionless) 

SO 

= So 

, Inlet strearawise coordinate (dimensionless) 

XBO 

= ^o 

, Real part of constant Op (dimensionless) 

YBO 

= Yo- 

, Imaginary part of constant Cp (dimensionless) 

zetao 

0 

II 

, Constant of integration (real) (dimensionless) 

ETAO 

= Ho 

, Constant of integration (imaginary) (dimensionless) 

Theory 




This subroutine solves for the location of the poles in the Schwartz-Christoffel 
transformation using the method described in Ref. (l) and (2). The following pro- 
cedure is used: 


1. Construct a table of arc length along tne walls and wall to each of the 
JL input points. 

2. An initial guess is made for the location of the poles using the approxi- 
mate solution from subroutine C00R1. 

3- The Schwartz-Chrisotf fel transformation is integrated along the walls to 
obtain arc length along the wall to each pole. 

U. A new guess for the location of the poles is made by comparing arc length. 


References 


1. Anderson, 0. L. : Calculation of Internal Viscous Flows in Axi symmetric- Ducts 

at Moderate to High Reynolds Numbers, Int. J. of Computer and Fluids, Vol. 8, 

No. 4, December I 98 O. P. 391-^11 • 

2. Anderson, 0. L.: User's Manual for a Finite Difference Calculation of Turbulent 

Swirling Flows in Axisymmetric Ducts with Slot Cooled Walls. Vols. I & II, 

usaamrdl-tr-74-50, 1974 . 


II-12G 



Subroutine C00P5 (XH, XT, CURVH, CURVT) 

Object 

Interpolate wall curvature at Station S 
Options 


None 




Var iables 




XH 


X^(S) 

, Arc length ID wall (dimensionless) 

XT 


X^(S) 

, Arc length OD wall (dimensionless) 

CURVH 


K^(S) 

, Curvature ID wall (dimensionless) 

CURVT 


K^(S) 

, Curvature OD wall (dimensionless) 

RM(3, J) 


x^iCJ) 

, Table of X^j (dimensionless) 

RM(4, J) 


%l(J) 

, Table of Kjj (dimensionless) 

RM(7, J) 



, Table of Xj. (dimensionless) 

RM(8, J) 


Kti(j) 

, Table of Xp (dimensionless) 

Theory 




The values 
values of K^CS), 

of ^(S) and X,j,(S) are input. Then the tables are searched and 
K.p(S) are calculated by linear interpolation. 
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Subroutine CPLXl (Xl ,Y1 ,XB1 ,YB1 ,N1 ,N2 ,L0P ) 


Object 

Evaluates real and imaginary parts, Schwartz-Christoffel transformation 
Options 


L0P=1 Wz = X + iy = dW/dZ 
=2 W 2 Z = X + iy = d2w/dz2 


List of Svmbols 





N1 

= Ni 

9 

Product index 


N2 

= N2 

9 

Product index 


N3 

= Nj_+1 

9 

Product index 


XI 

= X 

9 

X coordinate in Z plane 

(dimensionless ) 

Y1 

= Y 

9 

Y coordinate in Z plane 

(dimensionless ) 

XBl 

= X 

9 

X coordinate in Wz plane 

or Wzz plane (dimensionles 

YBl 

= Y 

9 

Y coordinate in Wz plane 

or Wzz plane (dimensionles 

Theory 






This subroutine evaluates the real and imaginary parts of the Schwartz-Christoffel 

transformation. 
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Subroutine CR(^SS1 (XI , Y1 ,X2 ,X2 ,XB1 , YBl ,XB2 , YB2 , S , T) 


Obj ect 


Find intersection of two straight lines 
Options 

If 0 ^ s - 1.0 lines cross within points on line A 
If 0 - t 1 1.0 lines cross within points on line B 

List of Symbols 


(XI, Yl), (X2.Y2) 
(XB1,YB1),(XB2,YB2) 
S 
T 


= (Xi.Yi),(X2,Y2) 
= (Xi,Yi),(X2,Y2) 
= s 
= t 


, Pair of points line A (dimensionless) 
, Pair of points line B (dimensionless) 
, Parameter line A (dimensionless) 

, Parameter line B (dimensionless) 


Theory 

The equations for the two straight lines can be written in parametric form 

X = (Xg-X, )■ S+ X| 


Y= (Y 2 -Y,) s-*-Y, 
X = (X2 - X,) - t + X, 
Y = (Yg - Y, ) • t + Y, 


These equations can be solved at the point of intersection for s and t. Then if 

0< S < I (2) 

line B crosses line A between the pair of points on line A. If 

0< t < I (3) 

then line A crosses line B within the pair of points on line B, 
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Subroutine CR0SS2(K,M0LD,RB1,ZR1,RB2,ZB2,SX,TX,MNEW) 


Object 

Search around coordinate rectangle for crossing point. 
Options 

If M = M0LD Skip check for crossing point 


0 i ss < 1 lines cross within 

0- < XT < 1 lines cross within 

List of Symbols 


J,K 




MNEW, 

M0LD 



Rl.Zl 

,R2,Z2 

= kl»Zx,R 

2*^2 

RB1,ZB1,RB2, 

,ZB2 = R]_, 

* ^2 * ^2 

RX,ZX 


= Rx, 


SS.TT 


H 

II 


SX,TX 


Sx> 

Tx 


Theory 

This subroutine searches around the 
of line B . 


points on line A 
points on line B 

, Corner point on grid rectangle 

, New, old crossing side 

, Points on line A - coordinate grid 
(dimensionless) 

, Points on line B (dimensionless) 

, Coordinates for intersection (dimensionless) 
, Parametric coordinates for lines A and B 
, Parametric coordinates for intersection 

coordinate rectangle for the intersection 
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Subroutine DEFFIL 


Obj ect 


Define random access data files 
Options 


None 

Theory 

This subroutine is used only by the IBM version of the ADD code to define 
random access data files using ANSI standard DEFINE FILE I/O. Details of the data 
files are described in Section 5.5 of the Manual. 
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Subroutine DLINE (L , SBl , ANBl ,XB1 , JBl , KBl) 


Object 

Store coordinates of output data line 
Options 
None 
Variables 


JBl, KBl 
L 

SBl, ANBl = s^, nj^ 
XBl = 


Index of corner point 
Index of input data point 
Coordinates of input data point 
Distance from inlet of input data point 


Theory 

This subroutine stores coordinates of input data line at the point L. 
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Function DROBRT (C, ETA, L0P) 


Object 

(Compute derivative of Roberts' transformation 


Options 

L0P = 0 wall - wall boundary 

L0P = 1 wall-f reestream boundary 

10P = -1 f reestream-wall boundary 


Input Variables 
C = C , 

ETA = n , 

L0P 

Output Variable 
DROBRT = 9n/3n 


Distortion parameter 
Input variable 
Option 

Output variable 


Theory 


The transform of the Roberts' Stretching for a distorted mesh is given by 


d 7 ) 
d n 



c in 


/ C-H/2 \ 
\ C - 1/2 / 


1-1 


( 1 ) 


(V-'^2)] 

where the options are 


^ }l<^p = o ^ = 

n=nJ” n=2n 


7} = 7j/2 


( 2 ) 


7j' = (\+7i)/2 . 

n = 2n'-i = 


We note that 


0 < 17 < 1.0 0 < n' < 1.0 


(3) 
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Subroutine DUMI0 (Arg. List) 


Obj ect 


Read/write to random access data files 


Options 


I«llTYPE = 1 
I0TYPE = 2 


Argument List 


I0TYPE 

UNIT 

RECSIZ 

ADDR 

JREC 

STATUS 


Theory 


Write to file 
Read from File 


Read/write option 

Logical unit number 

Record size in words 

Beginning address 

Number of logical records 

Switch indicating completion of l/O 


This subroutine reads or writes data blocks to random access files and is called 
only by subroutine BLKRED. It is used only in the IBM version of the ADD code as a 
transition between UNIVAC NTRAN I/O and ANSI standard DEFINE FILE I/O. 
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Subroutine EC0INP 


Object 

Write input data 
Options 

None 

Variables 

See subroutine REDINP 
Theory 

This subroutine prints and labels the input data according to the output 
format described in Section 4.3. 
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Subroutine ELINE (TTT,SXT,XXT,TTH,SXH,XXH) 


Object 

Extrapolate for end points on output data line 
Options 

None 

List of Symbols 


SXH, 

SXT 

” » Sfji 

, Streamwise coordinate ID and OD walls 

TTH, 

TTT 

= Tjj.Tt 

, Parametric variable ID and 

OD walls 

XXH, 

XXT 

= Xh,x^ 

, Arc length ID and OD walls 

(dimensionless) 


Theory 

The parametric variable T interpolates along the input data line such that 

T - !n~f' _ >-^1 

■ Sg-S, ■ Xg-X, (1) 


J = _ ^LL_AST-I 

^LLAST “ ^LLAST-I 


X -y 

= _ T ^LLAST-I 
^LLAST ~^LLAST-I 


( 2 ) 


Since Tjj, Xjj, S^, ^2 are known, Sj^ and Xj^ may be determined from Eq. (1) . 

Likewise, Sj^lAST’ ^LLAST determined from Eq. (2). 
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Subroutine ERPIN (II) 


Object 


Check 

error in normal pressure 

Options 


II = 1 

Initial plane 

II = 2 

J-1 plane 

II = 3 

J plane 

Theory 



This subroutine integrates the normal momentum equation at the Ilth plane to 
obtain the static pressure difference OD wall pressure minus ID wall pressure. 

This static pressure difference is compared to the input static pressure difference. 
If the difference between calculated and input static pressure difference is 
greater than one percent returns DIAGNOSTIC NO. 09 and the calculated static pres- 
sure difference required to correct the input data. 
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Subroutine EWL0SS (Arg. List) 


Obj ect 

Calculate end wall loss 
Options 

None 


Argument List 
AVRH, AVRT 

a, 

^^’S2 


9 

Axial velocity ratio 

C^RH, C0RT 



9 

Chord (dimensionless) 

DELSIH, DELSIT 

* 

^ih’ 

* 

^2H 

9 

Displacement thickness (dimensionless) 

RH0IH, RH0IT 

a, 

^IH’ 

a, 

^IT 

9 

Density (dimensionless) 

TEMH, TEMT 

ih’ 

IT 

9 

Temperature (dimensionless) 

UIH, UIT 

"iH’ 

^'iT 

9 

Total upstream velocity (dimensionless) 

U2H, U2T 

'\j 

^2H’ 

^2T 

9 

Total downstream velocity (dimensionless) 

Y 

Y 


9 

Distance from wall (dimensionless) 

YT 



9 

Duct height (dimensionless) 

ZEWL0S 

^EW 


9 

End wall loss coefficient 

Internal Variables 

ANl, AN2 

Ni, 

^2 

9 

Power law 

CF 



9 

Friction coefficient 

DELIH, DELIT 

^IH’ 

'^IT 

9 

Boundary layer thickness (dimensionless) 

DEL2H, DEL2T 

^2H’ 

'^2T 

9 

Boundary layer thickness (dimensionless) 

HI, H2 

Hi, 

«2 

9 

Shape factor (dimensionless) 

RETH 



9 

Reynolds number based on momentum 


thickness 
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Subroutine EWL0SS (Arg. List) (Cont'd) 


Internal Variables 


THETIH, THETIT 

®1H’ ®1T 

, Momentum 

thickness 

(dimensionless) 

THET2H, THET2T 

®2H’ ®2T 

, Momentum 

thickness 

(dimensionless) 


Theory 


The interaction of the blade boundary layer (profile loss) with the endwall 
boundary layer is treated by adding an endwall loss to the blade force calculation 
procedure. This endwall loss is based on empirical correlations developed by 
Papaliou (Ref. 1) and is described in detail in Ref. 2. 

Ref erences 


1. Papailiou, K. D. , R. Plot, and J. Mathieu: Secondary Flows in Compressor Bladings, 

ASME J. of ENg. For Power, Vol. 99, April 1977, pp. 211-224. 

2. Anderson, 0. L., G. B. Hankins, D. E. Edwards: Extension to an Analysis of 

Turbulent Swirling Compressible Flow for Application to Axisymmetric Small Gas 
Turbine Ducts. UTRC Report R81-915395-12, NASA CR 165597, February 1982. 
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Subroutine FAVER2 


0b,1 ect 

Solve for mass flow weighted average flow 
Options 


None 



List of Symbols 


AAR 

= IB. 

, Area ratio (dimensionless) 

AMA 

= M 

, Local Mach number (dimensionless) 

AMACHE 

= M 

, Area average Mach number (dimensionless) 

AMFH 

= (Gft/V)H 

, Wall bleed hub (dimensionless) 

AMFT 

= (gm/v)t 

, Wall bleed tip (dimensionless) 

AMG 

= %ID 

, Midpoint Mach number (dimai sionless ) 

AMM 


, Maximum Mach number (dimensionless) 

ASH 

" -^SH 

, Surface area hub wall (dimensionless) 

AST 

= Ast 

, Surface area tip wall (dirdensionless ) 

ATFH 

-= (GQ/V)h 

, Heat flux/length hub (dimensionless 

ATFT 

= (gq/v)t 

, Heat flux/length tip (dimensionless 

BWORK 


, Work input (dimensionless) 

CPC 

= P^/^i 

, Normalizing factor for Cp (dimensionless) 

CPC0MP 

= Cpc 

, Pressure coefficient compressible (dimensionless) 

CPIWC 

= Cpi 

, Pressure coefficient incompressible (dimensionless) 

DASm.,DASH2 


, Area increment hub (dimensionless) 

DAST1,DAST2 

= AAgT 

, Area increment tip (dimensionless) 


11-142 



Subroutine FAVER2 (Cont'd) 


DENTPl, 

DENTP2 

= 

AI 

DISP 

= 

i 

DPS Til, 
DPSTI2 

= 

A0 

DQSHl, 

DQSH2 

= 

AQh 

DQSTl, 

DQST2 

= 

AQip 

DTHEOl, 

DTHE02 

= 

<] 

ENTP 

= 


PIO 

= 

"^0 

PI 00 

= 

^oo 

PSI 

= 


PSIO 

= 

^0 

m 

= 

1/2 (PU^)] 

ftSH 

= 


QST 

= 


R0M 

= 

Pmax 

RUM 

= 


THEO 

= 


THEOO 

= 

^ oo 

UM 

= 

^MAX 


, Change in Entropy (dimensionless) 

, Dissipation function (dimensionless) 

, Increment in mass flow (dimensionless) 

, Increment in heat flow hub (dimensionless) 

, Increment in heat flow tip (dimensionless) 

, Increment in total temperature (dimensionless) 

, Entropy (dimensionless 
, Average total pressure (dimensionless) 

, Average inlet total pressure (dimensionless) 

, Mass flow (dimensionless) 

, Initial mass flow (dimensionless) 

, Freestream dynamic pressure (dimensionless) 

, Total heat flow hub (dimensionless) 

, Total heat flow tip (dimensionless) 

, Freestream density (dimensionless) 

, Maximum momentum (dimensionless) 

, Average total temperature (dimensionless) 

, Inlet average total temperature (dimensionless) 
, Freestream velocity (dimensionless) 
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Subroutine FAVER2 (Cont'd) 


Theory 

The mass flow weighted average quantity 0 is defined by 

— I r' GP'-’s j. ^ „ 

^ 

where 

GPUg 

0 


The area is given by 

A = dn (3) 

*b 


It is noted that the mass flow weighted averages satisfy certain conditions 



The entropy is related to the dissipation by 


_d_ 

dS 


/G^n 

/ GI^I \ 

7 

r/GQ\ 

/GQ\ 1 

V V /j ^ 

V V Ih 

7-1 

LW®/t 

\ V® 'H j 


/ 


r r gq d® 
"I L dn 


7 


V®‘ 


dn 


+ yMr f ^^2 [^ns^ns ^r] 


dn 


( 6 ) 
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Subroutine FAVER2 (Cont'd) 


Eq. (^4) states that the mass flow is the net flow crossing the wall boundary. 

Eq. (5) states that the total energy flux is the net crossing the boundary walls 
Thus, for an adiabatic wall, the mass flow weighted average total temperature is 
constant. Finally, Eq. (6) states that the entropy flux is the net crossing the 
boundary plus the change due to the dissipation function and heat fluxes. Then 
using the definition of entropy we have 


rio2 

^01 



( 7 ) 


and the loss coefficient is given by 


C 


PL 



' - Poa/ Pq. 


I - P,/ P, 


01 


I - ngg/Hoi 

I- n I /rioi 


( 8 ) 


If we now defined the mass flow weighted averages 


n 


0 


r‘ n 



y-l 

2 



dn 


(9) 




0 



V 


v-yw M 



y-l 

2 


M dn 


(10) 


^ = arf 


7T 


0 V 


M dn 


(11) 
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Subroutine FAVER2 (Cont'd) 


and assume that the dependent variables have a small error E defined by 

fi'= (1+ 

©'= (i + e„)® 

m' =(l 

fio 


We may determine the' values of E to correct the solution to satisfy the mass flow 
weighted average equations. Thus, if we substitute Eq. (12) into Eq. (9) through 
Eq. (ll) and keep first order terms in E, 


where 


'® 


-n 


Cp 

Eu 

Eu 




^00 

E - E 

‘-no 

2 2 ^ 0 

En - E 

^ _E 

2 

Eu sin a 


& 


Eu^ = E u cos a 


-no 


n'o/n -1 




'0b ' ®0 /®o 


= 'i'H - 1 




) 


(13) 


(14) 


Thus the prime quantities are calculated From Eq. (9) through Eq. (ll) and the 
unprimed quantities from Eq. (4) through Eq. (6). 
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Sutroutine FAVER2 (Cont’d) 


The pressure coefficient is defined by 


Cp 


n - n, 

^01 ^ I 


( 15 ) 


An effective area A can be defined as the geometrical area minus the blockage caused 
by the boundary layer. If we define the freestream as the point of maximum velocity 
across the duct we have by definition 

= m = {S'j- Ng ( 16 ) 

Hence, the blockage B is defined as 

B = I - A/Ag^f = I - m/(yD^U^A) ( 17 ) 


The area averaged (effective Mach number) may then be defined by the isentropic 
flow relations. Thus, 


^eff 

A 


Meff 

M 


y-1 

I + M 



1 + 



1 y* ' 

2 y-1 


(l8) 


The effectiveness T] of the diffuser is based on an ideal isentropic flow with the 
mass flow weighted average Mach numbers. Thus, the ideal pressure coefficient is 
given by 


C 


PI 


n - n, 

no.-n, 


( 19 ) 


where 





( 20 ) 
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Subroutine FAVER2 (Cont'd) 



JIl 

M 


1 + 


1 111 

2 y-l 


Then 


( 21 ) 


T} " 


( 22 ) 


The wall friction coefficient is defined as 


Cf 



( 23 ) 


The wall surface area and heat flow are determined by integrating the equations 


A 


s 


277-r 
•'o V 


ds 


(24) 


0, 



27TRQ 

V 


ds 


( 25 ) 


using the trapezoid rule. 



Function FC0LES (ARG. LIST) 


Object 

Compute Coles' 
Options 

None 

velocity profile 

ARG. LIST 

AK = K 

+ 

, Von Karmen constant 

AKPL = K 

, Rougness Reynolds No. 

API = 

, Coles' shape factor 

APLS = A"*" 

, Van Driest constant 

BPL = B"*" 

, Law of wall constant 

DELT = 6 

, Boundary layer thickness (dimensionless) 

Y = y 

, Distance from wall (dimensionless) 

YPLS = y''' 

, Universal distance 

List of Symbols 

FC0LES U'*’ 

Theory 

, Universal velocity 


This subroutine constructs Coles' velocity profile (see Ref. 1). For > 20 
Coles' profile is valid and is given by 


U+ = ;^{[ln Y+ +B-*-]+2nc Sin^(fy/S)} ( 1 ) 

where the constant of integration for the law of the wall is obtained from the 
rough wall model function BPLUSR. 

For y"^ < 20, the velocity in the sublayer is obtained by integrating the 
equation (see Ref. 2) 
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dU~^ ^ 2 

dY"*' I + [l + (2 /c Y + D)^]''^^ (2) 


The damping factor D for a smooth wall is given by Van Driest (Ref. 2) 


D = I - exp (- Y'*'/A'*') 


(3) 


and for a rough wall it is given by McDonald (Ref. 3) 


D = I - exp (- Y'^/A'^) + 




(A) 


References 


1. Coles, D.: Law of the Wake in the Turbulent Boundary Layer. J. Fluid 

Mechanics, Vol. 1, 1956. 

2. Van Driest, E. F. : On Turbulent Flow Near a Wall. J. Aeronautical Sciences, 

Vol. 23, No. 11, 1956. 

3. McDonald, H., and R. W. Fish: Practical Calculations ' of Transitional Boundary 

Layers. Int. J. of Heat and Mass Transfer, Vol. 16, No. 9, 1973. 
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Subroutine FC0RCT 


0b.1 ect 

Correct solution based on mass flow average 
Options 
None 


List of 

Symbols 


AMUW 

V>v/^r 

, Wall value of viscosity (dimensionless) 

DU 

= AU 

, Velocity difference (dimensionless) 

EC0R 

= (pU)/(pU)m 

, Mass flux ratio (dimensionless) 

EMB 


, Error in Mach number (dimensionless) 

EPB 


, Error in static pressure (dimensionless) 

ERB 

= E 

, Error in density (dimensionless) 

ETB 

“ Erp 

, Error in static temperature (dimensionless) 

EUB 

= % 

, Error in velocity (dimensionless) 

EUP 

" ^U0 

, Error in swirl velocity (dimensionless) 

EUS 

"" ^US 

, Error in streamwise velocity (dimensionless) 

EW 

= E 

, Strain 

Theory 




Corrections are applied to the flow variables according to the theory outlined 
in Subroutine FAVER2. In addition, the point-to-point instability described in 
Ref. (l) is minimized by recalculating the stresses and heat flux using central dif- 
ferences rather than centered differences. The corrections from FAVER2 are currently 
set to zero. 

Reference 


1. Keller, H. B. : A New Difference Scheme for Parabolic Problems. Numerical Solu- 

tions of Partial Differential Equations. II SYNSPADE 1970, Academic Press, N.Y. 
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Subroutine FCPLX 


0b,1ect 

Evaluates Schwartz Christoff el transformation 
Options 

L0PT=1 Compute and store functions and derivatives 
=2 Compute only derivatives 

List of Symbols (Note subscript notation for derivatives used) 


N1,N2 

= %,N2 

XS 

= S 

XN 

= n 

XSX 


XSY 

II 

m 

II 

XL 

= D 

XZETS 

“ 5s “ ^n 

XETAS 

“ ''^Is ~"5n 

XXS 

- ^s - Yn 

XYS 

= Ys = 

XBl 

= X = = ’'ly 

YBl 

= ? = Tl* 

XI 

= X 

Y1 

= y 

XDB 

= D 

XNZET 


XNETA 

= n^ 

XV 

= V 
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Subroutine FCPLX (Cont'd) 


xsxx 

= S = n = s 
XX XX yy 

XSXY 

" ^xy " ^ " ^xx 

XXSS 

~ ^ss ^ ^ns ^nn 

XXSN 

= ^sn = ynn =-y3s 

XB2 

~ ^ ~ §xx “ "'lyx 5x; 

YB2 

= Y =-§xy ="^yy "=^xx 

XZSS 

~ ^ss “ "^isn “~^nn 

XZSN 

~ ^sn~”'^nn “""'^ss 

XDBN 


XDBS 


XESDN 

= (11s/D)n 

XESDN 

= (^s/D)s 

XZSDN 

~ (?s/^^n 

XZSDS 

= (5s/D)s 

XYW 

II 

XVS 

= ^s 

XESS 

^ss 

XESN 

~ '^sn 

XDIS 

= Ds 

XDIN 

Theory 

= Dn 


The theory for evaluating the complex functions and all derivatives is derived 
in Ref. (l). With the use of orthogonality relations which are implicit in the 
theory of complex function, the functions and derivatives may be evaluated. It is 
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Subroutine FCPLX (Cont'd) 


noted that this subroutine was programmed to accept multiple sources in the z plane, 
although they are not used in the present calculation. The derived functions 
calculated in this subroutine are listed as follows: 


At r T 

S = I In [(X-bi)^ + y^l 

1 = 1 


( 1 ) 


n = I Aj. ton’'[y/(x-bi)] 
1 = 1 


( 2 ) 


Sx = I 


NS Aj (x-bi) 


i = > [(x-bj)"+y^] 


(3) 


Sv = I 


NS Aj y 


1 = ' [(x-bjf + y"] 


(4) 


^xx ~ ^ 


2(X-bj) j 

[(x-bif+y^] ( [(x-bi)"+y']S 


(5) 


Syy = I 


NS Aj 2y(X-bi) 


I- [(x-b,)%y^l 


( 6 ) 


D = - (Sx + Sy ) 


(7) 
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Subroutine FCPLX (Cont'd) 




( 8 ) 


Ys = - Sy/D 


( 9 ) 


Cs ■ Cx Cy 


( 10 ) 


Vs ~ Vx. ^s"*" Vy'^s 


( 11 ) 



(12) 


= -[ss, 


(SyxXs 


^yy^s) 


+ 2Sx(SxxXs + 


SxyY; 


'] 


( 13 ) 


-[as 


y ( ^yx^n 


SyyYn) + 2 Sx { Sxx^n'*" ^xy'^n ) ] 


(l4) 


X 



^XX^S*^ SxyYg 

D 


SxD 



(15) 
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Subroutine FCPLX (Cont'd) 


^sn = L D D 2 -1 


^ss * ^ss ^xx ^^xy ^ss"*” ^yy^s 

«sn= f,>‘sn-^«y''sn+l«„-C.x'\>‘s + C.ylxl-''s> 

'^s ■ ~ '^ ^ [^s ^ss '*’ ^ns] ^19) 

~ ^ ^[^s ^sn"*" ^n ^nn] (20) 


Numerical accuracy can be significantly improved by ordering the way in which sums 
and products are made. As an example, the first equation, Eq. (7.2.29.1), may be 
written. 


= I 


X ^bi 

y 




tyl > |X-bj 


Thus the square root of the sum of squares of 0(l) and S = 0(/X-bj/). This rule has 
been applied to all equations by extracting the order of magnitude of the term from 
each calculation and products and quotients allowed to cancel order of magnitude 
terms . 
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Subroutine FETA (B, ETA, AN, DEDN, D2EDN) 


Object 

Calculate distorted mesh to be used in Subroutine P0IS 
Options 

B = 0 Uniform mesh (no stretching) 

B > 0 Tanh stretching 

Variables 


B 


, Constant 

ETA 

n 

, Transformed Normal Coordinate 

AN 

n 

, Normal Coordinate 

DEDN 

8n/ 3n 


D2EDN 

S^n/ 3n^ 


Theory 




This subroutine calculates the distorted mesh that will be used in calculation 
by subroutine P0IS. The transformation is given by 

n = tanh [Bn] /tanh [B] B 

n = n B 


> 0 
= 0 


( 1 ) 
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Subroutine FINTG(IKL) 


Object 


Integrate equations of Schwartz-Christof fel transformation 
Options 

IKL = Number of streamlines 


List of Symbols 

IKL , Number of streamlines 


Theory 


Four simultaneous ordinary differential equations are integrated using a third 
order Runge-Kutta numerical integration method. These equations from the Schwartz- 
Christoffel transformation (Ref. l) are denoted symbolically as. 


dx 

ds 

dY 

ds 

K 

ds 

drj 

~ds 


Xg(X,Y) 
Yg(X ,Y) 
$s^x ,Y) 


( 1 ) 


The Runge-Kutta formulas applied to the first equation are 


B„ = ASXg(x,Y) 

B,2 = ASXglx + Bii/Z ,Y+ Bg/Z) 

B,3 = ASXs(x + 2B,g-B„ , Y + 2Bgg- Bg,) 


X(S + AS) = X(S)+(B,,+ 4B,g+ B,3)/6 


( 3 ) 


Reference 


1. Anderson, 0. L.: User's Manual for a Finite Difference Calculation of Turbu- 

lent Swirling Compressible Flow in Axisymmetric Ducts with Slot Cooled Walls. 
Vols. I. and II, USAAMRDL-TR-T!+-50 , 19T^- 
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Subroutine FINVIS (L0PT) 


Object 


Calculate upstream/downstream blade force variables 
Options 


L0PT 

= 

1 

Calculate upstream flow variables 


= 

2 

Calculate downstream pressure 


= 

3 

Calculate downstream flow variables 

I0PT1O 

= 

0 

Stator 


= 

1 

Rotor 

IDBG7 

= 

0 

No debug printout 


ss 

1 

Debug printout 

List of Symbols 


A(I,K) 

= 

^I,K 

, Coefficients of differential equations 

EPS 

= 

e 

1 *v+l 

, Maximum error in 4i 

. '■V 1 

ERR 

&: 

\r 

1 ^v+l 

-ip 1 , Error in stream function 

. 1 

ERRBKP 



“U „ /V , Error in streamwise velocity 

s 2 s 2 

ITER 

= 

V 

, Iteration counter 

PSIHT 

= 

4' 

, Upstream mass flow (dimensionless) 

PSIMAX 

= 

'*’max 

, Maximum mass flow (dimensionless) 

PSIT 

= 

iv 

, vth guess for V (dimensionless) 

PSIl 

= 


, vth guess for upper bound (dimensionless) 

PSI2 

= 

2 

, vth guess for lower bound (dimensionless) 

XL 

= 


, Lower bound on x 

XM 

= 

^IN 

, X for choked flow 

XU 

= 


, Upper bound on x 

Y-1 

XX 

= 

X (p/p ) y , Pressure ratio 

o 

XI 

= 

K 

, vth guess for upper bound 

X2 

= 

X- 

, vth guess for lower bound 
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Subroutine FINVIS (Cont’d) 


Theory 

The blade force is calculated using the blade element theory described in 
Ref. (1). Subroutine FINVIS calculates the flow variables just upstream of the 
blades (Station JLE) and just downstream of the blades (Station JTE) . This data 
is used by subroutine F0RCE to calculate the force applied to the flow by the 
blades . 

L0PT_=_1_ ^aJiCul^t^ Up^t^e^m_F^ow 


The inviscid flow field without blades is calculated by subroutine CALINV 
and stored on unit NCALV. The location of the upstream and downstream stations 
(JLE, JTE) are calculated by subroutine SLETE. Data from the JLE(th) station 
is read from unit NCALV by subroutine F0RCE. All the flow field variables may 
then be calculated using the isentropic flow relations and the equation of state. 
Velocities in both the relative (rotor) and absolute (stator) are calculated. 
Finally a first guess for the downstream density and streamwise velocity is cal- 
culated for use by the backpressure iteration in subroutine F0RCE. 

L0PT = 2 Calculate Downstream Pressure 


Subroutine CASC which is called by subroutine F0RCE calculates the total pres- 
sure loss and downstream flow angle. For stators these are calculated in the 
absolute frame of reference and for rotors these are calculated in the relative 
frame of reference. The work is calculated using the rotor speed and the upstream 
and downstream tangential velocities in the Euler equation. For stators the work 
input is zero. For rotors the work input is not known until the downstream tan- 
gential velocity is determined. This requires the backpressure iteration loop in 
subroutine F0RCE. With total pressure, total temperature, and flow angle known, 
the downstream flow field is solved using the radial momentum equation and the 
global continuity equation with the algorithm described in subroutine CKINPT. 

L0PT = 3 Calculate Downstream Flow Variables 


When the downstream static pressure is determined the flow variables are 
calculated using the isentropic relations and the equation of state. 

Reference 

1. Barber, T. J., P. Raghuraman, and 0. L. Anderson: Evaluation of an Analysis 

for Axisymmetric Internal Flows in Turbomachinery Ducts. ASME Winter Meeting - 
Flow in Primary Nonrotating Passages in Turbomachines, p. 107, December 1979. 
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Subroutine FLINE 


Object 

Find coordinates of output data line 
Options 


IWALL 

= 1 OD Wall 

= 2 ID Wall 




M0LD 

= 1 Side 

(JX, 

KX)„ (JX+1 

, JX) 


2 Side 

(JX+1 

, KX), 

(JX+1, KX+1) 


3 Side 

(JX+1 

, KX+1) 

, (JX, KX+1) 


4 Side 

(JX, 

KX+1) , 

(JX 

, KX) 

Variables 






JBH, JBT 




9 

Wall intersection index (see Fig. 1) 

JX, KX 




9 

Corner point index (see Fig. 1) 

LX 




9 

Input data point index (see Fig. 1) 

M0LD 




9 

Old intersection index 

MNEW 




9 

New intersection index 

(RBI, ZBl) 

= (Ri» Zi) 



9 

Cartesian coordinates of point LX-1 (dimensionless) 

(RB2, ZB2) 

— (1^2’ ^2 



9 

Cartesian coordinates of point LX (dimensionless) 

(RXH, ZXH) 

(Rjj, Zj^) 



9 

ID wall intersection (dimensionless) 

(RXT, ZXT) 

— (R'pj Zrj,) 



9 

OD wall intersection (dimensionless) 

(SBl, ANBl) 

= (Sj, iij^) 



9 

Streamline coordinates of point LX-1 

(SB], ANB2) 

= (S2, n^) 



9 

Streamline coordinates of point LX 

SX, TX 

= S, T 



9 

Interpolation parameters for coordinate line 
and input data line 

XBL 

= ^1 



9 

Distance from inlet LX=1 

XBL 

" ^LLAST 



9 

Dlatance from inlet LX=LLAST 

XB2 

= ^LX 



9 

Distance from inlet LX=LX 

XXH, XXT 

= Xjj. x^ 



9 

See Fig. 1 
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Subroutine FLINE (Cont’d) 


Theory 

A pair of points on the coordinate grid is selected according to the index 
MNEW and a straight line in parametric form is constructed with S as a parameter. 
Then a pair of points (LX, LX-1) is selected and a straight line is constructed 
in parametric form with T as a parameter. If the two lines intersect within the 
two pairs of points, then the input data line intersects the grid line. This cal- 
culation is done with subroutine CR0SS1. 

The search procedure is initiated by finding the intersection of the input 
data line with the walls using subroutine ALINE. Then starting at the ID wall, 
the three sides of the coordinate rectangle are searched using subroutine CR0SS2 
with the point (JX, KX) as the corner point (side M0LD is skipped) and the inter- 
section with side MNEW determined. Then the neighboring coordinate rectangle is 
searched (see Fig. 8). This procedure continues until a coordinate rectangle is 
found that contains the input data point LX. 

When a coordinate rectangle contains the point LX, the streamline coordinates 
(s 2 » 112 ) calculated by interpolation using subroutine CLINE. Subroutine 

DLINE stores the appropriate data.- 
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JX = JBT 




Z 


Fig. 8. intersection of Input Data Line with Coordinate Grid 
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Subroutine FL0WIN 


Obj ect 

Set up Initial flow conditions 
Options 

IDBG3 = 0 
= 1 

10PT1 = 3, 4, 7, 8, 9 

I0PT11 = 0 
= 1 

I0TP12 =0, 1, 2, 3 

N0PT4 = 1 
= 2 

N0PT5 = 0 
> 0 

N0PT18 = 0 
= 1 

N0PT19 = 1 
> 1 

DR0UGI = 0 
> 0 

RHS(10),RTS(10) = 0 
> 0 

List of Symbols 

AKAPAC = K 

AMACHE = M 

AMACHl = ii , 

AMIH, AMIT = Mg , 

AMUIH, AMUIT = » 


No debug 
Print debug 

Inlet flow option (see Section 4.2) 

Duct flow 
External flow 

Turbulence model option (see Section 4.2) 

ID wall boundary condition 

ID axis of symmetry boundary condition 

Continue calculating 

Fatal error - stop calculation 

Flag - turbulent flow 
- laminar flow 

Print DIAGN0STIC's 
Do not print DIAGN0STIC's 

Smooth wall 
Rough wall 

Adiabatic wall boundary condition 
Wall temperature boundary condition 


Prandtl constant 
Area average Mach number 
Mass average Mach number 
Edge Mach number (hub, tip) 
Edge viscosity (hub, tip) 
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Subroutine FL0WIN (Cont’d) 


AMUlffi, AMUWT 

= 


Wall viscosity (hub, tip) 

ANH, AWT 

= 

n , 

Power law 

APLUSC 

= 


Van Driest constant 

AKPLUS 



Rough wall Reynolds number 

APRESl 

= 

P 

Mass average static pressure (psf) 

APRSOl 

= 


Mass average total pressure (psf) 

ATEMP 

= 


Mass average static temperature (*R) 

A2CF 

= 

Cf/2 

Wall friction coefficient 

BBO 

= 

Bo 

Blockage 

BPLS 

= 


Constant in law of wall 

CFRS 

= 

^fR'^^fS 

Ratio rough/smooth wall friction 

DIM, DIT 

= 


Boundary layer thickness (hub, tip) 

DSH, DST 

= 

6^/rr 

Displacement thickness (hub, tip) 

DR0UGH 

= 

d/r , 

r * 

Equivalent sand roughness height 

H12H, H12T 

= 

Hi2 

Shape factor (hub, tip) 

IBL 


9 

Flag - inner loop iteration 

ITDEL 


9 

Inner loop iteration counter 

ITER 


9 

Outer loop iteration counter 

KH, KT 


9 

Index boundary layer edge (hub, tip) 

PIH, PIT 

= 

He 

Cole's shape parameter (hub, tip) 

RIH, RIT 

= 

Pe/Pr 

Edge density (hub, tip) 

R0WH, R0WT 

= 

Pvj/Pr ’ 

Wall density (hub, tip) 

RETH, RETT 

= 

Ree 

Reynolds numbers - momentum thickness 

SIGWH, SIGWT 

= 

» 

Wall stress (dimensionless) 

TIH, TIT 

= 

» 

Edge temperature (hub, tip) (dimensionless) 

TWH, TWT 

= 

> 

Wall temperature (hub, tip) (dimensionless) 

UIH, UIT 

= 

“e 

Edge velocity (hub, tip) (ft/sec) 

USH, UST 

= 

* 

u , 

Friction velocity (hub, tip) (ft/sec) 

YABS 

= 

Y , 

Distance from wall (dimensionless) 
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Subroutine FL0WIN (Cont'd) 


YPLUS = Y"'” , Universal distance 

GAMMA = Y » Ratio of specific heats 

Theory 

The general analysis is shown by the flow chart in Fig. 9 and takes place 
in the following steps: 

1. Calculate total pressure P^, total temperature Tj, swirl angle a, and an 
initial guess for the static pressure P from the input data according to the 
inlet flow option I0PT1. 

2. Calculate Mach number M, static temperature P, density p, and velocity com- 
ponents from the input data and the isentropic flow relations given by 

y 


^ = (1 

(1) 

'To = 1 + ^ 

(2) 

C = (y 

(3) 

C 

II 

2 

o 

(4) 

Ug = u cos a 

(5) 

= u sin a 

(6.) 


3. Find edge of boundary layer. On the first iteration the edge of the turbu- 
lent boundary layer is estimated by the displacement thickness and power law 
(Ref. 1). 


8 = (l +n) S* 


( 7 ) 


and for laminar flow using the Blasius profile (Ref. 1) 
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Subroutine FL0WIN (Cont'd) 


( 8 ) 


On successive iterations the edge of the boundary layer is estimated from Cole's 
profile (Ref. 2) as shown in Step 4. When 6 is known, flow variables at the edge of 
the boundary layer may be determined by interpolation of the core flow variables. 

4. Calculate boundary layer profiles. With the edge conditions known, the 

boundary layer velocity profile is constructed from the Blasius solution using 
a linear interpolation in function UBLAS. For turbulent flow the velocity 
profile is constructed from Cole's profile, for which the following relations 
hold: 





+ 2n, 


+ B 


Sin 


+ 2 


(fi)} 

"c} 



ue S* 

X TTT ~T~ = I + n. 




Ue 


(1 + Dc) - (2. + 3.2 Tic + 1-55 nc ) 


(9) 

( 10 ) 
(11) 
( 12 ) 


H,2 = B*/B 


(13) 


Equation (11) may be substituted into (10) to obtain 




+ B +2 



(14) 


With Ug, p^, and a guess for 6, Eq. (14) may be solved for U using function 
CFC0LE. Equation (11) may be substituted in Eq. (12) to obtain 
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Subroutine FL0WIN (Cont'd) 




(I + nj - [2 + 3.2 



( 15 ) 


where 


Up S* 

- I (16) 

and Eq. (15) may be solved for Cole's shape factor, Eg. Then the parameters 
for constructing Cole's velocity profile (Eq. 9) are known. A new guess for the 
edge of the boundary layer is obtained from Eq. (11). 

6. Check inner loop for convergence 




/s <<t> 


(17) 


7. Calculate new static pressure. The normal momentum equation (Ref. 3) is inte- 
grated in subroutine ERPIN. 

8. Check outer loop for convergence 


lp>^ + ' -pn/p‘' < 4 > 


(18) 


When the static pressure converges, the subroutine also returns the weight flow 
WFL0 for use in the weight flow iteration by subroutine WFITER. 

References 


1. Schlicting, H. : Boundary Layer Theory. 6th Ed., McGraw-Hill, New York, 1968. 

2. Coles, D. E. : The Law of the Wake in the Turbulent Boundary Layer. J. Fluid 

Mech, Vol. 1, 1956, p. 191. 

3 . Anderson, 0. L. : Calculation of Internal Viscous Flows in Axisymmetric Ducts 

at Moderate to High Reynolds Numbers. Int. J. Computers and Fluids, Vol. 8, 
No. 4, p. 391, December 1980. 
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Fig. 9. Flow Chart for Subroutine Flowin 
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Subroutine FN0RM 


Ob j ec t 

Normalize input variables 
Options 


None 




List of 

Symbols 



AM 

= 

P 

* 

Viscosity (slug/ft/sec) 

AMACH 

= 

M 

9 

Mach number 

APO 

= 

Po 

9 

Total pressure (psfa) 

ATO 

= 

^o 

9 

Total temperature (deg R) 

BELT 

= 

e 

9 

Momentum thickness (ft) 

BTHE 

= 

e/h 

9 

Momentum thickness/duct height 

ETA 

= 


9 

n coordinate at 6 

HI 2 

= 

%2 

9 

Shape factor 

PHI 

= 




RETH 

= 

Re0 

9 

Reynolds number 

R0U 

= 


9 

2 

Freestream mass flow/area (slug/ft /sec) 

TEST 

= 

Pt 

9 

Test on p coordinate 

XI 

= 

E 

9 

Roberts transformation parameter 

YPLUSl 

= 


9 

Universal distance to first mesh point 

ZETA 

= 


9 

? coordinate at 6 



6 

9 

Boundary layer thickness (ft) 



* 

6 

9 

Bisplacement thickness (ft) 



h 

9 

Buct height (ft) 


Theory 

This subroutine performs the following tasks: 

1. Ivlormalize all input variables list of variables in Section 8.0 

2. Set gas properties to default values if not specified 

3. Selects the mesh distortion parameter DBS if not specified. 
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Subroutine FN0RM (Cont'd) 


The mesh distortion parameter is selected in the following manner. From 
Ref. (1), the Roberts transformation can be written. 


„ . (c + l/2) -(c-l/2) 


With c given and n specified at uniformly distributed mesh points, n may be 
determined. The parameter c can be related to the mesh distortion at the boundary 
Thus , 


( 1 ) 

( 2 ) 


_j dn I 


4 >{ 0 ) 
i + <^ (0) 


2 


(c +1/2) in(- 


c + 1/2 \ 
c-l/2/ 


(3) 


For turbulent boundary layers, the first mesh point from the wall is set at y"*" = 1 
Hence , 


Y 


+ 


^ Vi 
u 



(4) 


where 


0 = S / H,2 


(5) 




(. 6 ) 


For laminar flow, 1/3 of the mesh points are placed in the boundary layer. Hence, 


^8 = 3 SVh 

= 0.3 


(7) 


Equation (7) may be substituted into Eqs. (11) and (2) to solve for c and into 
Eq. (3) to solve for DBS. Let 
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Subroutine FN0RM (Cont'd) 


C = 1/2 + € 


where e is generally a small number. Then 




b = l/(277g-l) 


e = \/{ 4 >° -I) 


DDS 


=[TT77nrjT-(W 


These equations have a minimum distortion DDS. Let 


n-r = (2-2'^^-l)/(l +2 


J/b, 


If n, > n„, e = 1. 
6 T’ 


( 8 ) 


(9) 

( 10 ) 
( 11 ) 
( 12 ) 


(13) 
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Subroutine F0RCE (L0PT) 


Object 

Calculate blade force/span 


Options 


LOPT = 1 
= 2 

I0PT2 = 1 
2 

3 

4 

I0PT6 = 0 
= 1 

I0PT1O = 0 
= 1 

N0PT5 = 0 
> 0 


Calculate blade force 
Print variables 

Calculate blade force from input data 
Not in use 

Calculate blade force from cascade data 

Calculate force from input loss and air angle 

Blade force is calculated 

Thickness effects only 

Stator 

Rotor 

Continue calculation 
Error - Stop calculation 


List of Symbols 
ALIND = a. 


, Induced angle of attack (deg.) 


ALPHS 

= “s 

9 

Stagger angle (deg.) 



AMACHI 

= M. 
1 

9 

Absolute Mach number 

(stator) 



9 

Relative Mach number 

(rotor) 


ANGI, ANG2= 

9 

Absolute flow angles 

(stator) 

(deg. ) 


^2 

9 

Relative flow angles 

(rotor) 

(deg.) 

CH0RD 

II 

n 

9 

Chord 



CX21 

- ”s2«sl 

9 

Streamwise velocity ratio 


ERRBKN 1 


9 

Error 



ITER 

= y 

9 

Back pressure iteration counter 

JCL 


9 

Centerline station 



JLE 


9 

Leading edge station 



JTE 


9 

Trailing edge station 
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Subroutine F0RCE (Cont'd) 


0MEGZ 

= r^n/uy 

, Rotor speed 

PHIC 


, Camber angle (deg) 

RH0CX 

= (PU32)/(PU,i) 

, Mass flow ratio 

SOLD 

= 0 

, Solidity 

THIKM 

= t/c 

, Thickness/chord 

THIKN 

= t/r^ 

, Thickness 

WIND 

= Wi 

, Induced velocity (dimensionless) 

Theory 




The analysis for calculating the blade force is described in Ref. 1. Three 
input options are in use. 

I^0PT^ = ^ _Cal£ul^a_t e_F^r£e_Fr^om ^n£ujt Da£a_ 

If the flow variables upstream and downstream of the blade row are known, the 
blade force can be calculated directly. 

^^_PT2^ ^ 2 ^ _C^^£'^l.a£e_Fo^r£e_F_r om ^a£C£d£ Da£a_ 


The flow just upstream of the blade is calculated by subroutine CALINV and stored 
on unit NCALV. This data is read from unit NCALV using subroutine BLKRED and used 
by subroutine FINVIS(l) to set up all the upstream flow variables. Subroutine CASC 
calculates the downstream flow angle and loss coefficient. Subroutine FINVIS(2) 
calculates the downstream static pressure, and FINVIS (3) calculates all the remaining 
flow variables. For stators (I0PT1O = 0) a back pressure iteration is not required 
if the RH0CX corrections are not used in subroutine CASC. 

^ _C£l£,u^a£e_F£r£e_Inpjat_L£ss, £n^ 

The loss coefficient and exit flow angle may be obtained from input data cards 
replacing the output from subroutine CASC. The remaining algorithm is the same as 
I0PT2 =3. 

I0PT = 2 Print Solution 


This code is self explanatory. 
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Subroutine F0RCE (Cont'd) 


Reference 


1. Barber, T. J. , P. Raghuraman, and 0. L. Anderson: Evaluation of an Analysis 

for Axisymmetric Internal Flows in Turbomachinery Ducts. ASME Winter Meeting - 
Flow in Primary Nonrotating Passages in Turbomachinery. p. 107, December 1979. 
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Subroutine F0RCL 


Object 

Compute blade force/volume 
Options 

I0PT6 = 0 Use blade force 

= 1 Use thickness effects only 
List of Symbols 

GAP = G , Gap (dimensionless) 

ZLE “ Z , Location of leading edge (dimensionless) 

LE 

ZTE = Z , Location of trailing edge (dimensionless) 

TE 

Theory 

The force per unit span is calculated by subroutine F0RCE. This subroutine 
calculates the force per unit volume by dividing by gap and projected chord length 
according to Ref. 1. 

Reference 


1. Anderson, 0. L. : User's Manual for a Finite Difference Calculation of Turbu- 

lent Swirling Compressible Flow in Axisymmetric Ducts with Struts and Slot 
Cooled Walls. USAAMRDL-TR-74-50, Vol. 1, 1974. 
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Subroutine GBLADE 


Object 

Interpolate blade parameters 

Options 

I0PT2 = 0 
> 0 

L0PT2 = 0 
= 1 

List of Symbols 

See C0MM0N /AC0NS0/ 

Theory 

This subroutine linearly interpolates the blade parameters 


No blade force 
Blade force see Page 1-13 
Calculate blade blockage 
Do not calculate blade blockage 
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Subroutine F0UTP 


Object 


Interpolation solution on output data line 
Options 


NLAST = 0 , 

Return 



> 0 

Interpolate Solution 

List of 

Symbols 



EPS 

E 

9 

Small increment in Y 

FIl(I) 

ffCn, Sj^l) 

9 

n interpolated solution at J+1 (dimensionless) 

FI2(I) 

fl(n,Sj+2) 

9 

n interpolated solution at J+2 (dimensionless) 

FI(I,L) 

f j(n,s) 

9 

Interpolation solution at n,s (dimensionless) 

GAP 

G 

9 

Blade gap (dimensionless) 

Jxl, Jx2 

Kxl ,kx2 

9 

Four corners of interpolation grid 

NLAST 

N 

9 

Number of points on data line 

THTS(K) 

0 

9 

Angle between streamline and Z axis (rad.) 

RHP , ZHP 


9 

ID coordinate of output data line (dimensionless) 

RTP,ZTP 


9 

OD coordinate of output data line (dimensionless) 

RWH, ZWH 

^>^WH 

9 

Intersection with ID wall 

RWT, ZWT 

^WT’^T 

9 

Intersection with OD wall 

Theory 





The coordinates of the end points of the output data line (RHP,ZHP,RTP,ZTP) 
are read and normalized. These points must lie outside the duct to resolve 
ambiguities in the interpolation procedure. The intersection of the output 
data line with the ID wall (RWH, ZWH) and the OD wall (RWT, ZWT) is determined 
using subroutine ALINE. The output data line between the ID and OD wall 
intersections is divided into NLAST-3 intervals for NLAST data points which 
include the end points. Since the second and next to last point must be 
inside the duct, the points (RWH, ZWH, RWT, ZWT) are moved by an amount E. 

The location of the NLAST points in the (n,s) coordinates is determined using 
subroutine FLINE. The output data point is then located within the rectangular 
grid bounded by (J+1, K+1; J+2, K+1; J+1, K+2; J+2, K+2) and the solution can be 
calculated using linear interpolation. 
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Function FTHIK (Z,IS,L0P) 


Obj ect . 

Calculate blade thickness 


Option 



L0P 

= 1 

FTHIK = t/c 


= 2 

FTHIK = d(t/c)/d(X/c) 

IS 

= 1 

NACA 5 digit series airfoil 


= 2 

Parabolic airfoil 


= 3 

NACA 5 digit series airfoil 


= 4 

Input thickness distribution 


= 5 

NACA 65A series airfoil- 

Theory 




The NACA series airfoil shapes are given in Ref. 1 and selected by the blade 
shape parameter IS. 

Reference 

1. Abbot, I. H. and A. E. Von Doenhoff: Theory of Wing Sections, Dover Pub., 1959. 
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Subroutine GDUCT 


Object 

Calculate duct shape 


Options 

I0PT3 


= 1 
= 2 
= 3 


Straight annular duct 

Arbitrary duct with evenly spaced axial stations 
Straight wall annular diffuser 


= 4 

Duct 

shape stored on file JDRUM 

= 5 

Arbitrary duct with arbitrary stations 

N0PT4 = 1 

No ID 

wall 

= 2 

ID wall 

List of Symbols 
AREAR 

2 

= r 

r 

2 

, Reference area (ft ) 

DUCTI (N) 


, Duct input parameters 

JL 


, Number of streamwise stations 

JLPTS 


, Number of input points 

KN0TS 


, Number of nodes for spline fit 

r(i,i,J),R(i.a,J) 

= 

, Output coordinates OD wall (dimensionless) 

R(2,1,J),R(2,4,J) 


, Output coordinates ID wall (dimensionless) 

RDII(J), ZDII(J) 

= 

, Input coordinates of OD wall (ft) 

RD2I(J), ZD2I(J) 


, Input coordinates of ID wall (ft) 

RADR 

= r 

r 

, Reference radius (ft) 


Z1 


, Duct axial length (ft) 
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Subroutine GDUCT (Cont'd) 


Theory 

For I0PT3 = 1 or 3, the duct shape is calculated for JL equally spaced axial 
stations using algebraic relations with DUCTI(N) as parameters. For I0PT3 = 4, 
the geometry and coordinates are read from data file JDRUM. For I0PT3 = 2 or 5, 
the least squares spline fitting and smoothing routine SM00TH is used. For I0PT3 = 2, 
the independent variable is the axial distance with evenly spaced axial stations. 

For I0PT3 = 5, the independent variable is the distance along the wall which may be 
arbitrarily spaced. 
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Subroutine IKIT 


Obj ect 

Set all C0I-1M0N block variables to zero 
Options 


None 

Theory 


This subroutine sets all C0MM0N block variables to zero which are defined in 
Section 6.2 in the Manual. Subroutine IZER0 sets all integer variable and subroutine 
ZER0 sets all real variables. 
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Subroutine INITF 



Variables 

FIV (1, L, K) ij; , Stream function (dimensionless) 

FIV (2, L, K) Ug, Streamwise Velocity (dimensionless) 

FIV (3, L, K) U^, Tangential Velocity (dimensionless) 

FIV (4, L, K) n , Static Pressure (dimensionless) 

FIV (5, L, K) I , Entropy (dimensionless) 

FIV (6, L, K) 0 , Static Temperature (dimensionless) 

FIV (7, L, K) p , Density (dimensionless) 

FIV (8, L, K) 

FIV (9, L, K) 

FIV (10, L, K) 

L =10 J-1 station 

L = 2 0 J station 

L =30 J-M station 


K = 1, KL streamlines 

3 

FIPARM (1) , Reference density, (slug/ft ) 

FIPARM (2) T^, Reference temperature, (“Rankin) 

FIPARM (3) P^, Reference pressure* (psf abs) 
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Subroutine INITF (Cont'd) 


Variables (Cont'd) 

FIPARM (4) = g , Gravitational constant (ft/sec) 

FIPARM (5) = , Reference viscosity (slug/ft/sec) 

2 2 

FIPARM (6) = C , Specific heat constant pressure (ft /sec /°R) 

P 

2 2 

FIPARM (7) = , Specific heat constant volume (ft /sec /°R) 

2 2 

FIPARM (8) = R , Gas constant (ft /sec /°R) 


FIPARM (9) = P^ , Prandtl number turbulent 

FIPARM (10) = Not used 

FIPARM (11) = u^ , Reference velcoity (ft/sec) 

FIPARM (12) = N0PT7 , Number of blocks stored inviscid solution 

FIPARM (13) = , Not used 

FIPARM (14) = I0PT15 , First station 

FIPARM (15) = I0PT16 , Last station 
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Subroutine INITFV 



Variables 

FIV (1, L, K) ij; , Stream function (dimensionless) 

FIV (2, L, K) Ug, Streamwise Velocity (dimensionless) 

FIV (3, L, K) U^, Tangential Velocity (dimensionless) 

FIV (A, L, K) II , Static Pressure (dimensionless) 

FIV (5, L, K) I , Entropy (dimensionless) 

FIV (6, L, K) 0 , Static Temperature (dimensionless) 

FIV (7, L, K) p , Density (dimensionless) 

FIV (8, L, K) 

FIV (9, L, K) 

FIV (10, L, K) 

L =10 J-1 station 

L = 2 0 J station 

L =30 J+l station 

K = 1, KL streamlines 

3 

FIPARM (1) , Reference density, (slug/ft ) 

FIPARM (2) T^, Reference temperature, (‘Rankin) 

FIPARM (3) P^, Reference pressure* (psf abs) 
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Subroutine INITFV (Cont'd) 


Variables (Cont'd) 


FIPARM 

(4) 

= g . 

Gravitational constant (ft/sec) 

FIPARM 

(5) 

' ^r * 

Reference viscosity (slug/ft/sec) 

FIPARM 

(6) 


2 2 

Specific heat constant pressure (ft /sec /*R) 

FIPARM 

(7) 


2 2 

Specific heat constant volume (ft /sec /*R) 

FIPARM 

(8) 

= R 

2 2 

Gas constant (ft /sec /*R) 

FIPARM 

(9) 

= Rv . 

rT 

Prandtl number turbulent 

FIPARM 

(10) 

= 

Not used 

FIPARM 

(11) 

= u , 

r 

Reference velcoity (ft/sec) 

FIPARM 

(12) 

= N0PT7, 

Number of blocks stored inviscid solution 

FIPARM 

(13) 

= N0PT8 , 

Number of blocks stored, viscous station 

FIPARM 

(14) = 

I0PT15 

First station 

FIPARM 

(15) = 

I0PT16 , 

Last station 

Theory 





This subroutine initializes the dependent variable array FIV which is stored 
on a disk file and sets all parameters FIPARM required by the calculation. 
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Subroutine INITQ 


Ob.ject 

Initialize data file parameters for Q array 
Option 


None 



Variables 



BLOCK (1) 

JGTEP 

, Block (record number) 

Q(1,K) 

R 

, Radius (dimensionless) 

Q(2,K) 

Z 

, Axial distance (dimensionless) 

Q(3,K) 

8R/3n 

, Derivative 

Q(A,K) 

9R/3S 

, Derivative 

Q(5,K) 

(cos6)^ . 
axi 

, Axisymmetric flow angle 

Q(6,K) 

V 

, Metric coefficient (dimensionless) 

Q(7,K) 

3V/3n = (K + AK ) 
s s 

, Curvature of streamline 

Q(8,K) 

3V/3s 

, Curvature of potential line 

Q(9,K) 

X 

, Distance along streamline (dimensionless) 

Q(IO.K) 

Y 

, Duct height (dimensionless) 

Q(11,K) 

Y/Y^ 

, Normalized duct height 

Q(12,K) 

A 

, Duct Area (dimensionless) 

Q(13,K) 

G 

, Gap (dimensionless) 

Q(14,K) 

3G/3n 

, Derivative 

Q(15,K) 

3G/3s 

, Derivative 

Q(16,K) 

3n/3n 

« Transformation of normal coordinate 
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Subroutine INITQ (Cont'd) 


Variables (Cont'd) 


Q(17,K) 


, Not used 

Q(18,K) 

n 

, Normal coordinate (dimensionlesi 

Q(19,K) 

n 

, Transformed coordinate 


K = 1, KL 


QPARM(l) 

r 

r 

, Reference radius, (ft) 

QPARM(2) 


, Not used 

QPARM(3) 

JL 

, Number of streamwise steps 

QPARM(4) 

KL 

, Number of streamlines 

Theory 




This subroutine initializes the independent variable array BLOCK which is 
stored on a disc file and sets all parameters QPARM required by the calculation. 

* Note that Q(7,K) stores either Kg or Kg+AKg. 
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Subroutine INTFRE 


Obj ect 

Initialize freestream conditions 
Options 

None 

List of Symbols 

See C0MM0N BL0CKS 
Theory 

Data read from file NDRUM is used to set up the freestream conditions for 
subroutine P01S. 
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Subroutine IZER0 (lADDR, NUM) 


Object 

Set integer variables to zero 
Options 


None 

List of Symbols 

lADDR Integer array name 

NUM Number of elements in array 

Theory 

This subroutine sets the first NUM elements of integer array lADDR to zero. 
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Subroutine MINVRT(A,B,N) 


Object 

Invert KxN Matrix 
Options 
None 

List of Symbols 


A 

= A 

, Augmented A matrix 

B 

= 1 

, Augmented B matrix (A“^) matrix 

N 


, Number of equations (rows) 


M , Number of columns 


Theory 

The A matrix is inverted using the Gauss-Jordin elimination procedure. First 
the augmented AN, M is formed including the identity matrix, 


A = (a I) 


( 1 ) 


Then the following revision formula is used 


b 




°l,J 


r l< I< N 1 
t l< J < M J 


( 2 ) 


b 


N, J- I 




1 < JM 


( 3 ) 


Note that the B matrix has one less column than the ^ matrix. Then the substitution 
is made 
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Subroutine MINVRT(A,B,N ) (Cont'd) 


°IJ = 


1 < T< N 
1 < J < M 


and then repeated until the B matrix is an NxN matrix or the A 


matrix. 
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Subroutine 0RTFUN 



List of 

Variables 


Q(J,K) 

= 


Coordinate functions at station 

G(1,K) 

= 



G(2,K) 

= 

[XY]k_i/2 


G(3,K) 

= 

[G(XV)]^_3^/2 


G(4,K) 

= 

3V1-J 

ou^K-1/2 


G(5,K) 

= 

fl aRi-J 

XR 3n K-1/2 


G(6,K) 

= 

rl 3 RtG 

XR 3S K-1/2 


G(7,K) 

= 

[L. (G)]'^ 

XG 3n V K-1/2 

G(8,K) 

= 

[V_ 3 (G)J^ 

XG 3n XY 

iX] 

9n 

G(9,K) 

= 

XG 3n XR 

3n k-1/2 

G(9+I,K) 

= 

G(I,K) 0 J-1, I 

=1,9 K=l, KL 

DF(1,K) 

= 


Streamwise blade force/volume 
(dimensionless) 

DF(2,K) 

= 

'"♦'“''k-1/2 

Tangential blade force/volume 
(dimensionless) 

DF(3,K) 

= 


Total pressure loss/volume 
(dimensionless) 

DF(4,K) 

= 

[X]-^ 

k-1/2 

Coordinate distortion 
(dimensionless) 

DF(5,K) 

= 


Streamwise blade force/volume 


(dimensionless) 
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Subroutine {SRTFUN (Cont'd) 


DF(6,K)’ 


Tangential blade force/volume 
(dimensionless) 

DF(7,K) 


Total pressure loss/volume 
(dimenionless) 

DF(8,K) 

■ 

Coordinate distortion 
(dimenionless) 


where 



X = dn/dn 



Theory 


The variables G and DF are calculated in this subroutine in order to be used 
in subroutine S0LVI. 
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Subroutine 0UTPUT 


Object 

Print title page 
Option 
None 

List of Symbols 
None 
Theory 

This subroutine prints the title page which records all modifications, dates, 
and references to changes incorporated into the ADD code.” 
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Subroutine P0IS (RESM.ITER) 


Object 

Solve Poisson equation 


Option 


IDBGP = 0 
= 1 
= 2 


No debug printout 
Printout residuals 
Print solution 


List of Symbols 


P(K,J) = 

>1^ 

, Stream function (dimensionless) 

F(K,J) = 


> Coefficient (1/PC) (dimensionless) 

PSI(K) 

'V 

4, 

, Iterative guess for J 

ITER 

V 

» Iteration counter 

RLX 


» Relaxation factor 

RESMAX = 

^MAX 

, Maximum residual accepted 

RESDM 


, Maximum residual/J station 

RESM 

e 

> Maximum residual/sweep 

Theory 




The solution algorithm is described in Reference 1. 

References 

1. Anderson, 0. L. and D. E. Edwards: Extension to an Analysis of Turbulent 

Swirling Compressible Flow in Axisymmetric Ducts, NASA Contract NAS3-21853, 
1981, UTRC Report R81-914720. 
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Subroutine P0ISCF 


Object 


(1) Set of coefficients of = 0 

(2) Set boundary conditions on tp 
(2) Set initial guess for 


Options 

IDBG17 


0 Compressible Flow 

1 Incompressible Flow 


IDBGT = 0 No debug test case 
= 1 Debug test case 


IDBGP 

=0 No debug test 
= 1 Debug printout 

printout 

Variables 



Q(I,K) 



, Coordinate functions 

FIV(I, 

L,K) 


, Dependent variables for inviscid flow 

P(K) 

= 

K 

.Stream function at station J 

A(K) 

= 

GP/V/ (dri/dn) 

.Coefficient for 3i|j/9ri 

F(K) 

= 

1/g/P 

.Coefficient of 

G(K) 

= 

V/g/P 

.Coefficient for velocity calculation 

R(K) 

= 

P 

.Density ratio (P/Pj.) 

T(K) 

= 

0 

.Temperature ratio (T/Tj.) 

V(K) 

= 

V 

.Metric coefficient dimensionless 

GAP 

= 

G 

.Gap (g/r^) 

VMET 

= 

V 

.Metric coefficient (dimensionless) 

RH0 

= 

P 

.Density ratio (p/p^) 

TEM 

= 

0 

^Temperature ratio (T/Tj.) 

uso 

= 

U 

so 

.Upstream constant velocity (uq/Uj.) 


11-197 



Subroutine P0ISCF (Cont'd) 


USINF 

U 

goo 

jFree 

stream axial velocity (’^g„o/^r) 

UPINF 

(j)<» 

,Free 

stream tangential velocity U(j,^/Uj. 

PSIKL 


jFree 

stream stream function (dimensionless) 

TEMINF 

0a, 

jFree 

stream static temperature ratio (T/Tj.) 

RHOINF 


jFree 

stream density ratio (P/p^) 

AMINF 


,Free 

stream Mach number 

PTINF 


jFree 

stream total pressure ratio (Pq/Pj-) 

BLK 


,See COMMON/SPCGD/ 

Theory 





This subroutine does the following steps 

(1) Reads coordinate Q file and solution FIV file 

(2) Interpolates the solution to the (n»S) grid 

(3) Calculates the coefficient F = 1/PG 

(A) Calculates coefficients BLK for streamline curvature calculation 

(5) Sets boundary condition on Y 

(6) Calculates initial guess for 'i‘ 

(7) Stores F, BLK, P on disk files 

The initial guess is given by the inviscid solution obtained from CALINV. 
boundary conditions are given by: 


V'(0,s) = 0 


V|/(i7,Sl) = Uso 


The 
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Subroutine P0IS0N 


Object 

Calculate axisymmetric streanline curvature 
Options 


I0PT7 = 0 

= 1 
IDBG15 = 0 

> 0 

List of Symbols 


No curvature corrections 
Curvature correction 
Use input KL streamlines 
Use IDBG15 streamlines 


IRH0 

, Density 

iteration counter 

ITERL 

, Maximum 

number of iterations 

KH0LD 

,No. ADD 

code streamlines 

KL 

,No. SCURVA streamlines 

RESMAX 

, Maximum 

residual for convergence 


Theory 

This subroutine is a calling subroutine for subroutines INTFRE, P0ISCF, 
P0IS, and SCURVA 
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Subroutine PRT0UT 


Object 

Calculate and print certain terms neglected in ADD code equations 


Options 


I0PT18 = 0 
= 1 

Do not 
Print 

print terms 
terms 


I0PTA 

Print 

every I0PT4 

station 

List of Symbols 




JFIRST 




JLAST 




KL 




DS 


As 


DETA 


An 


REY 


Re 

r 


F 




Q 




VARBL (1, J, K) 




VARBL (2, J, K) 


U 

s 


VARBL (3, J, K) 




VARBL (4, J, K) 


P 


VARBL (5, J, K) 


Z 

ns 


VARBL (6, J, K) 


^E 

Pr 


VARBL (7, J, K) 


R 


VARBL (8, J, K) 


Z 


VARBL (9, J, K) 


V 



First calculated station 
Last calculated station 
Number of streamlines 
Stepsize in streamwise 

Stepsize in transformed normal direction 
Reference Reynolds number 
Flow array 
Geometry array 

Stream Function (dimensionless) 
Streamwise Velocity (dimensionless) 

Swirl Velocity (dimensionless) 

Density (dimensionless) 

Streamwise stress (dimensionless) 
Effective Viscosity (dimensionless) 

Radius (dimensionless) 

Axial Distance (dimensionless) 

Metric Scale Coefficient (dimensionless) 
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Subroutine PRT0UT (Cont'd) 


VARBL (10, 

J, 

K) = G 

jGap (dimensionless) 

VARBL (11, 

J, 

K) = f 

, Streamwise Derivative of radius 

VARBL (12, 

J, 

^ 

» Curvature of potential line 

VARBL (13, 

J, 

K) = -U 

> Curvature of stream line 

VARBL (14, 

J, 

K) = y/h 

’Fractional distance across duct 

Output Variables 


Q(ll, K) 


(Y/Yt'' 

Fractional distance across duct 





UDUDS 




UN2V 


P t 


UUNV 

UDUNDS 


P U. U_/V 4^ 
in (Jn 

(3Un 


DUN2DN 


P u ^ 

dn 


UDUPDS 


_ 


UNUPR 


p Un U^/R § 

DSIGNS 


V/G ^ (G/V 

^ns ) 

DSIGNP 


> 

o 

> 

2s<^) 

UN 


-V/GP|f 


DSIGSS 


V/G ^ I'G/V 

^ss) 

Theory 





In the formation of the ADD code equations from the Navier-Stokes equations 
certain terms were assumed negligible. This subroutine calculates some of the 
neglected terms. 
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Subroutine QINTER 


Objec t 

Interpolate curvature from P(^IS mesh to S^LVI mesh 
Options 


None 

List of Symbols 


Q(J,K) Coordinate functions 

Theory 

After the curvature and flow angle has been calculated from the potential 
flow solution, this subroutine interpolates to obtain values at the numerical 
grid points which will be used in the S0LVI calculation. 
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Subroutine QPTMAX (QMAX,PTMAX,RUMAX) 


Object 

Calculate maximum total and dynamic pressure 
Options 
None 

List of Symbols 

PTMAX P^^ 

QMAX 

RUMAX puj. 

Theory 

At any given streamwise station, this subroutine searches through the 
solution for the maximum velocity. The total and dynamic pressure are calculated 
at the point across the duct where the velocity is a maximum. 


Maximum total pressure (psfa) 

Maximum dynamic pressure (psfa) 

2 

Maximum mass flow rate (slug/ft /sec) 
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Subroutine RDC0R(J) Coordinate Block 


Object 

Read J and J+1 
Options 
None 

List of Symbols 

J , Block number 

Theory 

Coordinate data is read into core using BLKRED. 
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Subroutine READPF(J,JJ) 


Object 

Read P and F files in NIST word blocks 
Option 

None 

List of Symbols 


J 

JJ 

N 

F 

P 

NST 

NBK 

NIST 

NFDRM 

NL 

NBIST 

NMOVE 


J 

JJ 

N 

F 

= 25 


= 23 


, Record number 
, Record number in block N 
.Block number 

2 

.Coefficients of V = 0 
.Stream function (dimensionless) 
.Number of records per block 
.Number of words to move pointer 
.Number of words per block 
.Coefficient file number F array 
.Last block number 
.Number of words for two records 
.Number of words to move pointer 


Theory 


The entire F and P arrays cannot be kept in core at the same time so that the 
I/O is arranged to keep fixed blocks in core (Fig. 1). Let (J.K) be a point on 
the computational mesh and (JJ.KK) the corresponding point in core. Let each 
record be the Jth line with the number of words in the record given by 

K = l.IST 

If there are NST records per block, then these are NIST words per block. 


NIST = NST X 1ST 
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Subroutine READPF(J,JJ) (Cont'd) 


Theory (Cont'd) 

The solution algorithm requires overlapping blocks as shown on Fig. 10. 
Hence we have the block number 

N = (J-2)/(NST-2) + 1 
and the JJ point in core is given by 

JJ + J-(N-l) X (NST-2) 

This subroutine is coded so that a new block is ready only when N = NL. 
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FIG. 10 I/O for Streamline Curvature Solution 
Overlapping Blocks 


Record No. 
J 

1 

2 

3 

4 

23 

24 

25 

26 

46 

47 

48 


70 

71 


93 

94 

95 

96 

97 

98 

99 
100 
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Subroutine READPG(J,JJ) 


Object 

Read variable for curvature calculation 
Options 
None 


List of 

Symbols 


J 

= J 

, Record number 

JJ 

= JJ 

, Record number in block N 

N 


, Block number 

G 


, Streamline coordinate data (see COMMON /SPCFD/) 

P 


.Stream function 

NST 

= 25 

.Number of records per block 

NBK 


.Number of words to move pointer 

NIST 


.Number of words per P block 

NGDRM 

= 25 

.Unit number of G array 

NPDRM 

= 24 

.Unit number for P array 

NL 


.Last P block number 

NBIST 


.No. words for 2 P records 

NGIST 


.No. words for 2 G records 

NGL 


.Last G record number 

Theory 




This subroutine reads the P file according to Subroutine READPF but 
reads only the Jth G record which is kept in core. 
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Subroutine REC0RD (ARG. LIST) 


Object 


Calculate record numbers 
Options 


I0PT14 = 0 , KDX = 1 


> 0 , KDX > 1 


JRECFI > N0PT7 
JRECFV > N0PT8 


List of Symbols 


J 

JRECFI 

JRECFV 

NFI 

NFV 

N0PT7 

N0PT8 


JRECFI = N0PT7 
JRECFV = N0PT8 


Jth coordinate station 

Corresponding inviscid flow record number 

Corresponding viscous flow record number 

Index inviscid flow 

Index viscous flow 

No. records on unit MDRUM 

No. records on unit NDRUM 


Theory 

See Section 5.5 
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Subroutine REDINP 


Object 

Read input parameters and variables 
Options 

Data cards are read according to input options 
I0PTN where N = 1,17. 

List of Symbols 

Common block variables 

Theory 

This subroutine reads all input variables 
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Function R0BRTS(C,ETA,L0P) 


Object 

Compute distorted mesh using Roberts' transformation 
Options 

L0P = 0 Wall - wall boundary 

= 1 Wall- free stream boundary 
=-l Free stream-wall boundary 

List of Symbols 


c 

= C 

.Distortion parameter 

ETA 

= n 

.Input variable (uniform mesh) 

L 0 P 

= 

.Option 


Output Variable 

R0BRTS = n .Output variable 
Theory 

The Roberts' transformation for a distorted mesh on both sides is given by 

(c+ 1/2) exp [2 in(|^^)(i7'-i/2)]-{c-l/2 
n' = ( 1 ) 

where 


0 < 17' < 1.0 0 < n‘ < 1.0 


For the different options we have 


V = •’7 

n = n‘ 


} L 0P = 0 


i = 7}/Z 
n = 2n' 



’q' = U-^-q)/2 
n = 2n' - I 



( 2 ) 

(3) 

(A) 

( 5 ) 
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Subroutine SCURVA (IDBU.KH0LD) 


Object 

Calculate curvature from potential flow solution 
Options 

IDBU = 0 Update density 

> 0 Print SCURVA solution, update curvature 

Variables 


KH0LD 


,No. ADD code streamlines 

KL 


,No. SCURVA streamlines 

Q(J,K) 


, Coordinate functions 

P(K) 

< 

, Stream function a station J (dimensionless) 

G(K) 

G 

, Coefficient for velocity calculation 

R(K) 

P 

, Density Ratio (p/p^) (dimensionless) 

US 

Us 

jStreamwise velocity (dimensionless) 

UN 

Un 

.Normal velocity (dimensionless) 

U 

U 

.Total velocity (dimensionless) 

COSTH 

cos^ (0) 

2 

. (Cosine) of flow angle 6 

CURV 

5V 

2n 

.Curvature of streamline (dimensionless) 


Theory 

Once the potential flow solution has been obtained from subroutine P0IS, 
this subroutine will calculate the flow angle and streamline curvature according 
to Ref. (1). 

References 

1. Anderson, O.L. and D. E. Edwards, Extension to an Analysis of Turbulent 
Swirling Compressible Flow in Asixymmetric Ducts, NASA Contract No. 
NAS3-21851, 1981, UTRC Report R81-914720. 
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Subroutine SLETE 


Ob j ec t 


Find indices for strut control surfaces 
Options 


If I0PT2 = 0, JLEDG = JTEDG = 101 


List of Symbols 


RHLE, 

ZHLE 

= 

^LE’ ^LE’ 

Coord inate 

RHTE, 

ZHTE 

= 

\te, ^HTE, 

Coordinate 

RHLE, 

ZTLE 

= 

%LE> ^TLE’ 

Coord inate 

RTTE, 

ZTTE 

= 

Rtte> Z^tE’ 

Coordinate 


Theory 


hub leading edge (dimensionless) 
of hub trailing edge (dimensionless) 
of tip leading edge (dimensionless) 
of tip trailing edge (dimensionless) 


The projection of the blade onto the (R,Z) plane is determined by calculating 
the coordinate points of the four corners of the blade using subroutines TRLETE, 
TRBLDE. Then a straight line is drawn connecting the leading edge points and a 
straight line connecting the trailing edge points is drawn. Subroutine ALINE is 
used to determine the intersections of the leading and trailing edges with the walls. 
For a unique solution, these four points must lie outside the duct. Then JLEDG is 
the station just upstream of the leading edge and JTEDG is the station just down- 
stream of the blade. 
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Subroutine SL0TA (L0P) 


Object 

Construct duct with slots 
Options 
ENTRY SL0TA 

LOP = 1 Read slot parameters 

= 2 Construct walls with all slots 

= 3 Normalize slot variables 

= 4 Arrange slots in order of S 

= 5 Write slot input data 

ENTRY SL0TGD (N0PT) 

Construct N0PT'th slot 


List of Symbols 


CSN 

= 

cos6„ 

DRSL0T 

= 

Hcos9^ 

DZ 

= 

AZ 

SHUB.STIP 

= 

\*^T 

ZHUB , ZTIP 

= 

^H’^T 

ZZHUB.ZZTZP 

= 

^ZH’^; 

ZZ 

= 

^SL 

Theory 




, Cosine of wall angle 

, Increment in slot height (dimensionless) 

, Increment in axial distance (dimensionless) 
, S coordinate location of slot 
, Axial location of walls (dimensionless) 

, Axial location of slots (dimensionless) 

, Axial location N0PT'th slot (dimensionless) 


A duct without slots is constructed on the first pass. Then on successive 
passes Hcos0^^ is added according to Fig. 11. 
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TIP RADIUS 



AXIAL DISTANCE Z 


FIG. 11. INTERPOLATION OF FLOW AT SLOT LOCATION 
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N04-M4-« 


Subroutine SLTFL0(L0P) 


Ob.iect 

Compute slot inlet flov 
Options 

L0P= Slot index 


List of Symbols 



ALST 

= a 

SLT 

3 

Slot injection swirl angle (rad) 

AMINF 

= Moo 

3 

Slot injection Mach number (dimensionless) 

AMUI 

- M-oo/l-lr 

3 

Freestream viscosity (dimensionless) 

AN 

= n 

5 

Boundary layer power law (l/7) 

ANUW 

= V 

3 

Wall kinematic viscosity (l/sec) 

A2CF 

= 2Cf 

3 

Wall friction coefficient (dimensionless) 

CCMH 

= 2.2/k 

3 

Constant for universal velocity (dimensionless) 

DEL 

= 6 

5 

Boundary layer thickness (ft) 

DELS 

* 

= <5 

3 

Boundary layer displacement thickness (ft) 

DELT 

= A 

3 

Boundary layer thickness (dimensionless) 

DSTR 

* 

= A 

3 

Boundary layer thickness (dimensionless) 

rWFL0 

= WsL 

3 

Weight flow in slot (dimensionless) 

DWFL01 

= WsL 

3 

Weight flow in slot (ib/sec) 

HSLT 

II 

3 

Slot height (dimensionless) 

H12 

= H 

12 

3 

Boundary layer form factor (6 /0) 

POINF 

“ ^00 

3 

Total pressure i^lot flow (dimensionless) 

RETH 

II 

w 

(D 

CD 

3 

Reynolds number, momentum thickness (dimensionless) 

RHUB 

= % 

3 

Hub radius (dimensionless) 

R0INF 

“ ^00 

3 

Freestream density (dimensionless) 
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Subroutine SLTFL0 (Cont'd) 


R0W 

RTIP 

SIGW 

TINF 

IW 

TOINF 

U 

UIN 

UINF 

USW 

WFL01 

WFL02 

Y 

YEAR 

YMID 

YPLUS 

YSL0T 


^■w 

= 

“ ®oc 
= 0w 

= 000 


= U 
= u 


- ^oo 

= u* 


o 



, Wall density (dimensionless) 

, Tip wall radius (dimensionless) 

5 Wall stress (dimensionless) 

, Freestream temperature (dimensionless) 

, Wall temperature (dimensionless) 

, Freestream total temperature (dimensionless) 

, Magnitude of velocity (dimensionless) 

, Magnitude of velocity (ft/sec) 

5 Freestream velocity (ft/sec) 

, Wall friction velocity (ft/sec) 

, Weight flow upstream of slot (dimensionless) 

5 Weight flow downstream of slot (dimensionless) 
, Distance from wall (dimensionless) 

, Distance from wall (dimensionless) 

, Midpoint of slot (dimensionless) 

, Universal distance (dimensionless) 

, Distance to slot lip (dimensionless) 


Theory 

The flow at the slot station is computed in two parts. First, the known 
upstream flow from the old Y grid is interpolated to the new Y grid from Rh ^^SL* 

In the second step, the flow entering through the slot is computed for YgL^YsR'i'. 
Since the total pressure and total temperature of the slot flow is given, and the 
static pressure inside the duct has been computed, the Mach number of the slot inlet 
flow can be computed. 
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Subroutine SLTFL0 ( Cont ' d ) 


M 


-= {-#r 


( 1 ) 


Then the remaining freestream properties are 




( 2 ) 


/ 7r« 


(3) 






The slot flow is assumed to be fully-developed, hence 


A = H/2 


(5) 


A^ = A/(l + 7) 


( 6 ) 


Coles’ velocity profile is assumed, and a friction velocity based on -Coles' profiles 
is computed using Function UC01ES. Then, the velocity is given by 


U = ~ /(lYl,Y^ 8) 
“ r 


(7) 


© 



( 8 ) 
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Subroutine SM00TH (X,Y,JX,XB,YB,JXB,JXK) 


Obj ect 

Least squares spline fit 
Options 


None 



List of Symbols 



A(I), B(I) 

= bj 

, Constants of integration 

CK(I,J) 


, Spline coefficients 

JX 


, Number of input points 

JXB 


, Number of output points 

JXK 


, Number of spline- nodes 

X(J), Y(J) 

= X.Y 

, Input curve 

X(I), Y(I) 

11 

, Output curve 

YPP(J) 

= Y" 

, Second derivative input curve 

XK(I), YK(I) 

II 

> 

, Node locations 

Theory 




The input curve is numerically differentiated to obtain the second derivative 
of Y(X). JXK equally spaced nodal points are assumed including the end points. 

The ISML routine returns spline coefficients of the form 

Yj"= Ci,3-D^ +Ci,2D2 + Cj_,D+y/' (1) 

where 

D=X, -X, (2) 
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Subroutine SM00TH (Cont'd) 


Equation (1) is integrated algebraically to obtain 


Y - 

‘ 20 


0 " + 


- 1.2 

12 


o'* + 


Cl. I 
6 


A. 

Y, 


D D + 


AjD 


+ Br + Y, 


(3) 


with the boundary condition 


that the end points of the curve remain unchanged . 
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Subroutine S0LVI 


Object 

Integrate equations of motion 
Options 


ALPHAK 


BETAH 


IDBG5 


IDBG8 


IDBG14 


= 1.0 Include convection terms 

= 0. Reyhner, Flugge-Lotz approximation 

= 0. Forward differences for convection 

= 1.0 Backward differences for convection 

= 0 No debugging printout 

= 1 Debug printout of viscous solution 

= 0 No debug printout 

= 1 Debug of matrix inversion algorithm 

= 0 No debug printout 

= 1 Debug printout of slot injection 


IFLAG 

IFSTEP 

IGDUM 

IVTEST 

I0PT11 


Used with IDBG8 

= 1 Test for optimum step size 

=2 Do not test for optimum step size 

= 0 Algebraic turbulence model when 10PT12 = 3 

> 0 Two equation turbulence model 

= 0 Do not print solution 

= 1 Print solution every station for separated flow 

= 0 OD wall boundary condition 

= 1 OD freestream boundary condition 
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Subroutine S0LVI (Cont'd) 


I0PT12 

I0TP13 

I0PT14 

I0PT15 

I0PT16 

I0PT17 

N0PT1 

N0TP2 

N0TP4 

N0PT5 


= 3 
= 0 
= 1 
= 0 
= 1 
> 0 
> 0 
= 0 
> 0 
= 0 
= 1 
= 2 
= 1 
= 2 
= 1 
= 2 
= 0 
> 0 


List of Symbols 


AA(I,J) = 

AB(I,J) = 

AC(I.J) = Cjj 


Algebraic turbulence models 
Two equation turbulence models 
No slots 
Slots 

No global iteration for separated flow 
Global iteration for separated flow 
Start calculation at station I0PT15 
End calculation at station I0PT16 
No restart 

Restart calculation at station I0PT17 

First station setup flow parameters 

First calculated station 

Continue calculation 

Adiabatic wall 

Nonadiabatic wall 

ID wall boundary condition 

ID axis of symmetry boundary condition 

No fatal errors 

Fatal error stop calculation 


, Element of A matrix 
, Element of B matrix 
, Element of C matrix 
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Subroutine S0LVI (Cont'd) 


List of Symbols (Cont'd) 


AD(I,J) 

ADI(I,J) 

AE(I,J,K) 

AQ(I) 

AZ(I,K) 

BA(I,K) 

BA(I,K) 

DDSTAR = A6* 

DDSTRM 

DEL0R 



= e 


IJ 



= z. 


, Element of D matrix 

, Element of 5“^ matrix 

rs 

, Element of E matrix 

, Element of Q vector 

, Element of Z vector 



Entropy increase when LL = 1 


Total temperature increase when LL = 2 




, Change in 
, Minimum change 
, Maximum step s 


* 

6 

in 

ize 


6 


* 


DSIG = AZ = (Z'^-Z'^ ^)/Z'^ ^ 
W W w 

DSIGM = (AZ)„ 

M 


Change in wall stress 
Maximum change in wall stress 


DSTAR 

EENTP 

EPRES 

ER0TH 

ER0US 

EUSUS 

EUSUSM 

EUPUP 


= 6 

= e( ) 

= E(n) 

= e( ) 

= u J 

= e(u2) 

= e(U2)^ 
s M 


= e(u2) 


FVJ(I,J,K) = F(I,J,K) 
FV(I,K) = FVJ(I,J+1,R) 
JJ 


v-1 

v-1 


, Displacement thickness 

, Truncation errors in 

, Truncation error in II 

, Truncation error in 

, Truncation error in ( U ) 

s 

, Truncation error in (lj2) 

s 

2 

f Maximum error in (U^ 

, Truncation error in (lj2) 

<)> 

, Dependent variables previous global iteration 
, Dependent variables backward differencing 
, Coarse grid station counter 
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Subroutine S0LVI (Cont'd) 


JKDS 

ENREF = 

PIREF = ii^ 

THREF = 

Theory 

^ol^u_ti^n_A^g^r^th^m_ 


, Fine grid station counter 
, Reference entropy for Mach transformation 
, Reference pressure for Mach transformation 
, Reference temperature Mach transformation 


The equations of motion are described in Ref. 1. Details of the finite difference 
equations and the solution algorithm using the method of block tridiagonal factoriza- 
tion is described in Ref. 2. The programmer may refer to Table I for the definition 
of rows and columns in the diagonal block. A description of the global iteration 
scheme for treating separated flow is given in Ref. 3. 

Truncation Errors 


The truncation errors calculated in this subroutine are those associated with 
linearizing the equations of motion by expansion in a Taylor series (see Ref. 2). 

Then for any two variables A and B we have 

b'-Va>' 

where A and B are mass flow weighted averages. 

_St^p_S^ze A 1 ^ 0 £i_thm 

The controlling factor in the selection of the streamwise step size is the 
thickness of the boundary layer rather than the length of the duct. Practical exper- 
ience has demonstrated that the streamwise step size should be of the order of mag- 
nitude of the boundary layer thickness. In some cases, for highly curved walls, the 
radius of curvature of the vcall may be a controlling factor. Since wall stress was 
found to be a sensitive indicator of wall curvature, a criteria based on the rate 
of change of wall stress was developed. In addition to these factors, the ADD code 
uses an analysis which linearizes the equations of motion using a Taylor series ex- 
pansion. Therefore, a third factor controlling the step size is this truncation 
error which should be kept below some maximum error. Finally, in certain cases, as 
equilibrium flow is approached, the boundary layer thickness grows slowly and the 
step size may be increased. 
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Subroutine S0LVI (Cont'd) 


The ADD code uses a fine grid interpolated from a coarse grid. Whereas the 
coarse grid is determined by the JL mesh points required to solve the potential 
flow through the duct, the fine grid (JLXKDS) is determined by the number of mesh 
points required by the boundary layer solution. At each grid point on the coarse 
mesh, the automatic step size algorithm calculates the fine grid parameter KDS 
using the criteria 


^J-H . sj* /AX \ 

KDS ■ ® I S* /m (2) 

where the J indices indicates the coarse grid step size. One integration step 
is taken and the following checks are made: 

AI < (AI)^ (3) 

€{Us^) < € (Us^)^, (A) 


If these conditions are satisfied, the integration continues for the KDS steps. 
If these conditions are not satisfied, KDS is doubled and the first integration 
step, is repeated. This process is continued until conditions (3) and (4) are 
satisfied. If KDS is doubled, IFSTL times the calculation is terminated. 

For cases approaching equilibrium flow 

If AS*< (AS*)^ , then (AX/S*)^, = 100 ( 5 ) 

and a much larger step size may be taken. 

Work-Loss Calculation 


Duct flow calculations require accurate determination of the work input and 
viscous losses. This subroutine integrates the following equations along a 
viscous flow streamline ^ = constant with the right hand side known. 


D (A©q) 
Dt 


= ( y - I) 


Vb 


H^/P 


( 6 ) 


D(AI) 

Dt 


n 


(I 


ns^ +/n/ } 


Re, + 


(7) 
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Subroutine S0LVI (Cont'd) 


where A6 q is the incremental work input produced by the blades, and AI is the 
incremental entropy rise due to viscous dissipation. 

References 


1. Anderson, 0. L.: Calculation of Internal Viscous Flows in Axisymmetric Ducts 

at Moderate to High Reynolds Numbers. Int. J. of Computers and Fluids, Vol. 8, 
No. A, pp. 391-411, December 1980. 

2. Anderson, 0. L.: User's Manual for a Finite Difference Calculation of Turbulent 

Swirling Flow in Axisymmetric Ducts with Slot Cooled Walls. Vols. I & II, 
USAAMRDL-TR-74-50, 1974. 

3. Anderson, 0. L. and D. E. Edwards: Extensions to an Analysis of Turbulent 

Swirling Compressible Flow in Axisymmetric Ducts. UTRC Report R81-914720-18 , 
1981. 
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TABLE I MATRIX ELEMENTS AA(IVAR.IEQ) 


Col. 


Row 


IVAR 

Variable 

lEQ 

Equation 

1 


1 

Heat flux equation 

2 

U 

s 

2 

Equation state 

3 


3 

S momentum equation 

4 

n 

4 

i(i momentum equation 

5 

I 

5 

Energy equation 

6 

© 

6 

Continuity equation 

7 

P 

7 

I stress equation 

ns 

8 

i 

ns 

8 

1 stress equation 

n(f) 

9 

1 

H({) 

9 

N momentum equation 

10 

Q 

10 

Entropy equation 
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Subroutine STRESI (L0P) 


Ob j ec t 

Calculate initial stress and heat flux 


Options 




L0P 

= 1 

Calculate F(L,2,K) 

L = 8, 10 


= 2 

Calculate F(L,1,K) 

1 = 8, 10 

I0PT17 

= 0 

Initialize stress at 

; first station 


> 0 

Calculate turbulent 

viscosity only 

List of 

Symbols 




See C0MM0N BL0CK variables 
Theory 

With the velocity distribution, turbulent viscosity, and turbulent conductivity 
known, the stress and heat flux can be calculated using the algebraic turbulence 
model described in Ref. 1. In addition, the total temperature and total pressure 
at the initial station are calculated for use by the work and loss calculation in 
subroutine S0LVI . 

References 


1. Anderson, 0. L. and D. E. Edwards: Extensions to an Analysis of Turbulent 

Swirling Compressible Flow in Axisymmetric Ducts. UTRC Report R81-91A720-18 , 
February 1981. 
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Subroutine SUBLAY (AK, AP, Y, Q, EPS) 


Object 

Initializes k, t turbulence model in inner layer 
Options 

None 

List of Symbols 
AK 
AP 
Y 
Q 

EPS 
AMUT 
Theory 

This subroutine initializes k and c in the inner region of the boundary layer 
by determining first the turbulent viscosity through the van Driest wall function and 
then determining the kinetic energy and dissipation through the relations described 
in Reference 1. 

References 

Anderson, 0. L. and D. E. Edwards, Extension to an Analysis of Turbulent 
Swirling Compressible Flow in Axisymmetric Ducts. NASA Contract NAS3-21853, 
1981, UTRC Report T81-914720. 


, von Karman constant 


, von Driest constant 


.+ 


K 


, Universal coordinate (p u /y ) 

w w 

2 

, Universal kinetic energy (k/u* ) 
, Universal dissipation (p/pc/u*^) 


, Universal turbulent viscosity (y_/y ) 

T w 
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Subroutine TPRINT (JNUM) 


Ob.i ec t 

Print CPU time 
Options 

ITIME = 0 Initialize time 

# 0 Continue time 

List of Symbols 

JNUM , Not used 

Theory 

This is a dummy subroutine which calls UNIVAC time routine SUPS. 
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Subroutine TRBLDE (X.Y.RP.ZP) 


Object 

Transformation from stacking plane to duct plane 
Options 

None 

List of Symbols 


RCLO, ZCLO = R^, 

CLO CLO 

RP, ZP = R, Z 


THCL 

®CL 

X,Y 

= X,Y 


, Origin of stacking plane (dimensionless) 

, Coordinates in duct plane (dimensionless) 

, Rotation of stacking plane (deg) 

, Coordinates in stacking plane (dimensionless) 


Theory 

The transform from the stacking plane to the duct plane is given by 



+ Y 

cos6„ 

+ X 

sin6„ 

(1) 

CLO 


CL 


CL 


z = 

+ X 

CD 

CO 

o 

CJ 

- Y 

sine 

(2) 

CLO 


CL 


CL 
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Subroutine TRLETE (KK,XLE,YLE,XTE,YTE) 


Ob.i ec t 


Find leading and 

trailing edge in stacking plane 

Options 

None 


List of 1 

Symbols 


ALP 

a 

, Stagger angle (deg) 

C0RD 

C 

, Chord (dimensionless) 

XCL, YCL 

= ^cl’\l 

, Coordinates of centerline (dimensionless) 

XLE, YLE 

= \e’\e 

, Coordinates of leading edge (dimensionless) 

XTE, YTE 

ss Y Y 

TE’ TE 

, Coordinates of trailing edge (dimensionless) 

Theory 

The 

projection of 

the blade to the (r, Z) plane is 


= X^^ - 1/2 C 

cos (a) 

LE 

CL 

Y • 

= Y 


LE 

CL 



= X + 1/2 C 

cos (a) 

TE 

CL 

Y 

= Y 


TE 

CL 



( 1 ) 

( 2 ) 

(3) 

(4) 
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Subroutine TURB 


Object 

Compute turbulent viscosity using algebraic turbulence model. 
Options 


I0PT12 = 0 
= 1 

= 2 


Original eddy viscosity model 

Eddy viscosity model with low Reynolds 

number corrections (not tested) 

Eddy viscosity model with curvature corrections 

List of Symbols 



AMUE 

= p.,/y 
E r 

, Effective viscosity (dimensionless) 

AMUER(K) 

= (K > 

Too r 

jFreestream turbulent viscosity (dimensionless) 

AMUM 

= y /y 

,Freestream molecular viscosity (dimensionless) 

AMUW 

= y /y 
w r 

jWall value of molecular viscosity (dimensionless) 

AMUWK 


.Molecular viscosity (dimensionless) 

AMUO 

U 

II 

.Maximum freestream viscosity (dimensionless) 

BETAC 

= 

.Eide & Johnston coefficient (6.) 

BETAD 


.Eide & Johnston coefficient (6.) 

BETAPB 

= 

.Bradshaw coefficient (6.) 

BETAS B 

= “s 

.Bradshaw coefficient (0.) 

DELO 

= A* 

.Displacement thickness (dimensionless) 

DU 

= An 

.Velocity finite-difference (dimensionless) 

DUDZ 

= dU/dZ 

.Velocity derivative (dimensionless) 

E 

= E 

.Rate of strain (dimensionless) 

EM 

= E„ 

.Strain freestream (dimensionless) 
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Subroutine TURB (Cont'd) 


Variables (Cont ' d ) 


EMH,EM^ 



ENP 

= 

^n(f. 

ENS 

= 

^ns 

EW 

= 

Ew 

FCURV 

= 

F 

PHI 

= 


RHOM 

= 


RC 

= 

^ic 

RICPB 

= 


RICSB 

= 


RP 



SIGWH 

= 

r 

^wh 

SIGWK 

= 

Ewk 

SIGWT 

= 

^wt 

TPLUSl, 

TPLUS2 

= 

T+ t‘ 
1 : 

UK(K) 

= 


UM 

= 

u 

K> 

USTARH, 

USTART 

- 

U*,U^ 
H ’ 

Y 

= 

Y 

YK 

= 

Y 

k-fis 


, Strain hub, tip, edge of inner layer 
(dimensionless) 

, Swirl rate of strain (dimensionless) 
,Streamwise rate of strain (dimensionless) 

,Wall rate of strain (dimensionless) 

Curvature correction to eddy viscosity 
, Integrand (dimensionless) 

, Density freestream (dimensionless) 

,Streamwise component of Richardson Number 
, Tangential component of bulk Richardson number 
,Streamwise component of bulk Richardson number 
, Tangential component of Richardson number 
,Wall stress hub (dimensionless) 
jWall stress inner layer (dimensionless) 

,Wall stress tip (dimensionless) 

, Universal wall stress (dimensionless) (x/x^) 

, Magnitude of velocity (dimensionless) 

, Freestream velocity (dimensionless) 

, Friction velocity hub, tip (dimensionless) 

, Distance across duct (dimensionless) 

, Distance across duct (dimensionless) 
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Subroutine TURB (Cont*d) 


Variables (Cont'd) 


YMHjYMT = Y ,Y .Distance to inner layer hub, tip 

(dimensionless) 

YPLUSl, = Y'*^,Y''’ .Universal distance (dimensionless) 

1 7 

YPLUS2 

Theory 

The algebraic turbulence models are described in Ref. 1. 


References 


1. Anderson, O.L. and D.E. Edwards, Extension to an Analysis of Turbulent 
Swirling Compressible Flow in Axisymmetric Ducts, NASA Contract No. 
NAS3-21853, 1981, UTRC Report T81-914720. 
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Subroutine TURB2Q 


Object 

Calculate turbulent viscosity using k, E turbulence model. 
Options 


NTURB2 = 0 


Set up initial conditions 


^ 0 


Calculate 2-equation turbulence model 


I0PT17 


Restart capabilities, IOPT17>0 calculation 

List of Symbols 




DHF(l,r) 


,Not used 


DHF(2,K) 

= 

.Effective Viscosity 


DPF(1,K) 


,Not used 


DPF(2,K) 

~ (wp/y_) 

Prt 

.Effective thermal conductivity 


TRB(1,1,K) 


2 

.Turbulence kinetic energy, (k/Uj."^) 


TRB(1,2,K) 

< 

.Turbulence kinetic energy (k/uj."^) 


TRB(1,2,K) 


.Turbulence kinetic energy (k/u^^) 


WRELAX 

0) 

.Relaxation parameter (.5) 


ITMAX 


.Maximum number of iterations (100) 


TRB(2,1,K) 

eJ-1 

K 

.Turbulence dissipation (Pr/yr ^ 

) 

TRB(2,2,K) 

E'^ 

K 

.Turbulence dissipation (pj./yj. ^ 

) 

TRB(2,3,K) 

v-1 

Tr 2 

.Turbulence dissipation (pj-/yj. ^ 

) 

AMUL(K) 


.Molecular viscosity 


AMUT(K) 

j-i/ 
y /y 

K 

.Turbulent viscosity 


AMUE(K) 

, J-1 

(y£/yr)K 

.Effective viscosity 



11-236 



Subroutine TURB2Q (Cont'd) 


List of Symbols (Cont'd) 


E(K) ,FM(K) 



, Matrix solver 

TP(K) 

= 

pJ 

K 

.Turbulence production (dimensionless) 

TRE(K) 


pJ-1 

^ek 

.Turbulence Reynolds number 

TRJ (K) 

= 

J-1 

\tK 

.Richardson number 

CDAMP 

= 


. .0115 

Cl 

= 


.1.35 

SIGE 

= 

*^E 

,1.3 

SIGK 

= 

a-K 

.1.0 

CMUO 

= 

C, 

1 ' 

..09 

CMU 

= 


,.09* i_e(-CDAMP*Y^) 

C2 

= 

^2 

.1.8* 1.- ^ * e -B 2/36 (i _.2 r. ) 
1.8 ^ 

C3 

= 

s 

e-5Y-^ 

YPIUS 

= 

Y+ 

.Universal distance 

SK(1,K) 

= 



SB(2,K) 

= 

^Be 


SD(1,K) 

= 

\ 


SD(2,K) 


^^E 


GAM(L,K) 




GAM(2,K) 

= 

^E 

,M+h^/0g 

TURBIF 

= 

u'u'/u2 

.Freestream turbulence (not tested) 
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Subroutine TRUB2Q (Cont’d) 


List of Symbols (Cont'd) 
QPLUS 

EPS E+ 

ERRREY 

ERRREM 

AKINF K 

00 

EINF E„ 

UINF U 

•O 


, Universal turbulent kinetic energy u'/(u*) 

.Universal turbulent dissipation 

.Relative change in R per iteration 

.Error tolerance (.001) 

.Freestream turbulence kinetic energy 
(not tested) u'u'/u^ 

.Freestream turbulence dissipation (not tested) 
.Freestream velocity 


Theory 

The two-equation turbulence model of Chien is solved using the method 
described in Ref, (1). 

References 


1. Anderson, O.L. and D.E. Edwards, Extension to an Analysis of Turbulent 
Swirling Compressible Flow in Axisymmetric Duc'ts, NASA Contract No. 
NAS3-21853, 1981, UTRC Report T81-914720. 
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Function UBLAS(Y) 


Ob.i ec t 

Calculate Blasius velocity profile 
Options 
None 

List of Symbols 

UBLAS = u/u , Velocity ratio 

00 

Y = y/ 6 , Distance from wall 

Theory 

The Blasius velocity profile is linearly interpolated from a table of values 
given in Ref . 2 . 

Reference 

1. Schlicting, H. : Boundary Layer Theory. 6th Ed., McGraw Hill, N.Y., 1968. 
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Subroutine UMACH (PT, TT, P, U) 

Obj ect 

Calculate inviscid velocity 
Options 

None 

List of Symbols 


p 

P/P 

r 

, Static pressure (dimensionless) 

PT 

Pi/p, 

, Total pressure (dimensionless) 

TT 


, Total temperature (dimensionless) 

U 

u/u 

r 

, Velocity (dimensionless) 


Theory 

The inviscid velocity U is calculated from the input P, PT, TT using the 
isentropic flow relations. 
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Subroutine VARFUN (JINDEX, INDEX, MINDEX) 


Obj ec t 

Store variables to be used by PRT0UT 
Options 
None 

Input Variables 


JINDEX 


Index of streamwise station 

INDEX 


Index of viscous flow quantities 

MINDEX 


Index of effective viscosity 

F 


Viscous flow array 

Q 


Geometry flow array 

KL 


Number of streamwise stations 

Output Variables 



VARBL(1,J,K) 

4 > 

, Stream function (dimensionless) 

VARBL(2,J,K) 

Us 

, Streamwise velocity (dimensionless) 

VARBL(3,J,K) 


Swirl velocity (dimensionless) 

VARBL(4,J,K) 

P 

, Density (dimensionless) 

VARBL(5,J,K) 

^ns 

, Streamwise stress (dimensionless) 

VARBL(6,J,K) 

p j* 

. Effective viscosity 

VARBL(7,J,K) 

R 

, Radius (dimensionless) 

VARBL(8,J,K) 

Z 

> Axial distance (dimensionless) 

VARBL(9,J,K) 

V 

, Metric scale coefficient (dimensionless) 
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Subroutine VARFUN (JINDEX, INDEX, MINDEX) (Cont ’d) 


Output Variables (Cont'd) 


VARBL(10,J,K) G 
VARBL(11,J,K) 

oS 

VARBL(12,J,K) 

dN 

VARBL(13,J,K) 

9S 

VARBL(14,J,K) y/h 


j Gap (dimensionless) 

Streamwise derivative of radius 

’ Curvature of potential line 

, Curvature of stream line 

Fractional distance across duct 


Theory 

Module PRT0UT requires information from three streamwise stations in order 
to perform the differencing of differential equations. This subroutine stores 
the information required in array VARBL. 
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Subroutine WBLEED 


Object 

Calculate perforated wall bleed 
Options 

10PT18 =0 No wall bleed 
= 1 Tip wall bleed 

= 2 Hub wall bleed 

= 3 Tip/hub wall bleed 

PCHEK > 1.0 Flow enters tunnel 

< 1.0 Flow leaves tunnel 

Input Variables 


AHAS 

- 

Ratio of hole area to surface area 

CDISH 

= c 

Discharge coefficient 

PTP 

= p 

TP 

Plenum total pressure 

TTP 

= 

T 

Plenum total temperature 

Internal Variables 


AMTU 

II 

Tunnel Mach number 

GAMMA 

= y 

Ratio of specific heats 

GASR 

= R 

Gas constant 

PS 

= p 

Static pressure (psfa) 

PT 

^T 

Total pressure (psfa) 

PSTU 

" ^TU 

Tunnel static pressure (psfa) 

PTTU 

= P 

TTU 

Tunnel total pressure (psfa) 

TSTU 

= T 

TU 

Tunnel static temperature (deg R) 

TTTU 

= T 

TTU 

Tunnel total temperature (deg R) 

RH0R 

= P 

r 

Reference density (slug/ft ) 

USR 

= u 

r 

Reference velocity (ft/sec) 

PRESR 

= P 

r 

Reference pressure (psfa) 


1I-2A3 


Subroutine WBLEED (Cont'd) 


Internal Variables (Cont'd) 


TEMPR = T 

r 

SGN i 1 
Output Variables 

RH(9,J) = (pVh 

RT(9,J) = (PE )„ 
n 1 

Theory 


Reference temperature (deg R) 

Sign convention 

2 

Mass bleed hub wall (slugs/ft /sec) 

2 

Mass bleed tip wall (slugs/ft /sec) 


If one treats a single hole in a perforated wall as an orifice, then the mass 
flow can be derived in terms of the plenum stagnation conditions and the local static 
pressure inside the tunnel Holman (Ref. 1). Then an expression for the mass flow 
added to the tunnel flow is given by 




As yr'itJj 



( 1 ) 


where P.p and Tj are the plenum conditions, P is the local tunnel static pressure, 
A^/Ag is the ratio of the hole area to surface area, and C the effective discharge 
coefficient which is a property of the perforated wall. If the tunnel static pres- 
sure is greater than the plenum total pressure, the mass flow bleed is out of the 
tunnel. Under these conditions, P'p and T.j. are taken from the wind tunnel condi- 
tions, and P is the plenum pressure which is assumed known. 

The mass flow bleed is related to the stream function by 

d'f G . 

a.S "V p^'j^ 

Equations (1) and (2) provide the boundary condition for a perforated wall relating 
two dependent variables jp and P in terms of the characteristics of the perforated 
wall and the plenum conditions. 

The program checks the options according to the table below. 
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Subroutine WBLEED (Cont'd) 


O.D. Wall 


PCHEK = > 1.0 set 

TP STU 


P = P 
T TP 

P = P 

STU 

T = T 
T TP 


PCHEK = P „/P 


TP STU 


< 1.0 set P„ = P 


STU 


TP 

T = T 
T STU 


SGN = -1.0 


SGN = 1.0 


Reference 

1. Holman, J. P.: Experimental Methods for Engineers. McGraw-Hill 

New York. 1966. 


I.D. Wall 


+ 1.0 


- 1.0 


Book Co . , 
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Subroutine WFITER 


Object 

Inlet weight flow iteration 


Options 

WFLI = 0 
> 0 

I0PT1 =3,7 
= 4,8 


No weight flow iteration 
Weight flow iteration 

Iterate on Mach number 
Iterate on ID static pressure 


List of 

Symbols 



AMSI 

= 

V 

M 

Mach number 

BINPUT(3 

,1,1) = 

V 

% 

ID static pressure 

B311 

= 

% 

Input ID static pressure (P„/P ) 

n r 

EPRPIN 

= 

Ep 

Error in normal pressure gradient 

ERRW 

= 


Error in weight flow 

IWFLl 

= 

V 

Iteration counter 

N0PT19 

= 

V 

Iteration counter 

WFLI 

= 

w 

Input weight flow (Ib/sec) 

WFL0 

= 

V 

w 

Output weight flow (Ib/sec) 

Theory 





Subroutine FL0WIN calculates an output weight flow WFL0. This subroutine uses 
an iteration based on Newton's method to adjust the weight flow keeping the tPtal 
pressure and total temperature fixed. 

For I0PT1 = 3 or 7, the Mach number is varied using 


For 10PT1 


- M*^ w/w*^ 


4 or 8, the ID static pressure is varied 

'dn, 


rr rr *' /^^h\ / 


( 1 ) 


( 2 ) 
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Subroutine WFITER (Cont'd) 


The error in weight flow is given by 


= |W-W*'|/w 


( 3 ) 


The error in normal pressure gradient is obtained from subroutine ERPIN. 
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Subroutine WRITEPF(JJ) 


Object 

Store updated potential flow solution. 
Options 

None 


List of Symbols 
JJ 
NST 
NIST 

NPDRK =24 
P 


,JJth station in core 
,No. records per block 
,No. words per block 
,Unit number 
.Stream function 


Theory 

The stream function array P(JJ,KK) is arranged in core as described in 
subroutine READPF, Fig. 12. When an iterative sweep of one block is 
complete, the new updated solution is written on a disk file. This occurs 
when JJ = NST-1. 
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FIG. 12 I/O for Streamline Curvature Solution 
Overlapping Blocks 


Record No. 
J 

1 

2 

3 

4 

23 

24 

25 

26 

46 

47 

48 


70 

71 


93 

94 

95 

96 

97 

98 

99 
100 
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Subroutine WRTBLD 


Object 

Write output for WBLEED 
Options 

None 

Variables 

See subroutine WBLEED 
Theory 

This subroutine prints the output from subroutine WBLEED according to the 
output format described in Section 4.3. 
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Subroutine WRTCKI 


Obj ect 

Write output for CKINPT 
Options 

None 

Variables 

See subroutine CKINPT 
Theory 

This subroutine prints the output from subroutine CKINPT according to the 
output format described in Section 4.3. 
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Subroutine WRTCAL 


Obj ect 

Write output for CALINV 
Options 

None 

Variables 


See subroutine CALINV 
Theory 


This subroutine writes the output from subroutine CALINV according to the 
output format described in Section 4.3. 
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Subroutine WRTF0U 


Object 

Print solution on output data line 
Options 

None 

Variables 


C0MM0N BL0CK Variables 


Theory 


"i3 = 


at the output data points are printed in 
direction cosines given by 

cose 

-SIN0 

0 


SIN6 

cose 

0 


0 

0 

1 



( 1 ) 


where 6 is the angle between the S coordinate and the z axis (R,Z) coordinates, 
or the angle between the S coordinate and the X axis (Y,X) coordinates. 
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Subroutine WRTGDC 


Ob.i ect 

Write output summary - coordinate calculation 
Options 

None 

List of Symbols 

C0MM0N BL^CK variables 
Theory 

The output from GDUCT is printed on page 1 and the output from C00RST is 
printed on page 2. 
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Subroutine WRTINP 


Object 

Write input parameters and variables 
Options 

The options are the same as the input read options 
List of Symbols 

Common block variables 
Theory 

The input data is printed as follows: 


Page 

1 

Input 

parameters 

Page 

2 

Inlet 

flow data 

Page 

3 

Strut 

data 

Page 

A 

Strut 

thickness distribution 

Page 

5 

Strut 

flow variables 
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Subroutine WRTS0V 


0b.1 ec t 

Print solution at selected stations 
Options 

I0PT4 Print solution every I0PT4th station 

IVTEST > For separated flow print every station 

I0PT4 < Extended printout every I0TP4th station 

List of Symbols 

C0MM0N BL0CK variables 
Theory 

The dependent variables are printed according to the input print option I0TP4. 
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Subroutine WRTSUM 


Object 

Write output summary pages 
Options 

I0PT15 First calculated station 

I0PT16 Last calculated station 

JST0 Number of stations stored on F file 

I0PT15 1 J 1 10PT16 Solution complete 

I0PT15 i J - JST0 + I0PT15-1 Solution not complete 
Variables 

C0MM0N BL0CK variables 
Theory 

This subroutine prints the mass flow weighted average flow, pressure and skin 
friction coefficient, and convective heat transfer data. 
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Function XH(J) 


Object 

Find distance on hub wall 
Options 

Evaluate at point J 
List of Symbols 
None 
Theory 

The distance along the hub wall is given by 

XH = 


R^(J) - R^(J-I)^ + AZ^ 


1/2 


Q') 
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Function XT(J) 


Object 


Find distance along tip wall 
Options 

Evaluate at point J 
List of Symbols 
None 
Theory 

The distance along the tip wall is given by 


= |[r^{j)-r^(j-i)]V [zT(j)-Z-r(J-i)]' 


MZ 


XT 


(2) 


Subroutine ZER0 (ADDR, NUM) 


Object 

Set real variables to zero 
Options 


None 

List of Symbols 


ADDR Real array 

NUM Number of elements 

Theory 


This subroutine sets the first NUM real variables in array ADDR to zero. 


11-260 



►^1 


8.0 LIST OF SYMBOLS 




B 


C 


f 


C 


f 


C 

P 




e 

E 

E ,E ^ 
ns nij) 

F 


g 

G 

h 

H 

I 

k 

K 


Coefficients used in matrix solver 

Van Driest constant (26.0) 

Constant used in coordinate stretching 

Convection terms 

Skin friction coefficient 
2 2 

Specific heat (ft /sec /*R) 

Empirical constants in Two Eq. 

Turbulence model 

Rate of strain (units) (1/sec) 

P r 2 

Dissipation Rate __r r 
dimensionless u^ 

Streamwise and tangential strain 
(dimensionless) (e/(p u ^)) 

Ratio mixing length to flat plate mixing 
length 

Ratio mixing length to wake mixing length 
Coefficient in Poisson equation 
Gap between walls (ft) 

Gap between walls (dimensionless) (g/r^) 
Duct height (ft) 

Height of duct (dimensionless) (h/r^) 

2 

Turbulence Intensity (u'u’/u ) 

2 

Turbulent kinetic energy (ft/sec) 

Turbulent kinetic energy (dimensionless) 

(k/u^^ 
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LIST OF SYMBOLS (Cont'd) 


K 

c 


n 


K 

s 


K, 


K 


I 

I 

o 

M 

M 

r 

n 


P 

P 

P 


r 

T 


P 

r 



q 


Q 


r 

n 


r 

r 


Streamline curvature correction 
(dimensionless) 

Curvature of coordinate potential line 

(dimensionless) (r /r ) 
r n 

Curvature of coordinate streamline 

(dimensionless) (r /r ) 
r s 

Curvature of the swirl flow (dimensionless) 

Freestream turbulent kinetic energy 

(dimensionless) (u'u' /u ^) 

•0 r 

Mixing length (ft) 

Mixing length flat plate (ft) 

Mach number 


Reference Mach number 

Normal coordinate (dimensionless) 

Static pressure (psf abs) 

Reference pressure (psf abs) 
Turbulence production (dimensionless) 



Prandtl number (laminar) (C y/X) 

P 

Prandtl number (turbulent) (C^y/X)^ 

Heat flux (Ib/ft/scc) 

Heat flux (dimensionless) q/(p u C T ) 

r r p r 

Curvature of coordinate potential line 
Reference radius (ft) 
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r 

s 

r 


R 


R. 

IS 


is 


it 


1 * 


R. 


R 

s 



s 


S 


S 


s 


s 


E 

k 

K 


max 


s n 


u 

<P 


Curvature of coordinate streamline 

Curvature of the swirl flow 

Radius, r/r^ (dimensionless) 

Richardson numbers for streamwise curvature 
(dimensionless) 

Turbulent Richardson number (dimensionless) 

Richardson numbers for swirl (dimensionless) 

Radius of curvature of potential line 
(dimensionless) 

Radius of curvature of streamline (dimen- 
sionless.) 

Turbulent Reynolds number 

Reference Reynolds number (p u T /p ) 

r r r r 


2 2 

Gas constant (ft /sec /*R) 


Streamwise coordinate (dimensionless) 


Dimensionless source in dissipation equation 

(rr^/(u^Ur^)(S^) 

Source Term in turbulenct kinetic energy 
equation (units) 


Dimensionless source is turbulent kinetic 

3 

energy Eq. \) 

Source term in dissipation equation (units) 
Magnitude of freestream velocity (ft/sec) 
Reference velocity (ft/sec) 

Velocity components (ft/sec) 
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u ,u ,u 

s n d 


a ,a^ 
s 4> 


s 4> 


k e 


Friction velocity ft/sec) 

2 2 

Reynolds stress (ft /sec ) 

Magnitude of velocity (dimensionless) 
(u/u^) 

Velocity components (dimensionless) (u/u^) 

Freestream velocity components (dimen- 
sionless) (u/u^) 

Universal velocity (u/u*) 

Metric scale coefficient (dimensionless) 

Distance normal to wall (ft) 

Distance normal to wall (dimensionless) 

(y/r^) 

Universal distance to wall (p u y/p ) 

w w 

Axial distance (ft) 

Axial distance (dimensionless) (z/r^) 

Swirl angle (deg) 

Empirical constants for Bradshaw mixing 
length correction 

Empirical constants for Eide-Johnston 
mixing length correction 

Ratio of specific heats (dimensionless) 

Coefficients in two eq. turbulence model 

Displacement thickness (ft) 

Curvature of viscous flow streamlines 
relative to coordinates 
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e 

e 

0 

©, 

n 

I 

K 

X 


y 


n 

p 

p. 


w 


°E’°K. 


I ^ 
ns n()> 


Turbulence dissipation (ft /sec ) 

Angle of streamline to axis (deg) 

Static temperature (T/T^) 

Total temperature (T^/T^) 

Transform normal coordinate- 
Entropy dimensionless (I/R) 
von Karmon constant (.41) 

Thermal conductivity (Ib/sec/deg R) 
Molecular viscosity (slug/ft/sec) 
Effective viscosity (slug/ft/sec) 
Reference viscosity (slug/ft/sec) 
Turbulent viscosity (slug/f t/sec) 
Iteration counter 

Static pressure (dimensionless) (p/p^) 

Total pressure (dimensionless) 

3 

Density (slug/ft ) 

3 

Reference deposits (slug/ft ) 

3 

Density, at wall (slug/ft ) 

Density (dimensionless (p/p^) 

Empirical constants in two eq. turbulence 
model 

Streamwise and tangential stress (dimen- 
sionless) (t/p^/u^^) 
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T 

X 

X 

V 

10 

(0 


10 ,(0 

n s 





2 

Shear stress (Ib/ft ) 

Clauser Constant (.0168) 

Normal coordinate transform 3n/Sn 

2 

Stream function (ft /sec) 

Stream function (dimensionless) 
Relaxation factor 
Wilcox dissipation rate 
Vorticity components (dimensionless) 
Relaxation factor 
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